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PERFORMANCE OF FINNED THERMAL CAPACITORS 


By 

William R. Humphries 
SUMMARY 


The objective of this study was to investigate the performance of 
typical thermal capacitors, both in earth and orbital environments. Techniques 
which were used to make predictions of thermal behavior in a one-g earth 
environment are outlined. Orbital performance parameters are qualitatively 
discussed, and those effects expected to be important under zero-g conditions 
are outlined. A summary of thermal capacitor applications are documented, 
along with significant problem areas and current configurations. An experi- 
mental program was conducted to determine typical one-g performance, and 
the physical significance of these data is discussed in detail. Finally, numer- 
ical techniques were employed to allow comparison between analytical and 
experimental data. 

A transparent test specimen was built for the experimental study to 
allow visual observation of the phase change front. The test specimen 


xxii 


consisted of two compartments, a phase change material (PCM) housing and 
a finned fluid passage. The PCM housit^ contained aluminum fins, which 
partitioned the PCM compartment into 1.9-cm (3/4-in.), 1.27 -cm (l/2-in. ), 
and 0.635-cm (l/4-in. ) cells. During all tests, nonadecane normal paraffin 
material was used as the PCM. By introducing Monsanto Coolanol-15 into the 
finned fluid passage, the PCM was either cooled or heated. Heating and 
cooling rates were varied by varying the coolanol inlet temperature from 
1103.4 j/m^-sec (350 BTU/ft^-hr) to 204.9 j/m^-sec (65 BTU/ft^-hr), and 
from 567.5 J/m^-scc (180 BTU/ft^-hr) to 126.1 j/m^-see (40 BTU/ft^-hr), 
respectively. 

Instrumentation consisted of temperature sensors located in the paraf- 
fin, on the fins, in the fluid, and on the plate separating the PCM from the 
coolanol compartment. In addition, low speed 16 -mm motion pictures were 
used to document freezing and melting rates. 

Using an explicit forward finite differencing technique, analytical models 
of the 1.9-cm (3/4-in.) and 0.635-cm (l/4-in, ) cells were developed. Experi- 
mental fin and plate temperatures were used as boundary values for analytical 
models and these models were used to generate PCM temperature profiles and 
interface positions. 

The melt model, which utilized a Rayleigh number correlation to 
incorporate convection, satisfactorily predicted the melt front position. 
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However, a pure conduction model used to predict the freeze front position, 


lagged the expected positions. It was postulated that this discrepancy between 
the analytical model and the test data was due to the model using a planar 
interface rather than the actually observed irregular interface. The irregular 
interface was a result of numerous dendrites which formed at the liquid/solid 
boundary. 

Zero-g investigations indicated that the most significant difference be- 
tween one-g and zero-g thermal performance was the absence of buoyancy 
driven convection at zero-g. However, appreciable convection can be induced 
by surface tension forces. The ullage bubble, if present, could inhibit heat 
transfer at zero-g. However, this is very unlikely when paraffin is the PCM, 
due to its good wetting characteristics. 
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CHAPTER I 


INTRODUCTION 

Definition of a Thermal Capacitor 

From an electrical analogy, any device which has the capacity to store 
quantities of thermal energy can be defined as a thermal capacitor (TC) . How- 
ever, for the purpose of this study, only a device that absorbs and rejects 
energy by using a material which undergoes a phase change will be considered 
a thermal capacitor. This definition is further restricted in this study to 
include only latent heat effects associated with the melting process and with 
solid-state structural changes of phase. Although these devices will be 
referred to as thermal capacitors, other names such as thermal flywheel, 
phase change device, and fusible mass device have been used to describe the 
same equipment concept. 

Fundamentally there arc three uses to which the TC has been applied: 
(a) Thermal damping of oscillatory outputs, (B) inhibiting thermal excusions, 
and (c) maintaining constant temperature. Thermal responses in these situ- 
ations arc illustrated in Figure 1. 

Typically a TC is a passive device, Having no mo\ ing parts, which is 
normally composed of three integral components: an external housing, an 
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EMPERATURE 


A. THERMAL DAMPING OF OSCILLATORY OUTPUTS 


TIME 



B. INHIBITING THERMAL EXCURSIONS 



W/OTC 
WITH TC 


C. MAINTAINING CONSTANT TEMPERATURE 



W/OTC 
WITH TC 


Figure 1. Responses of a Thermal Capacitor 


internally finned filler material, and the phase change material (PCM) . The 
device being thermally controlled is excluded even though the capacitor may be 
an inseparable part of this device (Fig. 2). The PCM is a material which 
experiences transition or phase change in the regime of operating temperatures 
of the device (for the purpose of this study, only liquid/solid and solid/solid 
transitions were considered) . 



Figure 2. A thermal capacitor device. 

Statement of Problem 

Since the advent of manned space flight, aerospace designers have been 
looking for better methods of thermally controlling flight vehicle systems. In 
certain applications, thermal capacitors are emergii^' as strong competitors. 
The ability of thermal capacitors for either heatii^ or cooling, or for storing 
thermal energy, indicates their thermal flexibility. The fact that they are pas- 
sive makes them quite attractive, when compared to less reliable active 
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systems; such as, gas blowers, sublimators, and liquid pumped systems, to 
name a few. Furthermore for low and intermediate heat storage applications, 
the capacitor can represent a weight savings over other thermal control 
techniques. In addition to being an independent thermal control device, thermal 
capacitors can be used to supplement other devices such as, heat pipes and 
space radiators. Thermal capacitors have recently been used on an earth sat- 
ellite, three moon vehicles, and Skylab. Possible future design applications 
include using capacitors on a manned space station and in future space shuttle 
flights. 

Most designers have employed pure conduction models for design anal- 
yses, without considering convection during melting, the multidimensional 
effects, and the dendrite formations at the interface during freezing. Asa re- 
sult, significant discrepancies have existed in some cases between analytical 
models and one-g test data. When this occurs, the test data repeatedly indi- 
cate enhanced thermal performance, and the designer will normally use the 
more conservative analytical data as the design criteria. Unfortunately, this 
conservatism often leads to an overdesigned capacitor. 

Study Approach 

The objective of this work has been to investigate typical thermal 
capacitor performance in both earth and orbital environments. The motive was 
to gain insight into the processes involved in the thermal performance of a 
typical space vehicle phase change device. The study was directed toward 
accumulating experience which may assist future capacitor designers. 
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A survey of previous PCM studies and sources for both paraffinic and 
nonparaffinic materials property data has been included in this study. Tech- 
niques for estimating one-g performance are discussed, and quantitive compar- 
ison of these techniques, with possible zero-g influences, is presented. To 
provide additional insight into the different design techniques, a survey of past, 
current, and possible future applications is given. The significant problems 
encountered in these applications, with which the author is familiar, are out- 
lined. In addition, physical desciptions of capacitors used in space vehicle 
applications are profiled. 

An experimental program was conducted to study the phase charge 
processes occurring in a t>"pical capacitor. The test item was fabricated from 
transparent material to allow visual studies of the internal processes that were 
occurring within the PCM housing. Using thin aluminum fins, the internal 
PCM housing was divided into 0.635-cm (l/4-in.), 1.27 -cm (l/2-in. ) and 1.9- 
cm (3/4-in.) cells. Test data included thermocouple temperature gage out- 
puts and motion picture film data. Temperature data of the fin and lower heat 
transfer surface were input into single-cell computer models to determine the 
phase change front position time history. The film data also was used to deter- 
mine the phase change front position histories. The computed and film derived 
position histories were compared. Physical interpretations of film and tem- 
perature data are given and the validity of analytical models is discussed, along 
with possible reasons for model output discrepancies. 



CHAPTER II 


LITERATURE SURVEY 
Introduction 

The literature survey considered only those studies that related 
directly to the thermal design of a finned thermal capacitor. Special emphasis 
was placed on space vehicle applications. Information of interest pertaining 
to phase change materials is discussed, followed by a survey of applicable one- 
g and zero-g studies, and concluded with a brief description of previous experi- 
mental programs. 

Analytical comments are abbreviated, as volumes of material are 
available pertainir^ to solutions of the freeze/melt problem. There was no 
intent in this research to cover the large body of literature relating to macro- 
scopic and microscopic aspects of phase chaise. These are normally included 
in the materials disciplines (e.g., phase diagrams, nucleation, supercooling, 
superheating, etc.). For information concernii^ these areas, an excellent 
text by Chalmers [l] is available. 

Phase Change Material 

In general, PCMs may be catagorized into two groups: paraffinic and 
nonparaffinic. The paraffins are the most widely used PCM in current 



aerospace applications, while the nonparaffins are being used in nonspace and 
future aerospace study applications. The nonparaffin group incorporate a large 
body of materials, including: water, hydrated salts, organics, acids, and 
synthetic mixtures. 

Some researchers [2, 3, 4] have defined the qualities that a good PCM 
should exhibit, which are as follows: 

1. high heat of fusion 

2. proper melting point 

3. high. thermal diffusivity 

4. high coefficient of thermal conductivity 

5. non-corrosive 

6. low coefficient of expansion and small volume change during phase 
and lattice changes 

7. stable 

8. high flash point 

9. good wetting characteristics 

10. low cost 

11. minimum of void formations 

12. relatively pure 

13. non-toxic 

14. readily available 
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Grodzka [2l screened a large field of nonparaffin candidate PCMs covering 

a range of freezing/melting points from 16.12^C (61°F) to 83,96°C (183T). 
She designated a salt hydrate and three nonparaffin organics as prime PCM 
candidates. 

Grodzka and Fan [5] also examined the limits of thermal stability and 
long term thermal cycling effects on certain nonparaffins considered suitable 
for use as PCMs. They found that long term thermal cycling can result in a 
buildup of impurities that may interfere with efficient operation of the PCM, 
They concluded that this also could cause eventual destruction of nucleation 
catalyst and result in stratification of impurities. 

Designers have found normal paraffins to be a good PCM for space 
applications. Normal paraffins are hydrocarbons whose generalized chemical 
formula is given by ^ selected list of normal paraffins that are 

commercially available, along with their published freeze/melt points, are 
given in Table 1, 

Paraffins with an even number of carbon atoms between 20 and 32, and 
those with any odd number of carbon atoms exhibit a lattice transition. The 
even numbered carbon atom paraffins exhibit this transition near their melting 
point, whereas odd numbered paraffins exhibit the transition in the solid state, 
as much as 16^C (30''F) below their freeze/melt temperature. Broadhurst [6] 
has shown that the energy associated with this transition is subtractive from the 
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TABLE 1 . NORMAL PARAFFINS 


Name 

Chemical 

Freezing Point 

Formula 

°F 

°C 

Undecane 

Cif H24 

-14. 1 

-25.6 

Dodecanc 

C(2 ^26 

14.7 

-9.6 

Tridecane 

C13 H28 

22.3 

-5.4 

Tetradecane 

C,4 H30 

42.6 

5.9 

Pentadecane 

^15 H32 

49.9 

10.0 

Hexadecane 

^IG 

64.7 

18.2 

Heptadecane 

^17 H3g 

71.6 

22.0 

Octadecane 

^18 H38 

82.8 

28.2 

Nonadecane 

^19 H40 

89.4 

32. 1 

Eicosane 

^20 

98.2 

36.8 

Hcneicosane 

C21 H44 

104.9 

40.5 

Docosane 

^22 ^^ 4 f. 

111.9 

44.4 

Tricosane 

C23 II .48 

117.7 

47.6 

Tetracosane 

^24 ^^50 

123.6 

50. 9 

Pentacosane 

^25 ^^52 

128.7 

53.2 

Hexacosane 

C 2 G H 54 

133.5 

56.4 

Hcptacosane 

^27 ^56 

138.2 

59.0 

Octacosane 

^28 H58 

142.5 

61.4 

Nonacosane 

C29 Hgo 

146.7 

63.8 


C;jO H(J2 


Triacontane 


150.4 


(>5.8 














normal energy absorbed or liberated due to phase change. Paraffins exhibitir^ 
this phenomena have latent heats of fusion that are 10 to 20 percent below the 
latent heats exhibited by their carbon atom neighbors (Fig. 3) . Shlosinger 
and Bentilla [7] have noted that the freeze/melt point of normal paraffins 
increases with the number of carbon atoms (Fig. 4) . 
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Figure 3. Energy associated with solid-liquid phase change and solid- 
solid transition of even-numbered and odd-numbered paraffins. 

Designers have considered mixing paraffins to achieve a freezing point 
which is different from that of a pure material. In such cases, the question 
arises as to the resulting mixture’s freezing point, latent heat of fusion, and 
other properties. Bentilla, et al. , [8] constructed a two constituent phase 
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Figure 4. Temperature associated with solid-liquid phase change and 
solid-solid transition of even-numbered and odd-numbered paraffins. 

diagram (Fig. 5) for octacosane/eicosane mixtures. These data indicate an 
undisturbed mixture will melt/freeze over a range of temperatures. Further- 
more, additional data from this source, using the same constituents, shows 
that there is some selective freezii^ at the individual phase change temperatures 
of the components (Fig. 6, 7, 8). Nagel‘ and Shelpul^ independently produced 
similar results with dodecane/tridecane and eicosane/docosane mixtures, 

^Nagel, R. : Unpublished Data. McDonnell Doiaglas Company, St. 

Louis, Mo. 1972. 

^Shelpuk, B. : Thermal Design of the Lunar Communications Relay 
Unit. Unpublished RCA Presentation, 1971. 
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COMPOSmON <MOtE FRACTION OF OCTACOSANEl 


Figure 5. Temperature-composition phase diagram for the octacosane- 
eicosane binary system. 

respectively. However, Prenger [9] has shown that a discrete phase change 
temperature proportional to the mass ratio of the constituents can be achieved 
by constant mixing of the liquid phase for a tridecane/tetradecane mixture 
(Fig. 9) . Shelpuk has observed that on mixing eicosane and docosane (Fig. lO) , 
both of which have latent heat of fusions of approximately 244 J/gm (l05 BTU/ 
lb) , the heat of fusion of the mixture was reduced as low as 139 J/gm (60 
BTU/lb) . The uncertainty of a mixture melt/freeze point, added to the 
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Figure 6, Heating -cooling curves for pure eicosane. 




Figure 7. Heating -cooling curves for (80% -50%) octacosane - (20% 
50%) eicosane binary system. 
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Figure 8. Heating-cooling curves for (20% -8%) octacosane - (80% 
-92%) eicosane binary system. 
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Figure 10. Heat of fusion for eicosane/docosane mixture, 
reported degradation in the latent heat of fusion, couple to make paraffin 
mixtures unattractive from a heat transfer viewpoint. However, other consid- 
rations discussed in Chapter 3 may override thermal considerations. 

Property data for PCMs can be found in various texts, most ostensibly 
in references 10, 11, and 12. Hale, Hoover and O'Neil [13] have tabulated 
some selected property values, as well as cost and compatibility information 
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for a group of paraffin and nonparaffin materials. The Boeing Company^ has 
tabulated some thermal properties (excluding the thermal conductivity values) 
of paraffins for a range of temperatures. A more comprehensive property 
tabulation was compiled by the Phillips Petroleum Company [ 14] . However, 
certain property values do not agree with those quoted in earlier references 
[7, 11] . This disagreement in property value was not surprising in light of 
recent findings [15] which indicate that paraffin properties can vary significantly, 
depending on the production techniques, the levels of impurity, and the pre- 
treatment. It must be noted that a single paraffin may be commercially 
available in a number of purity levels, which are not specified in most property 
data. Commonly, technical grade paraffin (> 95 percent pure) is used because 
of its lower cost. However, more expensive laboratory grade paraffins (> 99 
percent pure) are available. These higher purity paraffins offer a more 

predictable freezing point, as well as higher latent heat of fusions than techni- 
cal grade. 

The most comprehensive literature search and property measurement 
survey to date is given in reference 7. Unfortunately, this study considered 
only paraffins with freeze/melt points between 5.6 °C (42° F) and 65.6°C 
(150°F) . A comprehensive tabulation of paraffin properties excerpted from 
numerous documents is given in Appendix A. 

3 

The Boeing Company, unpublished internal memorandum. 
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One-g Performance 


For estimating the one-g performance of a closed cell phase change 
system, a determination of the relative importance of the three basic modes of 
heat transfer occurring internal to the capacitor must be made. Radiation for 
most applications at normal temperatures is negligible, leaving only conduction 
and convection as possible transfer modes. 

Causes for convective motion in one-g could be due to: 

1. buoyancy forces 

2. surface tension forces 

3. volume change forces 

4. external forces 

5. other less common effects, such as electric fields, atomic radia- 
tion, chemical reactions, etc. 

Neglecting the affects of the last two, only buoyancy, surface tension, 
and volume chai'^};e remain. 

In early studies of buoyancy driven convection. Lord Rayleigh [I6l re- 
lates the non-dimensional Nusselt number (Nu) to the Grashoff (Gr), Prandtl 
(Pr) number product. The significance and definition of these non-dimensional 
quantities are given by Knuds en and Katz ( 17l as 

hL / . \ 

Nu = dimensionless temperature ratio - 


1 !) 



Gr = (Bouyancy forces) (inertia forces) _ gi/3 [ Tw - T®) 

(Viscous forces)^ ^ ^ 

Pj. = Molecular diffusivity of momentum _ 

Molecular diffusivity of heat (^) 

Later, the product of the Grashoff and Prandtl number became known as the 
Rayleigh Number (Ra) . 

Numerous experimenters have conducted studies that were related to 
hydrodynamic instability caused by buoyancy as it effects the heat transfer 
process, including the studies cited in references 18 through 29. Silveston 
and O'Toole [30] verified the Nusselt number versus Rayleigh Number corre- 
lation for a fluid confined between two parallel plates (Fig. 11) . A number of 
authors, including Catton and Edwards [31 1 later showed the effect of L/D ratios 
for closed cells (Fig. 12) , where L is the cell height and D represents the 
fin cell spacings. 

The critical Rayleigh Number is defined as the value at which convec- 
tion begins. In the region below the critical value, the heat transfer is sub- 
stantially by conduction only, Nu - 1. The critical value, which depends on 
the boundary conditions, is 1708 for a fluid confined between two infinite 
horizontal, isothermal, conducting walls; and 720 for a non-conducting walled 
container [29], As shown in a bounded cell (Fig. 12a), this critical value 
tends to increase as the cell sides approach one another (i.e. , as L/D in- 
creases) , and also tends to increase as the walls become more conducting 
(Fig. 12b). 
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RAYLEIGH NUMBER, Ra 


Figure 11. Silveston's experimental results in the neighborhood of 
instability in various liquids. 
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Figure 12. Region of influence of L/D on heat transfer. 
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Although scientists have been aware of surface tension driven convection 
for some time, quantitative studies of this phenomena are scarce. The flow 
patterns created in Benard's classic experiment [32], which produced cellular 
circulation patterns in a very shallow liquid was initially attributed to buoyancy 
effects. However, Block [33], Pearson [34], Scriven [35], and Sternling [36], 
later proved that this phenomena was due to surface tension driven convection. 
The initial discovery of this phenomena is attributed to Marangoni [37] , and 
the commonly used term " Marangoni Flow" was coined to name this phenomena. 

The nature of the Marangoni flow, as discussed by Young, Goldstein 
and Block [38] , is that temperature variations across a free gas/liquid inter- 
face changes the shear forces of the surface. This is due to the dependency of 
surface tension on temperature, which is estimated by Gambrill [39] to be 
linear: 

<7 = + bT . (4) 

The coefficient "b" is negative, so that an increase in temperature at the sur- 
face is accompanied by a subsequent decrease in the surface tension. Hershey 
[40] has shown qualitatively that a depression in the surface at a local hot spot 
results, causing the liquid to flow away from the hot zone and toward the cold 
zone (Fig. 13). Subsequently, McGrow and Larkin [41] have photographed 
this effect for a number of configurations, producing dramatic verification of 


the phenomena. 


HOT FLUID 



Figure 13. Surface tension convection patterns. 

The level of surface tension driven convection is correlated by using 

the Marangoni number, Ma, which is given by: 

_ Surface tension forces 
Viscous forces 



Ra Ma 

Ra Ma 

cr cr 


( 6 ) 
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This assumes that buoyancy and Marangoni driven convection are additive. 
However, Grodzka [43] noted that experimental data on buoyancy and Marangoni 
surface effects indicate that an additive theory does not follow. Regardless, 
the convective currents caused by Marangoni flow in one-g are usually small 
compared to those caused by buoyancy. Pearson [34] has shown that for most 
fluids at normal temperatures, a liquid thickness of one cm or less must be 
attained before Marangoni effects overshadow buoyancy effects. 

For systems with small characteristic dimensions, the non-dimensional 
Bond number. Bo, is given by: 

Bo = gravity forces ^ pgl} . 

Surface tension forces o ' ' 

Appendix B outlines a typical calculation usii^ nonadecane. This number is some- 
times used to qualitatively evaluate the relative importances of Marangoni effects 
as compared to buoyancy effects. From equation (7), it follows that a low Bond 
number indicates a high degree of surface tension effects. 

Although surface tension effects exist at all vinlike interfaces, only the 
liquid/gas interface (as opposed to liquid/liquid and liquid/solid interfaces) is 
expected to produce appreciable resulting flows. However, no proof of this 
conjecture was noted in the literature surveyed. 

Volume change driven convection can be caused by the phase change 
process. During freezing, the new layer of frozen material at the interface 
tends to contract, since the solid density is usually greater than the liquid 
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density. Consequently, the liquid near the front will flow toward the interface 
to fill the volume shrinkage caused by solidification. During melting, the liq- 
uid at the interface tends to flow away to allow for the volume created by melt- 
ing. Tien and Koump [44] have stated that this effect will cause both the freez- 
ing and melting processes to be retarded. This is due to the ingress of warm 
fluid during freezing and the egress of warm fluid during melting (Fig. 14) . 
However, if circulation is created in the melting processes, warm fluid could 
be drawn into the interface thereby augmenting the melting process. 
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Volume change driven flow patterns. 


Reference 44 has also shown in a computational exercise for a ficti- 


tious system where the solid density is 25 percent greater than the liquid den- 
sity, that only a 10 percent reduction in the freezing rate occurs due to volume 
change effects. Since paraffins experience only a 5 to 10 percent volume in- 
crease on melting, it can be inferred that an even smaller effect can be 
expected in paraffins for similar conditions. 
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In summary, the only modes of heat transfer that are involved in a 
typical one-g static phase change process are conduction and convection. 
Convection is caused only by buoyancy driven currents for reasonable con- 
tainer sizes. Surface tension effects are negligible, except for very thin films, 
and volume change effects may be ignored. 

Zero Gravity 

When comparing environs of near earth orbital space with that of the 
earth, several differences are noted in the heat transfer process. In earth 
orbit, the reduced gravitational force is nearly balanced by the centripetal 
orbital force, creating an effective zero-g environment. Effects caused by 
reduced pressure, radiation field, meteoroid bombardment, and three- 
dimensional spacecraft maneuvers are also possible. 

Restricting the hypothetical capacitor under consideration to be a 
hermetically sealed container that is isolated from exterior thermal effects by 
insulation and anti -penetration shields and to be aboard a non-maneuvering 
vehicle, then reduced gravity remains as the only important alien effect. 

The primary effect of reduced gravity on the heat transfer mechanism 
is in the lessenir^ or elimination of buoyancy convection. Typical measure- 
ments of the net gravitational acceleration force on a spacecraft indicates 
values of the order of 1 x lo ^ g. Using nonadecane paraffin properties, a 
typical Rayleigh number at zero-g is: 

Ra = 0.01 L3 at (8) 
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Details of this calculation are given in Appendix B. From this, it is apparent 
that for reasonable container sizes (i.e., less than 15-cm (6-in.) cell depths), 
the temperature difference across the liquidus portion of the cell must be 
= 444°C (SSO'F) to produce buoyancy driven convection. As this temperature 
range is well in excess of normal operatir^ values, buoyancy stimulated 
convection may be considered negligible. 

At one-g, Marangoni flow driven convection is normally unimportant, 
however, at zero-g this is not necessarily true. As discussed earlier, liquid 
flow caused by surface tension at a liquid/vapor interface may occur. 

A group of recent experiments, carried out during the mission of 
NASA' s Apollo-14, revealed that the surface tension driven phenomena in zero- 
g is a reality and can produce significant convection [45, 46 ]. Using data from 
these experiments, Grodzka [46] has plotted the relation existing between the 
Marangoni number and the ratio of effective thermal conductivity to actual 
conductivity (Fig. 15). These data show a rapid increase in the convective 
level at a Marangoni number slightly greater than 300. The temperature 
difference in the Krytox test liquid was only 2.5°C (4. 1°F) in this instance. 

Close examination of this data stimulates some questions. Applying 
these data to nonedecane paraffin contained in a cell with a 15.24-cm (6-in.) 
characteristic dimension at a Marangoni number of 300, a of 12 is pre- 

dicted in Figure 15. This is a very high convective level. The temperature 
difference across a parallel plate system required to reach this level is only 
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Figure 15. Heat transfer characteristics of Benard cells. 

O.OOll'C (0.002°F). However, for this same system a Rayleigh number of 
4.32 X 10® is predicted at one-g. From this, a of only 2.5 is indicated, 

which is a convective level well below that predicted by the tentative surface 
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tension data. This conflicts with earlier investigators, who reported that for 
this condition at one-g, Marangoni driven convection is unimportant. However, 
this could be explained, at least partially, by the fact that the Apollo-14 data 
is taken for a free surface, whereas the Figure 11 correlations were taken for 
top and bottom bounded cells. 

Using nonadecane paraffin properties at the phase charge temperature, 
the Marangoni number is given in Appendix B as 

Ma = 2374 ATL , (9) 

This indicates that for reasonably sized capacitor cells, a temperature 
difference of only 0. 003° K (O. 006° F ) is necessary for the onset of Marangoni 
driven convection in a paraffin filled capacitor at zero-g. This fact indicates 
that Marangoni convec tion ca i be appreciable at zero-g. 

Since the magnitude of volume change effects discussed earlier were 
low for freezing (where buoyancy driven convection is negligible) as well as 
melting, it may be implied that the volumetric effects are negligible for pure 
conduction as well as convective processes. Consequently, the volumetric 
effects on the heat transfer process are negligible for zero-g operation. Heat 
transfer modes available to a thermal capacitor at zero-g are then reduced to 
Marangoni driven convection and pure conduction. 

A secondary effect of reduced gravity which can significantly alter the 
heat transfer process is the ullage gas position. In a typical rectangular cell 
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containing ullage volume, the ullage gas may configure itself in any of a number 
of possible modes. A number of possible cell ullage locations are given by 
Reference 16, those of more importance are shown in Figure 16. When heating/ 
cooling from the bottom, only configurations d and c (Fi,.,. IG) would alter 
the normal process; however, when heating/cooling from the top or sides, con- 
figurations a, b, d, and e would all reduce the rate of heat transfer due to the 
insulating effect of the ullage gas. Small bubbles occurring in the liquid could 
induce convective currents. In a zero-g field, this motion could be caused by 
the Marangoni flow phenomena causing the bubbles to migrate toward warm 
zones. Also bubbles could be entrapped in the freezing solid (configuration d 
of Figure 16) decreasing the apparent thermal conductivity of the PCM. Fortu- 
nately, paraffins have the property of being good surface wetters, which in zero- 
g tends to force the ullage to form in the center of the cell (configuration c of 
Figure ir^). However, insufficient quantitative data are available on these 
phenomena to determine affect on the heat transfer process under given 
conditions. 

Grodzka I 2] has examined effects of the space environment on the 
microscopic processes. She concluded that complex coupling effects between 
phase change kinetics and various [)ossible modes of convective motion cannot 
be predicted accurately without actual flight data. She also concluded that the 
magnitude of magnetic and electric fields likely to be encountered in the earth’s 
orbit are not expected to alter i)hase change bi^havior significantly from that 
observed on tiie earth. Finally, she stated that radiation fields encountered 

:u) 




Figure 16. Zero-g ullage configurations, 
in earth orbit are expected to have little effect, except perhaps in the case of 
organic PCM where long time exposures will result in buildup of impurity. 
Although the later statement applies to paraffin, no definitive information on 
this effect could be found in the literature surveyed. 

Although insufficient reduced gravity data are available to corroborate 
these findings, some observations have been made as to the mechanisms of 
heat transfer in a thermal capacitor. When a free surface is present and 
thermal conditions are proper, Marangoni convection may be present. In- 
sufficient data are currently available to determine quantatively the influence 
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of this phenomena. In the absence of surface tension effects^ conduction is the 
sole process governing heat transfer. However, the conduction process may 
be altered by the ullage gas location. Furthermore, secondary effects of 
property degradation due to radiation may occur if the PCM is organic. 

Experimental Investigations 

The publications on experimental studies of thermal capacitors are 
limited. Altman, et al. , [47], Fixler [3, 48], Bannister [4], and Kaye, et 
al. , [50] experimented with various capacitor designs, proving the feasibility 
of the capacitor thermal control approach. The most comprehensive experi- 
mental data published were by E. Bentilla, K. Sterrett and L. Karre 1 8] . 

Usii^ four different paraffins, they tested a variety of packaging techniques 
including aluminum wool, aluminum foam, copper foam and aluminum honey- 
comb. They concluded that the temperature rise of the paraffin container 
during melting was due to the insulating effects of the liquid paraffin. They 
also noted that the experimental data indicated a higher thermal conductivity 
than that specified for the property data. Freeze data indicated that the conduc- 
tion model which they used to attempt to match test data, repeatedly yielded 
lower freeze rates. Data from melting tests, showing maximum coldplate 
temperature rises above the melting point temperature versus time for varying 
heating rates in an aluminum honeycomb system, are shown in Figure 17. 
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MIN. TEMPERATURE RISE, BOTTOM OF MAX. COLD PLATE TEMPERATURE 

PACKAGE, ABOVE^MELT PT. (300‘’K) r,Se AB0>/E MELT PT. (300°K) 


MODEL 13 . MOD 1 

ALUM. HONEYCOMB FILLER MAT'S 

10.8% TOTAL VOID VOL. (0.000061m) 

SHEET THICKNESS .0024 in. (0.00104m) 

NOMINAL MAX. SHEET SPACING 0.041 in. (.104cm) 
PACKAGE THICKNESS. .5 in. (1.27 cm) 

WT. FUS. MAT. (0.075 Kg) 

WT. PACKAGE/WT. FUS. MAT. 78 2% ^ , 

NOMINAL COLD PLATE AREA 0.0935 ft^ (0.0086m'^) 
TOTAL WT / AREA 3.16 LB / FT2 (15.4 KG/m^) 
ADHESIVE AM. CYANAMID COMPANY FM-1000 
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Figure 17. Adiabatic experiment test results. 
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Various investigators have produced microscopic photographs of the 
freezing interface in paraffins [51, 52, 53] • Magnification photography (175 
power) by Bannister and Richard [52] shows the freeze front is composed of 
a multitude of dendrites of varying geometry (Figure 18), All paraffin 
freeze data surveyed contained references to dendrite formations,. Dendrite 
arms appear much like cilia or fine fur to the naked eye, and they tend to grow 
in the direction of heat transfer. Westwater and Thomas [53] , have shown 
that the physical shape of these dendrites are at least partially dependent on 
the rate of freezing. Other dependencies were not well defined, with little 
information available which would allow correlation of dendrite geometry to 
cooling rates. 
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Figure 18 -A. Freeze front dendrite formations. 
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Figure 18-B. Freeze front dendrite formations. 



CHAPTER 3 


APPLICATIONS 
Summary of Applications 

The thermal capacitor design concept has a multitude of possible appli- 
cations for both commercial and aerospace systems. A number of actual appli- 
cations (along with problems encountered during design), as well as possible 
applications currently under study are discussed. 

A list of these applications are given below: 

A. Commercial 

1. Actual 

a. children's kool aid freezer 

b. portable freezer chest 

c. thaw warning indicator 

2. Possible 

a. thermal insulation for homes 

b. air conditioner supplement 

c. large solar thermal energy collector 

d. laser heat pulse absorber 

e. hospital food cart heat storage supplement 
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B. Military/Aerospaee 


1. Actual 

a. Titan transtage suction line temperature control 

b. Pegasus III temperature reference 

c. Lunar roving vehicle electronics temperature control 
cl. Skylab coolant loop augmentation 

e. Skylab waste refrigeration 

f. Poseidon gimbal shipping container 

2. Possible 

a. space shuttle wing leading edge temperature control 

b. space station space radiator augmentation 

c. astronaut suit cooling 

d. extra planetary probe thermal control 

e. Air Force satellite 

f. F-4 aircraft 

g. short term thermal control of army shelters 

Commercial Applications 

The commercial applications include a plastic container surrounded 
with a hollow annulus which contains a PCM. Placing the container in the freezer 
compartment of the refrigerator cools the PCM so that a water suspended bev- 
erage, such as kool aid, may be frozen in a slurry much like the corner store 
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Koolee' or ’Icee", A freezer chest is available which has a PCM liner that 
may be removed and refrigerator cooled, to be later used to store drinks and 
perishable items for lot^ trips. 

The mushroom industry commonly freeze their product for shipment. 
Thawing causes permanent quality degradation. As a result, a thaw warning- 
indicator has been developed using a phase change device. A salt PCM is 
attached in a bag to the side of mushroom containers. The frozen salt con- 
tains a chemical additive which causes litmus paper, visible to container handlers, 
to change color if the PCM melts [53]. 

Telkes [55, 56, 57] has, for many years, conducted research on appli- 
cations of PCMs to residential dwellings to minimize thermal oscillations on 
their heating/cooling systems. Brine and other salt solutions have been typi- 
cally used in this research. Also under study are applications of PCMs to 
large solar collection systems, with surface dimensions upward to 10 square 
miles. Rapid energy conversion inherent in laser systems are prime targets 
for a thermal capacitor type application. Hospital carts required to store 
warm food for extended periods may be fitted with PCMs. 

Aerospace Applications 

The most widespread use of thermal capacitor techniques has been in the 
aerospace field. One of the earliest applications in this field was used on the 
Titan water jacket, applied around the engine attitude control system oxidizer 



lines to limit temperature excusions below 32° F, An Kicosane PCM was used 
on the Pegasus III flight [58] to act as a constant temperature reference 
junction for a surface coating experiment* 

Recently, moon flights of Apollo 15, 16, and 17 have transported 
thermal capacitor devices, which were mounted on various heat producing 
electronic components (Fig. 19) in the Lunar Roving Vehicle (LRV) . The 
wheel drive control electronics (DCF) and the navigational signal processing 
unit [SPUl, designed by Flliott and Paoletti [59], used paraffin filled capaci- 
tors for cooling. These two devices had sufficient masses to allow energy 
storage from their respective components during LRV sorties, so that compo- 
nent temperatui’e excursions were kept below allowable lex els. Fxtn more 
unusual than the thermal capacitor applications was the thermal strap techni- 
que of transferring heat from the capacitor to a set of second-surface mirror 
radiators. After each sortie, the LRV was parked so that these mirror 
radiators were in the shade. A thermostatically controlled shutter was then 
opened, automatically closing after the two capacitors were refrozen. 

'I’he most significant problem reported in dexelopment of these capacitois 
was due to paraffin property variations. The dispersion in property data was 
so great that after being received, the paraffins xx^ere subjected to repeated 
degassings under reduced pressure to stabilize piopeities. 

Shelpuk, designer of the LRV lunar communications relay unit (LCIU ) 
capacitor, also noted paraffin property dispersion. He found that synthetically 
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Figure 19. Space vehicle capacitor configurations. 




produced paraffins had better property stability than naturally produced paraf- 
fins. He used a mixture of two paraffins in the LCRU. This allowed softening 
of the solid paraffin at temperatures below melting, thereby minimizing stress 
buildups in the capacitor cells. 

One interesting aspect common to all LRV capacitors was that all were 
completely filled at temperatures above their maximum operating values, 
leavir^ no ullage. This technique assumes gas does not leak into the paraffin 
container after cooldown. This technique inherently relies on the quality qf 
workmanship, since significant gas in-leakage could cause a large pressure 
buildup during subsequent heating. A strong argument for this technique is the 
fact that flight data telemetered back from the moon indicated that all LRV 
capacitors functioned properly. 

There were fifteen thermal capacitors used on the first manned earth 
orbitir^ space station, Skylab. A total of five of these were located in two 
cooling systems, thermally controlling liquid temperatures. Two of these five 
capacitors were located on the primary Skylab space radiator fluid outlet, 
while the other three were located at the outlet of a smaller radiator used to 
refrigerate food and biological waste samples. Ten other capacitors were 
mounted in the bottom of trays used to store human waste samples while they 
are being transported from the Skylab Vehicle to earth. 

The most significant problem encountered during development of these 
capacitors was the large structural stresses which occurred in the paraffin 
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housing during rapid melting. In the initial design, even though ullage volume 
was provided, stresses were sufficient to cause local buckling of the housing. 
This overstress conditioning was caused by the inability of the paraffin mate- 
rial to expand into the ullage volume. This occurred only when the capacitor 
was oriented in such a manner that the ullage volume migrated to a position 
away from the melting cells. As a result, the non-sealed corrugated fin-cell 
arrangement initially used was replaced with a honeycomb fin-cell arrangement. 
The honeycomb cells could be hermetically sealed on top and bottom, restricting 
ullage migration within the small honeycomb cell regardless of orientation. 

This particular design by the McDonnell-Douglas Company incorporated 
a 20 percent ullage allowance in each cell. Although this eliminated the depend- 
ency on workmanship required in LRV capacitors, the large gas volume could 
interfere with the heat transfer. Additionally, this design required a heavier 
walled housing than was required by an evacuated concept. 

A tabulation of the NASA Space Vehicle Applications is given in Table 
2, itemizing the subsystem controlled and the capacitor functions. Pictorials 
of NASA space vehicle capacitor systems with their associated attachment 
arrangements are shown in Figure 19. Table 3 itemizes interesting design 
features of these capacitors. 

The capacitor concept was used also for temperature control of the 
Poseidon engine gimbals during shipment. This device used a mixture of 
organic acids, consisting primarily of Palmitic and Stearic acids, to control 
the gimbals at temperatures below 55°C (131°F). 
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TABLE 2. NASA SPACE VEHICLE CAPACITOR APPLICATIONS 


Vehicle 

Subsystem Controlled 

Function 

Pegasus in 

Surface coating contami- 
nation experiment 

To provide a constant 
temperature reference 
for a thermocouple 

Apollo 15 
Apollo 16 
Apollo 17 

l) Lunar roving vehicle 
drive control 
electronics 

To provide regenerable 
heat absorbing source 


2) Lunar roving vehicle 
signal processing unit 
electronics 

To incre‘ase thermal 
capacitance available 
to signal process unit 


3 ) Lunar roving vehicle 
Lunar communica- 
tions relay unit 

To provide regenerable 
heat absorbing source 

Sky lab 

1 ) Airlock coolant 
system 

To augment space 
radiator performance 
and limit inlet temper- 
ature to downstream 
heat exchangers 


2 ) Refrigeration coolant 
system 

Enhance space radia- 
tor performance to 
limit cold inlet temper- 
ature to downstream 
components 


3 ) Urine freezer 

To provide temporary 
cooling for urine 
sample when being 
transported between 
Skylab and earth 


44 
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Another possible space application includes the thermal protection 
system for the. space shuttle orbiter wir^ leadir^ edge [60] . This concept 
utilizes a PCM which melts near 176. 8° C (350°F) to limit the extreme temper- 
ature surges normally encountered on frontal areas of supersonic vehicles. 

Phosphonium chloride, a high latent heat PCM, has been studied for 
use in Astronaut Extra Vehicle Activity (EVA) cooling equipment. Schelden 
and Golden [61] , have studied the capacitor technique for application to 
thermally controlled devices in Jovian and Venusian extraplanetary space 
probes. The Air Force has examined this technique for classified applications 
to unmanned earth orbiting satellites, and as a special cooling system for modi- 
fied F-4 aircraft. This technique of thermal control is also beir^ investigated 
by the Russians (Vaselou, Kalisheva, and Telepin) for use on gravity measuring 
devices [62]. 

A salt hydrate type PCM has been studied for application to Army 
shelter refrigeration systems. This PCM material is designed to maintain a 
relatively constant temperature in the absence or failure of the active refrig- 
eration system. 
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CHAPTER 4 


EXPERIMENTAL APPARATUS 
Test Item 

The initial capacitor design goal was to produce a test item which 
closely simulated the heat transfer paths expected in a typical flight thermal 
capacitor while still allowing observations of the melt/freeze phenomena. 
Unfortunately these two goals were found to be incompatible, since the high 
resistance to heat transfer offered by most transparent materials disallows 
simulation of transfer paths normally available to an all metal capacitor. As 
a consequence, the final design stressed capacitor visibility, relegating geo- 
metrical simulation as a secondary goal. 

The resulting design incorporated a structure fabricated almost entirely 
of Plexiglas* with aluminum surfaces transferring heat to and from the paraffin 
into the coolant. The capacitor was fabricated in two parts: the fluid passage 
and the paraffin housing (Figures 20 and 21). 

The overall fluid passage envelope dimensions were 26-cm ( 10. 25-in. ) 
long by 14.61-cm (5.75-in.) wide by 2.54-cm (l-in.) depth. Three 0.635-cm 
* 

A Rohm and Haas tradename 
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Figure 20. Thermal capacitor test item, top view. 





(1/4 -in. ) diameter tapped holes were inserted parallel to the housing length 
(i.e. , in direction of fluid flow) and into both ends of the initially solid plexi- 
glas block. These holes extended 3.81-cm ( 1. 5-in. ) in from both ends of the 
housii^. These three circular passages entered a 1.27 -cm (l/2-in. ) long by 
1.91-cm (3/4-in.) deep mixing plenum, extendii^ the width of the housing. 

These inlet and outlet plenums were connected by 10 fluid channels extendir^ 
the remaining 15.875-cm (6.25-in.) of the housing length. The channels formed 
individual flow cross sections of 1.27 -cm (l/2-in. ) by 1.91-cm (3/4-in.), 
separated by 0.318-cm ( l/8-in. ) thick plexiglas partitions. After milling these 
flow passages out, a 0.635-cm (l/4-in. ) thick bottom and 0.953-cm (3/8-in.) 
thick sides remained, with the top over the 10 fluid flow channels left open. 

The upper paraffin housing was fabricated by attaching 0.02-cm (0.008- 
in.) thick 5052-T30 aluminum fins to a 6061-T6 aluminum bottom plate (Figure 
22). The fins were 13.34-cm (5.25-in.) long and6.35-cm (2.5-in.) deep, 
notched at both ends to accommodate 0.48-cm (o. 188-in. ) deep grooves cut 
in the plexiglas sides. The fins were vertically mounted on the bottom plate 
by first soldering a 0.24-cm ( 3/32-in.) long base along the length of each fin. 
Three 0. 08-cm ( l/32-in. ) diameter rivets were then driven through the fin 
foot into the bottom plate, while being sealed at the same time by an epoxy 
filler with a silver matrix coating between the foot and the bottom plate. 

Four 1.9-cm (3/4-in.), four 1.27 -cm (l/2-in.) and three 0.635-cm 
(1/4 -in. ) fin spacings were used in the paraffin housing. The spacings were 
* 

Emerson and Cummings Co. , Ecco Bond 56C 
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FIN SPACING 



Figure 22. Typical fin mounting arrangement, 
symmetrically arranged with two 1.91-cm (3/4-in.) spacing on each of the outer 
ends, two 1.27 -cm (l/2-in. ) spacing following these, and three 0.635-cm 
(l/4-in. ) spacing at the housing center section. 

Then 0.95-cm (3/8-in.) grooved sides and 0.635-cm (l/4-in. ) thick 
ends were bonded together and attached to the bottom plate, using epoxy.* 
Finally, the bottom plate of the paraffin container section was epoxied to the 
open top of the fluid housing (i.e. , on edges and top of partitions) so that the 
fin lengths were perpendicular to the direction of fluid flow. 

The top of the housing for paraffin was fabricated of 0.635-cm (l/4-in. ) 
thick Plexiglas, with two 0.95-cm (3/8-in.) diameter paraffin fill bosses, and 

*NARMCO 7343 was used for all plexiglas bonding 
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was mounted on top of the paraffin housing with epoxy. Screws 0.32 cm 
(1/8 in.) in diameter were used to secure the paraffin housing top, being 
installed into tapped holes which extended through the top into the vertical sides. 
At each corner underneath the fluid housing 2.54-cm (l-in. ) cubical plexiglas 
standoffs were bonded to minimize heat transfer from the supporting surface. 

After fabrication and instrumentation, two calibrated strips with 0.318- 
cm (l/8-in. ) indexing were attached to the face of the paraffin housing (as 
viewed by the camera) . A vertical strip was attached to the fluid outlet end, 
extending the full 6.50-cm (2.56-in.) height of the PCM housing. Along with 
the length of the housing, in the fluid flow direction and at the top, a 15. 24 -cm 
(6-in.) calibrated strip was also attached. These strips provided visual 
dimensional reference for filmed data. 

Instrumentation 

At various times during the experimental program, the system was 
instrumented with as many as 40 temperature probes, 6 differential temperature 
probes, an RPM counter (pump flow rate controller), a flow meter, and a 
pressure transducer. All thermocouple temperature transducers were fabri- 
cated from 0.025-cm ( 0.01-in.) diameter (No. 30 gage) chromel constantan 
wire with 0.08-cm (O. 031-in. ) diameter thermocouple beads. Both resistance 
temperature bulb* (RTB) and chromel-constantan type probes were used for 

* 

Rosemont platinum resistance bulb 
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differential temperature measurements. The flowmeter was a Potter 0.954-cm 
( 0.375-in.) turbine type which provided flow measurements over the range of 
zero to 1.26 x 10“^ mVsec (2 gal/min) . The pressure transducer was a 
Statham strain gage type. 

The test item was instrumented with 27 thermocouples, as shown in 
Figure 23. Three paraffin cells were instrumented, two 1.9-cm (3/4-in.) 
width cells and one 0.625-cm (l/4-in. ) width cell. The 1-9-cm (3/4-in.) 
width cells are located on each side of the capacitor, adjacent to the outermost 
cells. Each of these had one thermocouple attached to the paraffin/fluid 
separator plate, three thermocouples attached to the surface of the interior fin, 
and four thermocouples suspended into the paraffin. The 0.625-cm (l/4-in. ) 
cell was located in the mid-section of the capacitor and had three thermocouples 
attached to the upstream-side fin surface and three thermocouples were 
suspended into the paraffin. 

After the first series of tests, three externally mounted thermocouples 
were additionally attached. These thermocouples were centrally mounted on 
the outer paraffin housing width and height, located at the capacitor inlet/outlet 
and the capacitor center top. These measurements were used to determine the 
sensible heat absorbed/rejected by the plastic structure. 

In addition to the 27 thermocouples mounted on the capacitor, four 
thermocouples and two RTB's were installed in the fluid inlets and outlets. 
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Figure 23. Diagram of thermal capacitor instrumentation. 






Differential temperature readings were taken on these thermocouples between 
the inlet and outlet fluid. The flowmeter recorded the fluid volumetric flow- 
rate while a thermocouple mounted in the fluid at position adjacent to it re- 
corded the fluid temperature. All other measurements were used to monitor 
facility ii^jut parameters for control purposes and were not used for data 
reduction. 

All thermocouples which were suspended in the paraffin were covered 
with 0.24-cm (3/32-in.) diameter ceramic sheaths to rigidize them in their 
locations. The lower end of the sheaths were stripped away to allow a 0.32-cm 
( 1/8-in. ) ler^th of the thermocouple to extend out of the sheath. Figure 24 
shows a typical paraffin suspended thermocouple arrai^ement. 



FLUID/ PARAFFIN 
TRANSFER PLATE 


Figure 24. Typical thermocouple instrumentation. 
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All surface and fin mounted thermocouples were attached with epoxy. 


The fin surfaces were thoroughly cleaned and l^htly roi^hened with sand 
paper before application. The thermocouple beads were then flattened and 
pressed to the surface, using 1-mil. thick 0.954-cm (3/8-in.) by 0.954-cm 
aluminum tape patches. Epoxy was then applied over the entire patch/thermo- 
couple installation. Figure 25 shows a typical fin mounting scheme. 


FIN MOUNTING SCHEME 



Figure 25. Typical fin thermocouple mounting arrat^ement. 

The circuitry for all thermocouple installations inside of the paraffin 
housing were routed through the capacitor top using epoxy to seal the penetra- 
tions. The circuitry for thermocouple measurements of fluid temperature 
were routed through the bottom of the fluid passage and sealed with epoxy. 
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A chromel-constantan thermocouple accuracy study, performed by the 
Marshall Space Flight Center Test Instrumentation Group, indicates an abso- 
lute accuracy of 0.4"* C (O. 8'’F) or better. RTB^s were more accurate, but 
exhibited poorer response. RTB measurements were not used for evaluating 
data because of the transient nature of these tests. Thin wire chromel- 
constantan thermocouples were selected because it minimized heat leak into 
the paraffin housing, thereby maximizing temperature response, to give the 
largest possible voltage output per unit temperature rise (Fig. 26). 

An extensive study was conducted to minimize error in reading the 
differential fluid temperature across the capacitor. However, due to the low 
temperature difference 0.9°C (l.5^F), the accuracy of this measurement 
never reached an acceptable level to allow correlation of the differential 
temperature data. 

The raw data was routed to the digital computer for interpretation on 
the basis of calibration curves that had been programmed into the computer 
system. Tabulated printouts of data was available immediately during and 
after each test. After the first 15 tests, the computer program was reconstructed 
to allow computerized plotting of the temperature versus time test data. 

Test System 

During all testing, the test item and support hardware were located in 
a thermostatically controlled room. As a consequence, the external environ- 
ment for all tests was at temperatures between 21. 2 and 23. 9°C (70 and 75°F) . 
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Figure 26. Temperature -millivolt graph for thermocouples. 



Nonadecane, normal paraffin was used as the PCM during all tests. The 
nonadecane was procurred from the Phillips Petroleum Company in the techni- 
cal grade. 

The plexiglas test item was mounted in a closed fluid flow loop (Fig. 

27) . This allowed circulation of cold and hot fluid through the fluid passage at 
varying fluid flowrates. A motion picture camera was mounted in a position 
whereby paraffin phase change position with respect to time could be viewed 
and recorded (Fig. 28). 

The coolant flow system consisted of hot and cold heat exchangers, a 
small-fluid pump, and appropriate valves and tubing (Fig. 27). The fluid was 
Monsanto coolanol-15, a silicate ester, used because of its frequent applica- 
tion to space vehicle coolant loops. The pump was a Vikii^ 2.5 m^/sec x 10”^ 
to 3.8 m^/sec x io“® (4 to 60 gal/hr) variable speed constant volumetric 
flow device. The pump speed was electrically controlled so that a constant non- 
varying mass flowrate could be set as desired. 

The coolant temperature was controlled by routing the thermal capacitor 
test item outlet fluid to a parallel bank of a hot and cold ’’tube -in-shell" heat 
exchangers. The hot heat exchanger contained water which was heated with a 
3000 joule/sec (3 kW) chromalox calrod unit mounted internally to the heat 
exchanger shell. This imit was thermostatically controlled to maintain the 
heat exchanger shell fluid at a maximum of 82. 3“ C (180*F). The other leg of 
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Figure 27. Thermal capacitor test facility. 



Figure 28. Motion picture camera arrangement, 
the parallel bank was routed through the cold heat exchanger, containing a 
coolanol-15/dry ice slurry. The shell fluid temperature stayed at a relatively 
constant -79° C (-110°F). 

Hand valves were mounted on the heat exchanger outlets to control the 
ratio of hot to cold fluid mixed. Downstream of these valves, the tubing was 
joined. From this point, a common line was routed back to the thermal 
capacitor inlet. 
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All line routii^s were fabricated from 1.27 -cm (l/2-in. ) stainless 
steel tubing of 0. 127 -cm ( 0. 049-in. ) wall thickness. The three 0.635-cm 
(l/4_in. ) ports on the thermal capacitor inlet and outlet were manifolded 
together before the inlet and outlet tubing was connected. Prior to testing, 
flow checks were made to assure an even flow split between the three fluid 
connections. 

The output of a thermocouple mounted at the test item inlet was trans- 
mitted to a visual output meter. Using this output, hand valves were adjusted 
to maintain the desired fluid inlet temperature. 

The capacitor was insulated on both top and bottom with 5.08-cm (2-in. ) 
thick fiberglas bating. The interconnecting tubing was insulated with 2.54-cm 
(l-in. ) fiberglas wrap. 

A motion picture camera was mounted so that it viewed the capacitor 
paraffin housing from the sides where the majority of the instrumentation was 
located. The camera was situated so that its field of view included the entire 
capacitor plus a digital timer (Fig. 28) . During initial testing, the timer was 
located on the fluid outlet side of the capacitor. In later tests, it was relocated 
to the top of the paraffin housing. 

The camera used during initial testii^ was a Mitchell, which filmed 
real time at variable speeds between 6 and 24 frames per second. Due to the 
voluminous amount of data which was accumulated, with a real-time output, 
the camera was intermittently switched on and off. Later, a Cine Special 
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camera with a variable framing rate of 0.25 to 1 frame per second was used 
for intermittent filming. This allowed compact filming without compromising 
data acquisition. During all runs, color film was used to optimize film 
definition. 

The timer was a Cramer Controls Corporation device. The digital 
readout contained five digits: four whole digits and one decimal second digit. 

During initial testing, the capacitor was illuminated from the camera 
side of the test system. However, the lightii^ was later moved to the opposite 
side of the capacitor for better film definition. 


Testing 

Prior to starting the tests, the capacitor was filled with nonadecane by 
heating the paraffin and the capacitor separately to a temperature of 71.7° C 
( 160° F) . The housing was filled by simply pouring the hot paraffin into one 
of the two fill ports, while allowing the other to vent. 


The initial weight of the nonadecane required to fill the paraffin housing 
was 1.95 lbs. The paraffin weight was determined by using a specially con- 
structed load cell weighing arrangement. The paraffin weight was determined 
by weighing the full paraffin container plus the paraffin prior to fill; then 
weighing the container plus the remaining paraffin after fill. 


Wtj = 


Wt 


c+pf 


- Wt 


( 10 ) 
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The inlet and outlet fluid temperatures, alor^ with the test item struc- 
tural temperature measurements, were used to determine the heat lost or 
gained from the fluid during paraffin phase change. This was compared to the 
heat received by the paraffin during phase change plus the sensible heat change 
of the test item, assuming heat loss to be negligible. 

Nonadecane was selected because its melting point is slightly above 
normal room temperature, its physical properties were known, and the liter- 
ature search had showed that none of the previous experimenters had used 
this particular PCM . 

The band of inlet fluid temperatures was arbitrarily established over a 
range of ± 21. 1°C (70°F) which centered on the published nominal freezing 
point temperature. This required an inlet temperature range from -6.7 to 
71.1°C (20 to 160°F). The mass flowrate range from 12.6 to 37.8 g/sec 
( 100 to 300 Ib/hr) was used because it represented space vehicle typical 
coolant flowrates. 

Based on information given in Reference 70, the need for a guard vacuum 
around the capacitor to assure an adiabatic system was not considered neces- 
sary. This is primarily because of the good insulating qualities of the plexiglas 
walls [K - 0. 16 J/m sec “K (0.09 BTU/hr ft“R) 1 . However, all 
surfaces except those requiring exposure were well insulated during all runs. 
One test was conducted to verify the adiabatic assumption. During this run the 
entire capacitor was enveloped in insulating material. These test data were 
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compared to data obtained from a previous uninsulated test which was performed 
with identical input conditions. This comparison confirmed the adiabatic test 
item assumption. 

A typical test procedure is itemized in Table 4. The designation used 
for all thermal capacitor tests by the testing group was with the prefix 230. 

The total number of tests performed in conjunction with this study was 60 tests, 
230-1 through 230-60. However, of these runs, a large number of facility 
checkout tests were performed, during which unacceptable data was obtained. 
Asa result, only 20 of these tests are reported herein. Of these, 9 were melt 
tests and 11 were freeze tests. A summary of the test designation, along with 
pertinent input parameters for those tests deemed acceptable, is given in 
Table 5. Appendix D presents all pertinent raw test data. 

The primary problem associated with operation of the test apparatus 
lay in the control of the fluid inlet temperature. Using hand operated valvit^, 
the inlet temperature could be maintained no closer than 2.8°C (5°F), with 
short excursions reachii^ a variation of 5.6°C (10°F) from the nominal 
target value. These excursions stemmed from the inherent time lag that 
existed from the time that hot and cold fluid was introduced into the system, 
until the time when the temperature sensor began to pick up the perturbation. 

For example, when the operator received a low temperature reading on the 
inlet sensor, he would increase the flow with the hot side valve. There would 
be a time interval before the hot pulse could be detected by the sensor. 
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TABLE 4. TYPICAL TEST PROCEDURE 


Test Time 
Reference 


Procedure 


Freeze Test Only 


Common 


Pretest 


heat paraffin to 
a uniform tem- 
perature above 
melt point, i.e,, 
35 to 48.9^C 
(95°F to 120"F) 


Pretest 


Using the pump 
speed control 
establish the de- 
sired constant 
mass flow rate 


Pretest 


by causing the 
fluid to bypass 
the capacitor 
establish the flu- 
id inlet tempera- 
ture 


T 0 


Simultaneously, 
open the capac- 
inlet valve and 


• turn on the mo- 
picture camera 


and 


• initiate instru- 
mentation data 
sampling 


During entire test 


maintain a con- 
stant inlet fluid 
temperature to 
capacitor 


Test end 


deactivate mo- 
tion picture cam 
era data sam- 
pling, and turn 
off pump 
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TABLE 5. SUMMARY OF TEST INPUTS 


Test Designation 


Nominal Fluid Inlet 
Temperature 

Nominal Fluid 
Flowrate 

Type Test 

» F 

“C 


kg/ sec 

230-5 

Melt 

160 

; 71.2 

300 

0.038 

230-6 

Melt 

130 

54.7 

300 

0.038 

230-7 

Freeze 

20 

-6.7 

300 

0.038 

230-8 

Melt 

110 

, 43.4 

300 

0.038 

230-9 

Melt 

150 

65.6 

300 

0.038 

230-10 

Freeze 

50 

j 10.0 

300 

0.038 

230-11 

Freeze 

70 

21.1 

300 

0.038 

230-15 

Freeze 

30 

-1.1 

300 

0.038 

230-49 

^ Melt 

150 

65.6 

100 

0.013 

230-50 

Melt 

130 

54.7 

100 

0.013 

230-51 

Freeze 

20 

-6.7 

100 

0.013 

230-52 

Freeze 

1 

30 

-1.1 

30 

0.0038 

230-53 

Freeze 

50 

10.0 

200 

0.025 

230-54 

Freeze 

70 

21.1 

100 

0.013 

230-55 

Freeze 

70 

21. 1 

100 

0.013 

230-56 

Freeze 

30 

-1.1 

100 

0.013 

230-57 

Melt 

160 

71.2 

200 

0.025 

230-58 

Freeze 

30 

-1.1 

200 

0.025 

230-59 

Melt 

110 

43.4 

200 

0.025 

230-60 

Melt 

130 

54.7 

300 

0.038 
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probably due to slow fluid velocity, poor mixii^, and the large sensible mass 
of the system. Asa result, if the operator reacted too quickly or overacted to 
a temperature pertubation, it would induce a cyclic temperature, with a large 
excursion band. This was magnified by the control band placed on the thermo- 
stat control in the hot heat exchanger, which was quite large 
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CHAPTER 5 


ANALYTICAL TECHNIQUES 
Introduction 

The energy equations which must be solved in order to determine such 
values as interface position and temperature profile variation with time are 
formulated in this section. The techniques used to solve these equations and 
the equation input values are discussed. 


General 


The first law of thermodynamics for a system may be written in 
differential form as 


de 

dF 


dq 

dt 


dw 

dt 


( 11 ) 


For an open, unsteady system (neglecting radiation, nuclear, and electromag- 
netic contributions). Bird, Stewart, and Lightfoot [63] have expressed the 
first law as the following energy balance: 


of \ 
internal and \ 
kinetic energy - 
out by / 


/. 


Rate of 


( Rate of \ / Rate 

accumulation! f intern 
of internal + 

and kinetic / y uuL II ill uy 

^ energy j \ convection / \ convection 


internal and 
kinetic energj 
in by 


/ i. X \ / oet rate of \ 

"of heat ] ' 


addition by] 
conduction/ 


the system 
on I 
y \^urroundingsy 
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for a fluid, this reduces to a vector tensor form: 

(a) (b) 

^ ip((i + V, + 1 V. pv (a + '4 vni 

(c) (d) (e). (f) 

= _ (V- Q) + p(V-'g) - (V- PV) - [V- (r* V)] . (12) 

where: 

(a) is the rate of gain of energy per unit volume 

(b) is the rate of energy input, per unit volume, by convection 

(c) is the rate of energy input, per unit volume, by conduction 

(d) is the rate of work done on the fluid, per unit volume, by gravity 

forces 

(e) is the rate of work done on the fluid, per unit volume, by pressure 

forces 

(f) is the rate of work done on the fluid, per unit volume, by viscous 

forces 

For a typical cell undergoing phase charge, the kinetic and potential 
energy contributions are small, eliminatir^ the last term in a and b and the 
whole of d. Assuming the pressure forces are not allowed to build up appreci- 
ably, e may be neglected. The limited fluid velocities encountered for normal 
melting processes eliminates f. Thus equation (12) reduces to: 

l-(pfi) = - [V- (pVO)] - (V.Q) , (13) 

ot 

For a fluid at constant pressure, the internal energy differential is given by: 
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= -pdv + cpdT 

Neglecting viscous dissipation and assuming constant properties and an 
incompressible fluid, equation ( 13) may be written: 


Since Q = _ KV^ T, then from equation (15): 

0 ^ ^ ^ 

V^T-tV-VT) 

9t pc \ V j. , 

P 


( 15 ) 


( 16 ) 


for the fluid phase. For the solid phase, 

and (16) reduces to the familiar Fourier 

9T „ 

— = aV'T 


the convective term does not exist 
Conduction equation 

(17) 


Using equation (l7) , a formulation of a one-dimensional phase change 

problem may be constructed (Fig. 29). Consider the liquid phase to be at a 

constant temperature, at time t = 0, and the surface temperature at x - 0 

to be step changed from T = to 0 and maintained there for all t > 0. For a 

pure conduction process, differential equation (17) for each zone is 

9T 92t 

s _ s 

9t " "s 

and 


9T 92t^ 

Li L 

— — = Q/ 

9t L 9^2 

with initial conditions: 


(IS) 
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( 19 ) 


and 


If 


at 


and 


T = T = Tatt=0 
Los 


T = 0 at X = 0 (for all t) 
s 


p = p = p » then 
L s 


9T oT ,g 

K — ^ = K + p Ah — 

9x L 8x dt 


X = S(t) ; therefore, 


T = T = 

L s Fr 


T = T as X 
L o 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


The first published discussion of this classical problem was by Stefan [64l , and 
the first solution attributed to Neumann [65l . These solutions were discussed 

by Carslaw and Jaeger [65l . 


T = T, 


LIQUID 

<\. kL, Pl 


SOLID 

“s' •‘s' 

S(t) 

- . ^ * 




Figure 29, One— dimensional freezing. 
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Numerous other solutions for various boundary and initial conditions, 
for pure conduction and convection processes are available in technical publi- 
cations. Sunderland and Muehlbaur, [67], gave an excellent bibliography of 
techniques used to solve the phase change problem through 1965. Additional 
work has been done since then on specific conduction problems, some of the 
more notable being cited in References 62, 68, and 69. The more important 
convection problems are cited in References 49 and 70. The work on freezing 
and melting done under the tutelage of Golden at the Colorado School of 
mines is worthy of mention [51, 71, 72, 73, 74] . These studies attempt to 
solve the fluid motion equations in conjunction with the energy equation by 
using a numerical scheme. Although this technique is commendable, as yet, 
it has not proven to be sufficiently flexible to allow application to a broad bank 
of real physical geometries. This is because the flow pattern must be known 
before a solution can be obtained. 

The typical capacitor cell presents problems which are considerably 
more complicated than those formulated in previous studies. In a typical cell 
( I’ ig. 30) , the problem is two-dimensional, involving convection and conduc- 
tion in the one-g case, with varying boundary conditions on a minimum of three 
sides; where T^^ = T (y, t) , T^ = T(t) , and y = S(x, t) constitute the phase 
change position occurring on the moving boundary. 

Assuming again a pure conduction process, the differential equations 
for the liquid phase are: 




( 24 ) 



' ’’b = “ 

and the boundary conditions are: 
Tfi = Tg, at y = 0 ; 

= Tg, at y = S (x, t) 

and 


(25) 

(26) 

(27) 

(28) 


=0, at y = H (29) 

From Rathgen and Jiji [75] , the second boundary condition at the interface is 
given by 



(See Appendix E for the derivation of this boundary condition) . 
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Due to the complexity of this problem it has not been solved in a closed 
form. In non-idealized problems, such as this, a numerical technique is 
typicallj'^ utilized. 

Single Cell Description 

After analyzing film data (see Chapters 6 and 7) and a report from a 
previous study [76l , it was concluded that proper selection of the cells in the 
capacitor housing allow the cells to be considered thermally independent. 
However, since this adiabatic assumption could have introduced error, the cells 
were also isolated from each other by specifying fin and plate temperatures as 
boundary conditions. Computer models of the 1/4 inch and 3/4 inch cells were 
constructed. These single cell models were used to predict variations in 
temperatures and interfacial positions with time. 

The model used an explicit forward finite differencing routine. Programs 
were initially written in a simplified computer language known as CINDA 
(Chrysler Improved Numerical Differencing Analyzer [77l). Later, a more 
specific Fortran program was written which allowed operational convenienee 
in varying the parameters of the phase change study. Using identical inputs, 
the programs were compared and found to hav'C nearly identical outputs. Asa 
result, these two programs were used interchangeably to analyze test results. 

A description of these programs is given in Appendix C. 
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Initial computer runs used to provcrthe CINDA and Fortran models were 
made on an IBM 230 computer. All other runs were made on the Univac 1108 
and IBM 7094 high speed computers. Run times varied from 15 minutes to 1 
hour, depending largely on the time step required for stability and the duration 
of the experimental test. During most computer runs, a time step of 0. 0001 
hours sufficed. 

For analytical purposes the cells were subdivided into the nodal 
arrangement shown in Figure 31. The general calculation flow chart scheme 
is shown in Figure 32 and the notation used in the Fortran program is shown 
in Figure 33. 

The computational procedure started with specification of initial 
temperatures for all nodes and calculation of all physical and thermal charac- 
teristics. Typically the following’ properties were used in all runs. 

Paraffin properties: 
p =47.2 Ib/ft^ (755.7 kg/m^) 

K = 0.087 BTU/Ilr ft ^^(O.HO J/m sec °K) 

Cp = 0.5 BTU/lb°R (2092 J/kg-"C) 

AH = 73.4 BTU/lb (1.70 X lo'’ j/kg) 

AH^ = 22. 11 BTU/lb (5. 14 X lo^ j/kg) 

= 0.00045 1/R° (0.00081 l/°C) 

M = 14.3 Lbm/Hr-ft (5.9x to~^ Newton sec/m^) 
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= 73.04“F (22.81°C) 

= 89.8°F (32.13°C) 

Fr 

Metal properties: 

p =171 Ib/ft^ (2.736 X io=* kg/m=*) 

K =93 BTU /Hr ft° R (l59. 2 J/m-sec°K) 

C = 0.22 BTU/lb°R (919.9 J/kg“C) 

P 

Plexiglas properties: 

p =72.5 Ib/ft^ ( 0.94 X lo* kg/m*) 

K =0.09 BTU/Hr ft“R (0. 154 j/m-sec-“K) 

Cp = 0.33 BTU/lb°R (1380.7 j/kg-°C) 

In addition to the above property data, cell geometry data and initial 

temperature data were input. Although the model had the capability of using 

temperature dependent properties, constant values were used. An initial study 

had proven that the variation in temperature and distance outputs was not highly 

sensitive to small variations in p, C and K. Moreover, the accuracy of 

P 

available property data was questionable. Therefore, the additional computer 
time required to incorporate a varible property routine and to change a forward/ 
backward finite differencing approach was considered unnecessary. 

Node positions were located as close as possible to the fin and plate 
experimental points of measurement, and the measured values of temperature. 
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Figure 32. Computer model flow chart. 
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(i,j+1 ) 


>RV(iJ+1) 

— sAAAAA/ 

RH(iJ) 

>RV<U) 

Ofu-i) 

Figure 33. Typical node arrangement, 
as a function of time, were input for these nodes. This technique precluded 
error in the side walls, due to the adiabatic cell assumption. Due to the small 
capacitance and the resulting small computer time step required to achieve 
stability, when using a transient for analyses of the thin fins, the steady state 
Laplace equation was used to solve for the fin temperatures. Also the nodes 
were arranged in two dimensions, assuring no heat transfer down the cell. 
During most tests, this assumption was verified by the fact the interface height 
remained relatively constant across the test unit’s width. 

The heat transfer rate of each paraffin node was computed by using the 
following two dimensional numerical approximating technique: 
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( 31 ) 


Rirrini — 


^ T(i,j-1) - T(i,j) ^ T(i,j 1) - T(i,j) 
RV(i,j) HV(i,j-l) 


This rate was stored and used in subsequent calculations. The temperature of 
the node at the next time increment, t+At, was found by 

T'(i,j) = T(i,j) + (32) 

The computation was then iterated by computing q'(i,j) using equation 
( 31 ) and the new temperature found using equation (32). The average of q(i,j) 
and q’(i,j) was then used in equation ( 31 ) to predict a more accurate T'(i, j) 
value. This procedure was repeated at each time step. 

Appropriate modification to equations (3l) and (32) were made for 
nodes located on or near a boundary. Since the boundary condition at the 
interface changes due to the absorption or rejection of the latent and 
transitional heats, a special technique had to be used at the phase change and 
transitional temperatures. The method used herein forced the temperature to 
be constant until sufficient heat had been absorbed or rejected by the node to 
change phase (Fig. 34). This was accomplished by monitorir^ the total energy 
stored in each node so that when the nodal temperature reached the transition 
or fusion temperature, it was maintained at this value until sufficient energy 
was absorbed or rejected to balance the transitional or latent heat of fusion for 
the node. The fraction of the node melted or frozen was determined at each 


time step based on the energy accumulation in the node. From this, the melt 


height was determined for each column of nodes. 



TEMPERATURE 

Figure 34. Heat stored versus temperature. 


After initial runs, the basic conduction model was modified for melt 

analysis. During melting, the thermal conductivity values were modified in 

the liquid region to account for convection. The method used was to determine 

the Rayleigh number by using an average melt height, given by 
m 

^ Xn (33) 

i=l / 

The driving temperature difference was the absolute difference between the 
melting temperature and the lower plate temperature. 

at = (34) 

Using these variable quantities along with appropriate physical properties, the 
Rayleigh number was determined at each time step. 



Initially the bounded cell Nusselt number versus Rayleigh number 
correlations of Catton and Edwards [31] were considered; however, since the 
initial point at which the melt layer exceeds the critical Rayleigh number is 
extremely small, the infinite flat plate correlations of Silveston and O^Toole 
[30l were used instead. These correlations are given for the four Rayleigh 
number regimes, as 

Regime I Ra < 1700 K .VK = 1 

eff 

Regime II 3,500 ^ Ra ^ 1700 = 0.00238 Ra^'^^^ 

Regime III 1 x 10^ ^ Ra ^ 3500 "" 0.229 Ra^'252 

Regime IV Ra > 1 x lo^ ^eff"^^ " 0/l04 Ra0-305pj.0-084 

Although correlations are actually for Nusselt number, the ratio 

is the same since, K = hL. 

eff 
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Node Sensitivity Study 


To avoid excessive computer time usage without compromising 
modeling accuracy, a study was made to determine optimum node spacing. 
Although the study included both the modified and unmodified conduction models, 
no significant differences were noted. 

Half of a 1.9-cm (s/4-in.) cell was modeled. The cell was assumed 
to be isothermal along its width, so that heat was allowed to transfer only in 
the cell height and length directions. The nodes were arranged in square, 
patterns. The node sizes examined ranged from 0.16-cm (l/l6-in.) to 0.95- 
cm (3/8-in.) . The melt front position versus time data was used for 
comparison. 

Figure 37 shows the convective model runs. From this data, it can be 
seen that there was a significant difference in data output when comparing 0.48 
-cm (3/16-in.) node data with 0.32-cm (l/8-in. ) node data. However, when 
the 0. 16-cm node data was run, the interface position matched that of the 0. 32- 
cm node model, within one percent for the range of values run. Stability runs 
using the sensitivity model also showed that a stable time step of 0.001 hours 
was required to avoid unstable data outputs. 

As a result, all runs were made with nodes arranged in 0.32-cm square 
patterns, requiring 21 vertical nodes for the single cell model. In the case of 
the half model of the 1. 9-cm cell, the node matrix was 3-by-21, and for the 
0.64-cm (l/4-in.) cell it was l-by-21. 
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MELT FRONT HEIGHT -in. (cm) 





CHAPTER 6 


DATA REDUCTION 
Film Data Reduction 

All motion pictures were screened on a Data Instruments Film Analyzer. 
The apparent phase change position versus time was recorded for the center 
of paraffin cells and at the innermost fins of the three representative cells. 
Early tests demonstrated that the plate which separated the coolant fluid and 
paraffin acted as an isothermal surface. For this reason, the melting of the 
1.98-cm (3/4-in.) and 1.27 -cm (l/2-in.) paraffin cells were nearly symmet- 
ric from inlet to outlet sides of the capacitor. Some edge effects were noted 
on the outer 1.9-cm cells, so that only the visual data from the innermost 1.9- 
cm cell on the fluid outlet was used for data reduction. The innermost 1.27 -cm 
cell on the fluid outlet side of the capacitor and the center 0.r)4-cm cell were 
also used in visual data reduction (Fig. 36). 

The film analyzer projects an enlarged image of each film frame on 
the moveable grid surfaces of the film analyzer (Fig. 37) . Two grids allow 
sensiiig of both horizontal and vertical positions. These grids are electrically 
integrated with the interpreter section of the analyzer so that it automatically 
senses these positions. Each cell image is calibrated in the analyzer as the 
ratio of the known height, H, of the cell to the image height. The moveable 
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1/4 in. (0.635 cm) 1/2 in. (1.27 cm) 3/4 in. (1.9 cm) 




CELLS USED 
FOR DATA 
REDUCTION 


Figure 36. Visual data reduction cells, 
horizontal grid was positioned over the phase change vertical position and the 
height was recorded by the interpretir^ electronics. These values, along with 
corresponding timer readings from the film, were fed in by the operator 
forming a bivariant height/time array. These data were automatically key- 
punched on computer cards for later use. 

During melt tests, the phase change front position was obvious because 
of the well defined liquid/solid interface. However during freezing, a dark 
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Figure 37. Data instruments film analyzer. 
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liquid/solid dendritic front progressed into the liquid phase, followed by a 
completely solid region characterized by a chalky appearance. Since the 
dendrite front progressed well in advance of the chalky region, some confusion 
as to the actual phase charge front position resulted. Asa result, a study was 
made to determine what line defined the front. The results of the study, dis- 
cussed in Chapter 7, indicated the forward most dendrite penetration represents 
the phase change front. This definition is shown pictorially in Figure 38. 


LIQUID REGION 



Although onlj' one individual transcribed most of the film data, two 
different individuals were used to take readings on selected runs to determine 
the possible reader error. Comparisons of the independent melt height read- 
ings showed that paired data varied no more than 0. 127 cm (0.05 in.) . How- 
ever for the freeze test, the variations were as large as 0.38 cm (O. 15 in. ) . 


00 


These paired readings were compared only for melt heights above 0.51 cm 

(0.2 in.) because the epoxy used to seal the paraffin housing to the fluid 

passage partially obscured readings below this level. 

The data points on the computer cards output from the film analyzer 

were programmed, using an existing root mean square smoothing routine. 

Another routine was used to make linear approximation of the phase change 

front velocity by first fitting the smoothed height versus time data to a first 

degree polynominal, x - a + bt, and then differentiating this to get a constant 

clx 

velocity approximation ^ ^ b. This was justified by the fact that in most 
cases the heioht-versus-time data could be approximated by a straight line, 
with little error. 


Freezing Temperature Data 

The temperatvire probe profiles during all freeze runs were similar. 

A typical profile of temperature variation of the paraffin probes with time is 
shown in Figure 39 for a 1.9-cm (3/4-in.) inch cell, with numerical designa- 
tions of significant events. F\ ent number one designates the point at which the 
entire cell was at a uniform temperature. This point also coincides with inci- 
pent freezing at the cell bottom, which indicates that during freezing tests the 
initial cooling was devoted almost entirely to removing sensible heat from the 
liquid mass, cooling the entire cell to the paraffin freeze point. The tempe- 
rature did not then drop significantly until the freeze front had passed through 
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Figure 39. Typical 1.9-cm (3/4-in.) cell freezing temperature history, 
the thermocouple positions which are designated as events two, three and four. 
From visual data, it was apparent that the dendrite front passed by the thermo- 
couples just before their temperature started to drop from the plateau. The 
plateau temperature ranged from 29.5°C (85°F) to 32.3°C (90°F). In early 
tests, the plateaus occurred at temperatures near 90° F, however the plateau 
tended to occur at lower temperature during later tests. 

A typical 0.635-cm (l/4-in.) cell temperature profile during freezing 
is given in Figure 40. In contrast to the 1.9-cm (3/4-in.) cells, in which all 
thermocouples reached the freeze plateau simultaneously (Fig. 39), the 0.635- 
cm cell freeze front passed consecutively through each respective thermocouple 
location (Fig. 40) . This indicates that sensible heat was still being removed 
from the upper cell liquid while the lower cell had already frozen. This shows 
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fins to be the dominant cooling path in the small cells, whereas the lower 
plate is dominant at wider spacings. This effect was partially caused by the 
higher cooling rates due to the closer fin spacing in the smaller cells. Events 
one, three and five in Figure 39 represents the time at which the liquid had 
dropped into freeze temperature zone. Events two, four and six represents 
the freeze front passing the respective thermocouple positions. Point seven 
is at the freeze plateau level. The freeze plateau exhibited by the two cells 
were approximately equal. 



Figure 40. Typical 0.635-cm (l/4-in. ) cell freezing temperature history. 
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I nterpretation of Melting Temperature Data 
A melting history comparison is shown in Figure 41, comparing the 
characteristics of melting front progression for low and high heating rates. 
Although the 1.9-cm (3/4-in.) cell melting temperature profiles were similar, 
there were differences, depending on the heating rates. Point one corresponds 
to a temperature level at which a marked increase in the rate of temperature 
rise begins. From observations, the melt front appeared to pass the thermo- 
couple locations just after the temperature rise occurred and well before, point 
two was reached, at temperatures slightly below the published melting poi nt. 
Point two corresponds to a characteristic overshoot, followed either by a 
temperature leveling at point three or a slight drop. After maintaining a pla- 
teau (represented by point four) for a period of time, the fluid temperatures 
dropped to a level represented by point five, prior to the cell being completely 
melted. Also prior to complete melting, temperatures recorded in the liquid 
exhibited unstable oscillations between points two and six levels. At six, the 
cell was completely melted and all thermocouples rose above the pseudo pla- 
teau level. 

For the low heating rate experiments, the temperature rise which had 
occurred in the high heating rate cases was either greatly diminished or 
completely absent. The temperatures slowly rose to a plateau at point seven 
only slightly above the published melting point, with the phase change occurrii^ 
at point eight as each thermocouple reaches the plateau temperature. Again 
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Figure 41. Tj^ical 1.9-cm (3/4-iti. ) cell melting temperature history, 
after the cell was completely melted, all liquid teinperatures started to rise 
sharply, as indicated by point nine. Although unstable temperature oscillations 
were present between eight and nine, they were diminished and of a higher 
frequency than for the high heating rate cases. 

The 0.635-cm (l/4-in. ) cell data is not worthy of separate discussion 
since the primary difference as compared to the 1.9-cm (3/4-in.) cell data, 
was the delay time at the plateau temperature. In the case of the 0.635-cm cell, 
the plateau consisted only of a slight curve inflection near the values of the 1,9- 
cm (3/4-in.) cell plateau. In the fast melt cases, the 0.635-cm (l/4-in. ) 
cell melted so rapidly, that in some cases temperatures in these cells reached 
a second higher plateau corresponding to the liquid equilibrium value. 

Test Observation 

Figures 42 through 49 show views of the thermal capacitor cells 
during the phase change process for four arbitrarily selected times. Low and 
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Figure 42. Initial stages of slow melt phase change process. 
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Figure 43. Final stages of slow melt phase change process. 
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Figure 44. Initial stages of fast melt phase chaise process. 
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Figure 45. Final stages of fast melt phase change process. 
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Figure 46. Initial stages of slow freeze phase change process. 
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Figure 47. 


Final stages of slow freeze phase change 


process. 
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Figure 48. Initial stages of fast freeze phase change process. 
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Figure 49. 


Final stages of fast freeze phase change 


process. 
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and high heating rate tests (230-59 and 230-49) and low and high coolii^ rate 
tests (230-54 and 230-58) are presented, respectively. It appears from these 
photographs that the zero point on the vertical scale is below the plate level. 
However, as discussed earlier, the bottom of the cells are obscured by an 
epoxy coating used to seal the bottom plate to the coolant passage. The rear 
epoxy line appears as a lighter horizontal line in these photographs. The 
irregular dark vertical lines appearing in the 1.9-cm (3/4-in.) and 0.635-cm 
(l/4-in. ) cells are thermocouple wires of which the wider ones are the freely 
suspended sheathed wires and the narrower ones are the wires leading to plate 

and fin thermocouples. 

The absence of photographs of elevated phase change interface heights 
near the top of the cells was caused by two limitations. The first was due to 
the prohibitive length of time required to run a complete eapacitor phase change 
cycle, especially in low heat flux tests. The second, which occurred in the 
case of melting only, is the fact that when the cells had melted to within 1.27 
cm (1/2 in.) of the top, the unmelted solid tended to slide down into the liq- 
uidus region. Either of these reasons was sufficient to terminate a test. 

In the case of melting (Figures 42 through 45) some disparity is seen 
between the melt heights of the 1.9-cm (3/4-in.) cells on the right and left 
sides of the capacitor. This disparity occurred only in a few tests and was 
discernable in these photographs. The nebulous appearance is due to fur-like 
dendrite formations at the interface, shown magnified in Figure 18A. 
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Also the nearly planar front evident in the melting photographs gives 
way to a sharp saw tooth freeze front profile. The saw tooth profile is less 
pronounced for the closer fin spacing because the freeze front moving up the 
fins and the cell center are merged. 

Comparison of the meniscus formed in the upper portion of each frame, 
between the ullage and the liquid, indicates a decreasii^ level with time. This 
decrease is caused by freezing contraction occurring in the paraffin, and was 
most pronounced in slow melting tests as shown by comparison to the fast 
melting and freezir^ photographs. Since the melting heights were repeatedly 
higher on the coolant inlet side (i. e. , left side in photographs) , the bottom 
plate was apparently not exactly isothermal, but sustained a temperature 
gradient, increasing from inlet to outlet sides. 

In both rapid and slow cases it is obvious that the 1.9-cm (3/4-in.) and 
1.27 -cm (l/2-in. ) cell melting rates are nearly equal. However, a much 
faster rate is seen in the 0.635-cm cell. 

In Figures 46 through 49 (representing slow and fast freezing), the 
liquid/solid interface, which was well defined in melting, appears nebulous. 
Although not well defined in the photographs, the rough dendrite interface is 
visible as a dark region. For those frames in which the freezing of the cell is 
nearly complete, the chalky zones are visible. 

It is obvious that the freeze rates are much slower than the melt rates. 
This is caused by loss of convection in the freezing modes. 
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CHAPTER 7 


DISCUSSION OF EXPERIMENTAL AND ANALYTICAL RESULTS 
Phase Change Temperature 

Because the apparent temperature plateau achieved during melting 
tests did not correspond to the phase change temperatures exhibited during 
freezii^, only freeze tests were used to estimate variations in the phase change 
temperature. From discussions to be given later, it will be apparent that the 
phase change temperature was not constant. Figure 50 shows a plot of the 
measured freeze temperature as a function of time. It can be seen from this 
Figure that the apparent freeze temperature decreased with time. This was 
accompanied by a ''yellowing" of the solid paraffin, which had been initially 
chalky white in the solid state. 

To account for this apparent change in the freeze temperature in the 
numerical work, the phase change temperature used in the numerical compu- 
tations was altered to agree with experimental freeze data. However, since 
numerous investigators have reported that deviations between freezing and 
melting temperatures are common, referred to as phase chaise hysteresis, 
the phase change temperature was not altered for numerical predictions 
corresponding to the melt tests. However, for a single melt test, computer 
runs were made with varying phase change temperatures to illustrate the effect 
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Figure 50. Measured freeze temperature variation. 



of alterir^ this property. The influence of changing the fusion temperature 
on numerical predictions for the selected melt runs is illustrated in Figure 51. 

From this Figure it can be seen that small changes in the melting 
temperature drastically alter the predicted melt height velocity. As this 
temperature is depressed, the onset of increased interface velocity occurs 
earlier. This is explained by the fact that for the same boundary temperature 
in both cases the temperature difference in the Rayleigh Number expression 
is proportionately higher for the case with lower melting temperature. This 
reduces the height at which the onset of convection occurs and thereby causes 
the melting velocity to increase at an earlier time. No such drastic differences 
were noted for the predictions for freeze runs since contributions due to 
convection were absent and accordingly omitted from the numerical model. 

Apparent Heating/Cooling Rates 

The apparent experimental heating/cooling rates for tests were deter- 
mined by using appropriate paraffin property data along with the interface 
velocity information discussed in Chapter 6. 

Since the 1.9-cm (3/4-in.) cell freeze and melt fronts approached a 
planar front, a uniform horizontal front was assumed. Even though the front 
was not exactly planar, especially in the freezing case, this is an acceptable 
approximation since the phase change front velocity, which was used to deter- 
mine the experimental heating/cooling rates, did not alter shape with time. 
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Figure 51. Comparison of test data to predicted data. 


Also, since the sensible heat absorbed/rejected by the paraffin and its housing 
were small compared to the heat absorbed/rejected by the phase front, all of 
the energy transfer was considered to produce phase charge. Finally, to 
express the heating/cooling rates on a unit area basis, the lower plate was 
considered to be the base area through which the heat was transferred. Visual 
data revealed that the assumption of a planar interface was not valid for 0.635- 
cm cells. As a result, heating and cooling rates and fluxes are given only for 
1.9-cm (3/4-in.) cells. 

Using the approximations outlined above, the followii^ equations were 

used to generate experimentally determined heating/cooling rates, 

• ds 

Qap “ P ^ V Ah 

and 

(<5/A)ap = 

These values are given in Tables 6 and 7. 

Comparison of Experimental and Analytical Data 
The experimental data were compared to analytically generated inter- 
facial positions at varying times. Although the analytically predicted temper- 
ature profiles with respect to time and distance are discussed, it should be 
emphasized that the primary objective in formulatii^ the numerical model was 
to match the interface position data, and consequently the rate and flux data. 

It was realized that the artificial manner in which convective effects were 
incorporated into the numerical model would not provide accurate predictions 
of the temperature field in the liquid region. 
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TABLE 6. 1.9-CM (3/4-INCH) CELL COOLING RATES AND FLUXES 
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Freeze Temperature Profiles 


In the case of freezing, the analytically produced temperature profiles 
versus time and distance matched quite well with experimental results. Typi- 
cal comparisons are shown in Figures 52 and 53, 



Figure 52. Freeze test temperature versus distance. 

From Figures 52 and 53, it is apparent that the largest discrepancy is 
in the solid region. This is due to the method by which the heat transfer 
surface area was artificially increased for all nodes to account for the aug- 
mented interface velocities, thereby causit^ an apparent increase in the 
thermal conductivity of the paraffin. In the liquid region, the temperature 
profile tended to yield a better match for the experimental data with augmented 
transfer area, indicating the thermal diffusivity properties used for the liquid 
were either low or that some convection was occurring in the liquid during 
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Figure 53. Freeze test temperature versus time, 
freezing. This was possible, to a degree, since the fins were typically cooler 
than the paraffin in the center of the cell, giving rise to a doublet with flow 
from the fins toward the center, see Figure 54. 



Figure 54 


Possible doublet flow during freezing 


The temperature plateau in the liquid portion of the cell formed 
quickly and was maintained as the freeze front passed through the entire cell. 
This indicated that most sensible heat was removed from the paraffin before 
latent heat removal started. The plateau apparently was caused by the liquid 
region being surrounded by a nearly isothermal freeze boundary. Although 
not apparent on filmed data, visual inspections and the freeze models indicated 
a thin freeze film occurred on the fin and plate wall almost immediately. This 
is consistent with the rapid drop in fin and plate temperature. On reaching the 
freezit^ temperature, the fin temperatures leveled out until the freeze front 
had passed through the vicinity. 

As noted earlier, the plateau temperature decreased below the published 
phase change temperature with time. From inspection of the test data, it was 
apparent that this depressed plateau was the temperature level at which phase 
change was actually occurring. Since most freeze tests were slow, it could 
be easily determined that the liquid dendrite front had passed the thermocouple 
at about the same time the thermocouple indicated a ’’decline” from the plateau. 
Since phase change had to be occurring on the dendrite zone or the chalky zone, 
and the chalky zone passed through at temperatures well below plateau values, 
the freeze front was taken to be the line of demarcation between the liquids and 
dendrite zones and the corresponding plateau was then considered to be the 
freeze temperature. 
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Apparently the depression which occurred in the measured fusion 


temperature could be a result of contaminants ingested during testing, or al- 
ready present in the newly procured paraffin. The cause for this freezing point 
depression was not investigated. 

Melt Temperature Profiles 

The predicted and experimental results for melt tests exhibit more of 
a discrepancy in temperature versus time and distance data, see Figures 55 
and 56. The primary reason for this is the artificial way in which the effective 
thermal conductivity of the liquid was used to account for convective effects. 

The analytical predictions, based on a pure conduction model with the thermal 
conductivities modified in the liquid, do not reveal a constant temperature 
plateau. The experimental data, although exhibiting some oscillations, remains 
relatively constant at a plateau value for a period of time inversely proportional 
to the heatit^ rate. 

The plateauing is immediately preceded by a temperature jump. This 
jump may be due to the film coefficient which occurs at the liquid/solid inter- 
face as a result of convection in the liquidus region. The convection is a result 
of buoyancy induced convection in the fluid cell. In the case of melting, the 
mixing action of the convective currents caused the pseudo-plateau, as 
opposed to the isothermal boundary causing the plateau in the freezing case. 

The unsteady oscillations in the liquidus plateau region are caused by these 
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convective currents. Since these tests were run, Griggs [78l has run tests on 
individual cells which verify the occurrence of these oscillations. The plateau 
level is dependent upon the mean temperature between the melt temperature 
value at the solid/liquid interface and the plate temperature. Typically a 
simple averaging of these two values at any time approximated the plateau 
level. The plateau level is maintained in the cell until the entire cell is melted, 
which is followed by a rapid rise in all measured liquid temperatures. A 
second plateau is reached in some cases. This plateau is the limiting temper- 
ature of the cell, slightly lower than the plate temperature. 

As the heating rates were increased and for narrow fin spacings, the 
plateau became less pronounced, tending to disappear with very high heating 
rates in the 0.635-cm (l/4-in.) cells. Reference 77 verifies these findings, 
showing that the plateau duration in time and absolute level is a function of cell 
wall material as well as heating rates. 

Freeze Front Position 

Comparison of measured transient freeze front position histories with 
corresponding numerical predictions are shown in Figures 57 and 58 for slow 
and fast cooling rates, respectively. From these figures it is apparent that a 
pure conduction modeling technique under-predicts the transient freeze inter- 
face position as well as the freeze front velocity. 

A literature research reveals this disparity has also been noted in 
References 8 and 74. Reference 8 gave no explanation for this underprediction. 


118 



CELL fl PHASE CHANGE HEIGHT - In. (cm) 



Figure 57. Comparison of freeze front position test data and predicted 
data for a slow cooling rate. 

Although Golden [74] attributed this disparity to erroneous property data, evi- 
dence shows this not to be the case in this instance. 

Possible causes of underprediction are liquid phase convection, aniso- 
tropic or erroneous property data resulting in higher thermal diffusivities or 
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Figure 58. Comparison of freeze front position test data and predicted 
data for a fast cooling rate. 

latent heat values than used in the model, and augmented heat transfer surface 
area at the interface due to dendrite formation. Bailey* has shown however, 


^Verbal communication between the author and Dr. Bailey of N.C. 
State University. 
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that the deviation in latent heat of fusion of the paraffin used in these tests 
from the values quoted in the literature is negligible. 

Convection in the bulk liquid or local interdendrite convection could 
produce this underprediction. However, since freezing is from below the 
former is not expected to be a large effect, and the latter effect produced by 
the volume change phenomena has been shown in Chapter 2 to be negligible. 

To examine the approximate level of paraffin thermal conductivity 
change required to match test data, a highly simplified technique was used. 

The vertical conductors (i.e. KA/AX values) in the conduction model were 
artifically increased by applying multiplying factors to the conductors in the 
liquid and solid phases. These factors were varied until a good interface 
position match was attained. 

Plots of fast and slow cooling rate tests showing the effect of the opti- 
mum match are shown in Figures 59 and 60, respectively. Applying this tech- 
nique to a number of tests, multiplying factors ranging from 4 to 10 resulted. 
Although the technique used to reproduce test data is highly simplified, it 
serves to show the magnitude of thermal conductivity coefficient augmentation 
required is much greater than any expected increase due to anisotropic effects 
or property measurement inaccuracies. 

Using the author’s data, Griggs [77] has shown that by altering the 
vertical conductors at the interface only, good approximations of test results 
can also be achieved (Fig. 61). Extrapolation of these data indicates multi- 
plication factors on the order of 100 to 200 are required for best results. 
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Examination of microscopic photographs indicates these magnitudes are not 
inconsistent with the dendrite geometry at the interface. 



Figure 59. Comparison of freeze front position test data and predicted 
data for fast cooling rate. 

From the foregoing discussion it can be concluded the exclusion of the 
heat transfer surface area augmentation at the interface is the most probable 
cause of underprediction. However, combination effects of augmentation with 
other effects discussed is also a possibility. 
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Figure 61. Comparisons of predicted interfacial positions at center of 
cell for MSFC in-house run 230-7. 

Melt Front Position 

The pure conduction model did not predict the melt front position with 
any degree of accuracy (Fig. 62). Estimations of natural convective levels, 
clearly indicated that this was due to the augmentation created by the convective 
currents. As a result, the basic conduction program was modified, as dis- 
cussed earlier, to provide an estimate of convective enhancement. An esti- 
mating approach was taken rather than try to solve the convection problem, 
because of the complicated nature of the problem when the energy is coupled 
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with the momentum equations. No attempt was made to form convection co- 
efficient conductors because insufficient film coefficient data were available 
to estimate these coefficients at a melting solid/liquid interface. 



Figure 62. Conduction mode sensitivity for 1.9-cm (3/4-in.) cell. 

Comparison plots of experimental data with analytical data are shown 
in Figures 63 through 71 for 1.9-cm (3/4-in.) cells and in Figures 72 through 
76 for 0.635-cm (l/4-in.) cells. From these plots, it is obvious that in most 
cases the modified model satisfactorily represents melting rates and conse- 
quently the heating rates. The analytical data in most cases slightly under- 
predicts the experimental results, although the lag is only slight. 
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Figure 63. A comparison of experimental and analytical data for test 230-5. 

A constant melt temperature of 89. 8T was used in the model for all 
computer sums. The experimental results of freezing tests indicated that phase 
change temperature appeared to decrease with time. The later predictions 
could have been improved by decreasing the model melt temperature to the 
experimental level. 
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Figure 65. A comparison of experimental and analytical data for test 230-8. 
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Figure 66. A comparison of experimental and analytical data for test 230-9. 
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Figure 67. A comparison of experimental and analytical data for test 230-49. 
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Figure 68. 


A comparison of experimental and analytical data for test 230-50. 
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Figure 69. A comparison of experimental and analytical data for test 230-57. 
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Figure 70. A comparison of experimental and analytical data for test 230-59. 


33 





CELL PHASE CHANGE HEIGHT - in. (cm) 



1/4 INCH (0.635 cm) CELL 
TEST 230-5 


O 


O TEST DATA ^ 

— CONVECTION MODEL 



TIME IN SECONDS 


Figure 72. A comparison of experimental and analytical data for test 230-5 
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Figure 73. 


A comparison of experimental and analytical data for test 230-57. 
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Figure 74. A comparison of experimental and analytical data for test 230-59 
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Figure 75. A comparison of experimental and analytical data for test 230-8. 
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Figure 76. A comparison of experimental and analytical data for test 230-56. 
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From the data is it obvious that the phase change position versus time 


for melting runs can be approximated by a straight line. This phenomena can 
be justified analytically. In Appendix B the Rayleigh Number is given as 
Ra = 1.68 X 10^ AT 

For these tests the range of temperature differences in the liquid region was 
2.8°C (5°F) < AT < 13.9°C (25°F) 

Using these limits with the above equation and the Rayleigh Number 
regimes defined by Silveston and O'Toole [30l , a plot of Rayleigh Number versus 
melt height may be generated. (See Figure 77). 

From this figure, it is apparent that most of the cell is operating in 
regime 4. Although this is basically a transient problem, a steady-state 
approximation of the velocity profile may be made. Noting that the Nusselt 
Number is given by 


and combining this with the regime 4 Nusselt Number expression, the heat 
transfer coefficient may be approximated by 


. K „ 0.305 ^ 0.084 

h = 7 ( 0 . 104) Ra Pr 


Again assuming constant properties and a constant temperature differ- 
ence, the following relation holds 


h 


Cl L 


-0.085 
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neglecting sensible heat and an energy balance on the melt interface yields 


hAT = AHpV 


or 


V 


CzL 


- 0. 085 


This indicates that velocity is relatively constant, only weakly dependent on the 


melt height. 


142 



CHAPTER 8 


SUMMARY AND CONCLUSIONS 

1. In a one-g environment, a pure conduction model cannot be used to 
determine freeze or melt front position accurately in paraffin. The freeze 
and melt front positions are both underpredicted. In the case of melting, the 
actual rate is faster because of the augmenting effect of natural convection. 

For freezing, the inability to reproduce experimental results is thought to be 
due to the influence of dendrites. Howev'er, this suggestion is not proven herein. 

2. A convection model accurately predicted experimental results for 
melting. However, a notable fallacy is that this model cannot accurately re- 
produce temperature profiles in the liquid phase. This model utilizes an 
effective thermal conductivity variation with Rayleigh Number to simulate 
conductive effects. Although the model has been demonstrated to give good 
results for only the geometry and heat flux ranges discussed in the text, the 
model is expected to apply to larger cell sizes and for lower and higher heat 
fluxes up to nucleate boiling initiation. Because of the reversal trend predicted 
by Jones and Smith [79l , for cell sizes below 0.635-cm (l/4-in.) (Figure 78) 
alterations are necessary. These alterations require the Silveston and O'Toole 
correlations be replaced by the Catton and Edwards closed cell correlations. 


143 





FIN SPACING, S, IN. 


Figure 78. Typical effects on fin spacing on average heat-transfer coefficients. 

3. During freezing of paraffins the phase change boundary is made up of 
numerous dendritic arms. This phenomena is not present in melting. Although 
this effect is well known in material sciences, it was not considered by a 
number of authors investigating heat transfer during solidification. 

4. The freeze front position cannot be matched without making unverified 
assumptions. It is apparent that dendrite formations augment the surface area 
available for heat transfer at the interface during freezir^. Augmentation of 
this area will improve conduction model data matching. It is hypothesized that 
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failure to include this augmented area effect is the major cause of the reported 
underprediction. Further study of these effects is needed before positive 
conclusions may be reached. 

5. Insufficient information is available to determine if the unexpectedly 
rapid freeze front progression noted in one-g will also occur in zero-g. How- 
ever, assuming the augmented rate is due to dendrite formations, as hypoth- 
esized, the paraffin freeze rate should again be underpredicted by a pure con- 
duction model. This is true since there is no evidence in the literature that 
dendrite formations are affected by the gravity level. 

6. In a zero-g environment, conduction modeling techniques are expected 
to be applicable to determining melt front positions and temperature profiles 
unless significant surface tension driven convection is present. However, 
there are insufficient zero-g data to corroborate this finding. 

7. From the literature, surface tension driven convection can be appreci- 
able for normal paraffin capacitors in zero-g only in systems contain- 
ing adjacent liquid/gas phases. Insufficient data are currently available to 
allow quantitative estimates of convective levels created by this phenomena. 

8. Excluding convection due to orbital maneuvers, the only other affect 
expected to be important during phase change is that of ullage bubble location. 
The ullage location, in the case of paraffins, will be favorable for maximum 
heat transfer from metal surface, because of good wetting characteristics of 
paraffins. 
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APPENDIX A 


PARAFFIN PROPERTIES 
Introduction 

A tabulation of paraffin property data is given in Table A-1. These 
data were compiled from the authors listed below in references A-1 through 
A -10. Items A -25 through A -59 give a list of property references compiled 
by the Thermophysical Properties Research Center at Purdue University. 
These authors in-turn referenced those items listed for A -11 through A -24. 

In most cases the authors did not specify the paraffin grade, however, it is 
assumed that all data was acquired by using pure or research grade paraffin. 
Only those data considered reliable are presented herein. 

Data for twenty normal paraffins with carbon atom chains ranging from 
eleven to thirty are presented. These paraffins cover a freezing point temper- 
ature range from -25.6°C (-14°F) to 65.6°C (150°F), encompassing the 
normal interest of the thermal capacitor designer. 

The reader is referred to references A-4, and A-5, for more specific 
information of property variations with temperature. 

Available paraffin cost data is given in Table A -2. 
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7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Name 
Undecane 
Dodecane 
Tridecane 
Tetradecane 
Penta decane 
Hexadecane 
Hepta decane 
Octadecane 
Nona decane 
Eicosane 
Heneicosane 
Docosane 
Tricosane 
Tetracosane 
Pentacosane 
Hexacosane 
Heptacosane 
Octacosane 
Nona cos ane 
Triacontane 


Chemical Molecular 
Formula Weight 


C 11 H 24 

156.302 

C 12^26 

170.328 

C 13^28 

184.354 

C 14 H 30 

178.380 

C 15^32 

212.406 

C^16^34 

226.432 

C 17^36 

240.458 

00 

U 

254.484 

C 19 H 40 

268.510 

C^20^42 

282.536 

C 21 H 44 

296.562 

^22^46 

310.588 

^23^48 

324.614 

^24^50 

338.640 

^25^52 

352.666 

^26 ^54 

366.692 

^27^56 

380.713 

^28^58 

394.744 

^29^0 

408.770 

^30^2 

422.796 


Phase Change 
Temperature 


®F 

®C 

-14.1 

-25.6 

14.7 

-9.6 

22.3 

-5.4 

42.6 

5.9 

49.9 

10.0 

64.7 

18.2 

71.6 

22.0 

82.8 

28.2 

89.4 

32.1 

98.2 

36.8 

104.9 

40.5 

111.9 

44.4 

117.7 

47.6 

123.6 

50.9 

128.7 

53.8 

133.5 

56.4 

138.2 

59.0 

142.5 

61.4 

146.7 

63.8 

150.4 

65.8 


Transition 

Temperature 


°F 

®C 

-32.8 

-36.0 

none 

none 

-0.4 

-18.0 

none 

none 

27.9 

-2.3 

none 

none 

50.9 

10.5 

none 

none 

73.0 

22.8 

97.2 

36.3 

90.5 

32.5 

109.4 

43.0 

104.9 

40.5 

118.6 

48.1 

117.6 

47.6 

127.9 

53.3 

127.4 

53.0 

137.4 

58.6 

136.8 

58.3 

1 143.6 

62.0 


^Values taken at 158® F (70® C) 
**Values taken at lOTT (41.7°C) 











Density at 

Change Temperature 


Latent Heat 

1 Heat of 

Total Heat 



of F 

'usion 

Tran 

sition 


Solid 

Btu/lb 

j/g 

Btu/lh 

> j/g 

Btu/lb 

J/g 

Ib/ft^ 

Kg/m3 

Ib/ft^ 

Kg/m 

57.9 

134.6 

18.4 

42.7 

76.3 

177.3 


- 

- 

- 

93.0 

216.1 

none 

none 

93.0 

216.1 

- 

_ 

- 

- 

66.5 

154.5 

17.9 

41.6 

84.4 

196.1 


- 

- 

- 

97.7 

227.1 

none 

none 

97.7 

227.1 

48.4 

774.4 

50.8 

812.8 

70.5 

163.8 

18.6 

43.2 

89.1 

207.0 

- 

- 

- 

- 

101.3 

235.4 

none 

none 

1 

101.3 

235.4 

48.4 

774.4 

52.1 

833.6 

72.4 

168.3 

19.6 

45.6 

92.0 

1 

213.9 


- 

- 

- 

104.7 

243.3 

none 

none 

104.7 

243.3 

48.5 

776.0 

53.3 

852.8 

73.4 

170.6 

22.1 

51.4 

95.5 

222.0 

48.2 

771.2 

52.4 

838.4 

107.0 

248.7 

none 

none 

107.0 

248.7 

48.6 

777.6 

51.2 

819.2 

69.1 

160.6 

22.4 

52.1 

! 

91.5 

212.7 

47.4* 

758.4* 

- 

- 

69.1 

160.6 

39.6 

92.0 

108.7 

252.6 

47.6* 

761.6* 

- 

- 

71.5 

166.2 

28.8 

66.9 

100.3 

233.1 

47.7* 

763.2* 

- 

- 

69.7 

162.0 

39.7 

92.3 

109.4 

254.3 

48.02* 

768.3* 

- 

- 

76.7 

178.3 

31.8 

73.9 

108.5 

252.2 

48.01* 

768.3* 


- 

69.8 

162.2 

37.8 

87,8 

107.6 

250.0 

48.2 

771.2* 

- 

~ 

68.2 

158.5 

32.7 

76.0 

100.9 

234.5 

48.7 

779.2 

- 

- 

70.4 

162.8 

38.6 

89.7 

109.0 

252.5 

48.4* 

774.4* 

- 

- 

69.5 

161.5 

31.2 

72.5 

100.7 

234.0 

- 

- 

- 

- 

108.0 

251.0 

none 

none 

108.0 

251.0 

- 

- 

- 

- 


Specific Heat 
at Phase Change 
Temperature 


Btu/lb-°R 


J/Kg-"K 


0,50 


0.51 


0.52 


0.53 


2092 


2134 


2176 


2218 


0.52 


2176 
















TABLE A-1. PARAFFIN PROPERTY DATA 


Coefficient 
Conductivi 
Change T< 

of Thermal 
ty at Phase 
3mperature 

Boiling Point at 
14.7 psia 

Heat of Vaporiza- 
tion at 14.7 psia 
and Boiling Point 

Vapor Pressure 

Coefficient of 
Expansion at 
60“F (15.6°C) and 
14.7 psia 

(1.013 X 10^ N/m^) 

Btu/hr-ft-°F 

j/sec -m-°K 


'^C 

Btu/lb 

j/g 

Psia 

@°F 

©-c 

l/°F 

l/'C 

- 


384.60 

196.05 

114.2 

265.4 

41.6 

220 

104.5 

0.00056 

0.00101 

- 

- 

421.30 

216.45 

110.2 

256.1 

- 

- 


0.00055 

0.00099 

- 

- 

455.79 

235.63 

106.3 

247.0 

41.5 

139 

59.5 

0.00052 

0.00094 

0.087 

0.149 

483.43 

251.00 

103.4 

240.3 

1.03 

166 

74.5 

0.00051 

0.00092 

- 

- 

519.13 

270.84 

100.0 

232.4 

- 


- 

0.00050 

0.00090 

0,087 

0.149 

548.23 

287.07 

97.7 

227.0 

1.03 

222 

105.6 


- 

- 


575.28 

302.06 

95.0 

220.8 

- 

- 

- 

- 

- 

0.087 

0.149 

601.02 

316.38/ 

92.6 

215.2 

1.03 

347 

175.1 

- 


- 

- 

625.5 

329.99 

90.5 

210.3 

- 

- 

- 


- 

0.087 

0.149 

648.9 

343.00 

87.8 

204.0 

1.03 

388 

197.9 

0.00047 

0.00085 

- 

- 

671.2 

355.39 

- 

- 

- 

- 

- 

- 

- 

- 

- 

692.6 

367.29 


- 


- 

- 

- 

- 

- 

- 

712.9 

378.58 

- 

- 

- 

- 

- 

- 

- 

- 

- 

732.6 

369.53 

- 

- 

- 

- 

- 


- 


- 

751.5 

400.04 

- 

- 

- 

- 

- 

- 

- 

- 


769.5 

410.05 

- 

- 

- 

- 


- 

- 

- 

- 

786.9 

418.61 

- 

- 

- 

- 

- 

- 

- 

- 

- 

803.7 

429.07 

- 

- 

- 

- 

- 

- 

- 

- 

- 

819.9 

439.07 

- 

- 

- 

- 

- 

- 

- 

' 

- 

835.5 

446.74 

- 

- 

- 

- 

- 

- 

- 
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Surface Absolute Viscosity 

Tension of Liquid at 68“ F 

at 14.7 psia (20“C) and 14.7 psia 

(1.013 X iq 5 (1.013 x lo^ N/m^) 


dynes/cm Centipoise N-sec/m^ 


1.185 0.001185 

1.503 0.001503 

1.880 0.00188 

2.335 0.002335 

2.863 0.002863 

3.474 0.003474 

4.196 0.004196 


4.29** I 0.00429** 


Refractive Index 
of Liquid at 68° F 
(20“C) and 14.7 psia 
(1.013X iqS N/m2) 


1.4173 


1.4216 


1.4756 


1.4289 


1.4355 


1.4368 


1.4389 


1.4408 


1.4425 


1.4420 


1.4247 


1.4260 


1.4276 


1.4286 


1.4302 


1.4310 


1.4321 


1.4330 


1.4340 


Heat of Formation 
of Gas at “R 


Btu/lb-mole J/gm-mole 


-88434 -2. 05 X lo® 

-94986 -2.20 x lo^ 

-101556 -2.36 X lo^ 

-108176 -2.51 X 105 

-114678 -2.67 x iqS 

-121230 -2.82 X 105 

-127782 -2. 97 x io5 

-134352 -3. 12 X lo® 

-140886 -3.27 x io5 

-147456 -3.43 x io5 


1.4348 








Net Heat of Combustion at 
77®F (25°C) and 14.7 psia 
(1.013 X 10^ N/m^) 


Btu/lb-mole j/gm-mole 


2968.2 


3230.3 


3492.6 


3754.7 


4016.8 


4278.9 


4541.0 


4803.5 


5065.4 


5327.8 


10553.3 


11163.3 


11265.5 


CRITICAL CONSTANTS 

Flash Point Pressure [remperature \ Specific Volume 

^C 1 ftVlb I w?/Kg 


Psia 1 N/m^ 


6898.1 149 65 

7507.2 160 71 

8116.8 175 79 


282 1.94 X 10^ 6 94 36 8 0.06 7 6 0.0042 

262 1.81 X 10^ 728 387 0.0676 0.0042 

250 1.72 X 10^ 76 1 4 0 5 0 . 06 6 7 0.0042 


8725.9 250 121 235 1.62 X 10® 791 422 0.0667 0.0042 

9335.8 - - 220 1.52 x 10® 818 417 0.0667 0.0042 

9944.2 - - 206 1.42 X 10® 844 451 0.0667 0.0042 

10553.3 - - 1.32X10® 869 465 0.0667 0.0042 

176 1.21 X 10® 890 477 0.0667 0.0042 

162 1.12X10® 912 489 0.0667 0.0042 

162 1.12X10® 932 500 0.0667 0.0042 


12381.8 















A-1 

A -2 
A -3 

A -4 

A -5 

A -6 

A -7 
A -8 

A -9 

A-10 

A-11 

A-12 
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TABLE A-2. PARAFFIN COST DATA* 



Technical 

Grade 

Pure 

Grade 

Research 
Grade * 



Cost/lb 

Cost/lb 

Cost/lb 

Cost/lb 


Paraffin 

> 90% Pure 

> 95% Pure 

> 99% Pure 

> 99.8% Pure 

1 

Undecane 

_ 

$ 6.40 

$11.30 

$ 86. 40 

2 

Dodecane 

- 

5. 75 

.08 

107. 30 

3 

Tridecane 

- 

10. 10 

20.00 

106. 40 

4 

Tetradecane 

- 

8.60 

9. 20 

- 

5 

Pentadecane 

- 

13.60 

20.00 

- 

6 

Hexadecane 


14.30 

8.40 

- 

7 

Heptadecane 


13.60 

22.00 

- 

8 

Octadecane 

- 

13.60 

9.30 

- 

9 

Nonadecane 

- 

12,70 

24.00 

- 

10 

Eicosane 

$13.70 

1. 60 

1 . 30 

- 

11 

Heneicosane 

- 

- 

- 

- 

12 

Docosane 

- 

- 

30.00 

- 

13 

Tricosane i 

- 

- 

- 

- 

14 

Tetracosane 

- 


30.00 

- 

15 

Pentacosane 

- 

- 

- 

- 

16 

Hexacosane 

- 

- 

- 


17 

Heptacosane 

- 

- 

- 

- 

18 

Octacosane 

' - 

- 

30. 00 

- 

19 

Nonacosane 

- 

- 

- 

- 

20 

Triacontane 

- 

- 

- 

- 


* Lowest 1972 prices from Eastman (Kodak) Organic Chemicals, Rochester, 

N. Y, ; Humphrey Chemicals, North Haven, Conn. : and Phillips Petroleum 
Co. , Bartlesville, Okla. - based on quoted quantity prices nearest one pound. 

> 99. 9 percent purity sold for 90 per 5 ML sample. 
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APPENDIX B 


NON-DIMENSIONAL NUMBERS 
Introduction 

When estimating thermal performance in zero-g and one-g environments, 
the Rayleigh Number, Marangoni Number and Bond Numbers often appear. The 
Rayleigh Number can be used to estimate the convective level due to natural 
buoyancy in a one-g field. The Maragoni Number may be used to correlate the 
Nusselt Number at zero-g. Finally, the Bond Number may be used to assess 
the relative magnitude betv»een buoyancy driven and surface tension or cap- 
illary driven convection. 

Rayleigh Number 
The Rayleigh Number is given by 
Pa = GrPr = ^ (Tg- ^ . 

where 

_ k 

pep ’ 

g = gravitational constant = 32.2 ft/sec^ (9.8 m/sec), 

/3 = coefficient of Volumetric expansion, 

L = characteristic fluid dimension — in this case L is taken to be 
the average height of the fluid phase in a cell, 
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Tg = boundary temperature of the hot fluid temperature at the transfer 
plate, 

= is the equilibrium temperature at the phase front = 

fJL 

V = — the kinetic viscosity, 

k = the coefficient of thermal conductivity in the liquid phase, 
p = the average liquid phase density, 

Cp = the specific heat of the liquid phase, and 
^ = the absolute viscosity of the liquid phase. 

Assuming the following constant values for these properties, the one-g 
Rayleigh number for melting may be found in terms of the fluid temperature 
difference and characteristic dimension: 

)3 = 0.00045 (l/°R) [0.00081 (1/°C)] , 

k = 0.087 BTU/HR-FT-°R (0. 149 J/sec-m-“K), 
p =47.2 Ib/ft^ (755.2 kg/m^), 

Cp = 0.5 BTU/lb m-°F (2092 J/kg-°K), 
p = 14.3 Ibm/hr-ft ( 0. 0059 n j/m^) , and 
givir^ for one-g 

Ra = 1.68 X 10® tJ AT; 

where 

AT= T - T 
S F 

For a zero-g determination, the Rayleigh Number is zero because of the 
gravity term. For the case where the net equivalent gravity acting on the 
melting cell is reduced to 1 x io“^ g, the zero-g Rayleigh Number is given by 
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Ra 


16.8 L3 at 



Substituting a surface tension value of 28.7 dynes/cm at 20‘^C (68°F) , a critical 
temperature of 490"C, a density of 0.79 gm/cm'^ and a molecular weight of 
268.51 for nonadecane [ 14] , the value of a is found to be 3.0. 

Now using equation (B-3) , we have 
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da 

dT 


-0.062 


dynes 

cm-“C 




SO that the Marangoni Number may now be determined by using 


(B-4) 



p va 


(B-5) 


Assuming a linear temperature profile, the — , l} product becomes ATL. 

Q J_j 

Substituting the liquid properties for nonadecane given earlier, the gravity 
independent Marangoni Number relationship is found to be 

Ma = 1684 ATL (B-6a) 

for At in “C and L in centimeters, and 

Ma = 2374 ATL (B-6b) 

for AT in °F and L in inches. 


Bond Number 


The Bond Number is given by 


Bo = — — 
a 


(B-7) 


Again using the nonadecane liquid properties given earlier, the one-g Bond 
Number is given by 

Bo = 26.6 L^ (B-8a) 

where L is in centimeters, and 

Bo = 171.6 l 2 (B-8b) 

where L is in inches. 
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Reynolds, et al. , [81] have shown by a graphic representation similar to Figure 
B-1, the hydrostatic regimes for typical liquids. The gravity dominated 
regime being above the paraffin diagonal in this Figure and the capillary or 
surface tension dominated regime being below. Operatir^ points falling in a 
near proximity of the diagonal are not clearly dominated by either phenomena. 



Figure B-1. Hydrostatic regimes for typical liquids. 
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APPENDIX C 


COMPUTER PROGRAMS 
Introduction 

The computer programs used in the foregoii^ study are detailed in this 
Appendix. Some of the information given herein is excerpted from Reference 
78. Copies of the FORTRAN programs for a typical melting run with inclusion 
of convective effects (Table C-l), a typical freezing run, written in FORTRAN 
(Table C-3) , and a CINDA freezing run (Table C-4) are given, along with the 
corresponding notations. The details of each step is discussed for the melting 
FORTRAN model in Table C-2. All programs, excluding the CINDA program, 
are written in FORTRAN V. Programs utilize explicit forward finite differ- 
encing techniques. 

Since the CINDA model uses techniques similar to those employed in 
the FORTRAN model no explanation is given and the reader is referred to 
Reference 77 for explanations of listed sub routines. A skeleton flow chart for 
the melting model is given in Figure C-l. 

It was later discovered that the FORTRAN program for melting would 
not run for the case of m = 1. The statements causing this incompatibility are 
indicated in the melting pregram listing by an arrow placed at the left of the 
appropriate statement. The FORTRAN freezing program incorporates the 
changes necessary to allow runs for m = 1. 
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A J 

AM 

B 

BETA 

BJ 

C(I,J) 

CL 

CP 

CPI 

CP2 

CP3 

DAVG 

DDOT 

CELT 

DEN 

DENI 

DEN2 

DEN3 

DFLO* 

DFIN 

DIF 

DMID 


NOMENCLATURE 


1-3, used in computing interface location 

value of integer M converted to floating point 

length of section, ft, 

volume expansivity of wax, R"’ 

J-2, used in computing interface location 

thermal capacitance of node Btu/F 

constant pressure specific heat of wax, Btu/lbmF 

constant pressure specific heat of wax, Btu/lbmF 

constant pressure specific heat of bottom plate, Btu/lbmF 

constant pressure specific heat of fin, Btu/lbmF 

constant pressure specific heat of top plate, Btu/lbmF 

average height of liquid based on amount melted, in. 

interfacial velocity for nodes adjacent to centerline (I=MM), in/hr. 

absolute value of temperature difference betv/een bottom plate 
and interface, F. k 

wax density, Ibm/ft^ 


bottom plate density, Ibm/ft 
fin density, Ibm/ft^ 


top plate density, Ibm/ft^ 


approximate interfacial location for nodes adjacent to fin 
based on amount melted being equal to F2JM, in. 


(1=2) 


interfacial position for nodes 
amount being melted, in. 


adjacent to fin (1=2) based on any 


temperature difference used 
during iteration when sol 
fin temperatures , °F. 


in comparing new and old temperatures 
ving steady state equations for unspecified 


interfacial position for node adjacent to centerline, in. 
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NOMENCLATURE (Continued) 


DM0 

DT 

EPS 

ERROR 

F(1,J) 

FAC 

F2JM 

G 

H 

HBOT" 

HMELT 

HTOP^ 

HTR 

I 

J 

JOE 

KCHK 

KCOUNT 

M 

MCOUNT 

MFIN 

MM 


DMID evaluated at previous time, in. 
time increiDent, hr. 

arbitrarily set small number used as a comparator 

percent error in computed energy balance based on transfer 
rates, percent. 

mass fraction of node (I.J) which has undergone phase change 
since start of process 

time ratio used in linearly interpolating specified fin temperatures 
at a particular time in terms of bracketed data values. 

fraction of S which corresponds to 1/32 inch (arbitrary) 

acceleration of gravity, ft/hr^ 

PCM section height (See Figure 31), ft. 

heat transfer coefficient between external fluid and bottom plate, 
Btu/hr-ft^-F 

heat of fusion, Btu/lbm 

heat transfer coefficient between external fluid and top plate, 
Btu/hr-ft^-F 

heat of transition, Btu/lbm. 

integer designation of vertical column in which a mode is located 
(See Figure 31) 

integer designation of horizontal row in which a mode is located 
(See Figure 31) 

counter used in refining the heat transfer computation before 
progressing in time 

integer used to control printing of results at desired times 
(See definition of KCOUNT) 

integer counter used to print our results at times when 
KCOUNT = KCHK 

number of wax nodes in a horizontal row 

counter used in determining unspecified fin temperatures 

maximum value of MCOUNT which when exceeded causes program to stop 

M+1 (See Figure 31) 
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NOMENCLATURE (Continued) 


n 

ii) 

NDP 

NI 

NJ 

NN 

PR 

Q(I.J) 

QBT 

QBW 

QFTR 

QMELT 

QRAT(I,J) 

QS(I,J) 

QSIN 


number of wax nodes in a vertical column 

number of data points for measured fin temperatures 

ND-1 

N+4 

N+3 

N+2 

Prandtl number 

unnecessary variable - replaced where needed by QS(I,J) 

Instantaneous heat transfer rate through bottom, Btu/hr 

instantaneous heat transfer rate through bottom to wax only, 
Btu/hr ^ 

instantaneous rate of heat transfer to fin, Btu/hr. 

energy which acounts for amount of wax melted at any time, Btu. 

the instantaneous net rate of heat transfer to node I,J, Btu/hr. 

the energy stored by node I , J above TREF for wax and above 0 for 
metal nodes, Btu 

the energy stored by node I,J above TREF correspond! ng to initial 
temperature throughout network, Btu. 


QSW 

QTOP 

QTTR 

QTW 

QWAX 

Q1 

Q2 

Q3 

Q4 


instantaneous heat transfer rate from fin to wax, Btu/hr. 

instantaneous heat transfer rate out of top of section, Btu/hr. 

instantaneous rate of heat transfer to top plate, Btu/hr. 

instantaneous heat transfer rate from top plate to wax, Btu/hr 

net energy transfer to wax since start, Btu. 

energy stored by wax node above TREF corresponding to start 
of phase transition, Btu. 

energy stored by wax node above TREF corresponding to end of 
phase transition, Btu. 

energy stored by wax node above TREF corresponding to start 
of melting, Btu. 

energy stored by wax node above TREF correspond! ng to end of 
melting, Btu. 
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PA 

KATIO 

RH(1 ,J) 

RV(I,J) 

S 

51 

52 

53 

TAMB 

TAU 

TIME 

TIN 

TK 

TKL 

TKR 

TMELT 

TREF 

TTR 

Tl(I.J) 

T2(I,J) 

T3{1,J) 

TAU2 


NOMENCLATURE (Continued) 

ins tcinianeous rate of heat transfer to node I ,J based on temperatures 
obtained from QRAT(I,J) and then used to correct temperature 
predi ctions , Btu/hr 

Rayleigh number for liquid wax 

ratio of instantaneous heat transfer rate from fin to wax to that 
from bottom plate to wax 

horizontal thermal resistance between node I-1,J and node I ,J , 
hr-F/Btu 

vertical thermal resistance between node I ,J and I,J-1, Hr-F/Btu 

wax node width , ft. 

bottom plate thickness, ft. 

fin thickness , ft. 

top plate thickness, ft. 

temperature of environment external to top plate, F 
limiting time value to stop program, Hr. 
instantaneous value of time. Hr. 
initial temperature of all nodes, °F 

wax thermal conductivity (artificially allowed to vary in liquid 
to account for convection), Btu/hr- ft-F 

thermal conductivity of liquid, Btu/hr-ft-F 

ratio of effective thermal conductivity to thermal conductivity 
fusion temperature of wax, °F 

arbitrary reference temperature (should be less than TIN), F 
transition temperature, °F 
temperature of node I ,J at time t, ®F 
temperature of node I,J at time t+At, ®F 

temperature of node 1 ,J in fin at beginning of^each iteration 
step used in finding steady state solution, F 

arbitrarily defined time value used in print-out control. Hr. 
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NOMENCLATURE (Concluded) 


TKl 

TK2 

TK3 

TIM(L) 

TMl(L) 

TM2(L) 

TM3(L) 

TM4(L) 

VIS 

VI 

V2 

W 


thermal conductivity of bottom plate. Btu/hr-ft-F 
thermal conductivity of fin, Btu/hr-ft-F 
thermal conductivity of top plate, Btu/hr-ft-F 

time value corresponding to input data of measured fin temperatures. 

measured bottom plate temperature (input data), ®F 
first measured fin temperature (input data), °F 

second measured fin temperature (input data), °F 

third measured fin temperature (input data), “F 

viscosity of liquid, Ibm/hr-ft 
volume of wax melted at time t, ft^ 
volume of wax melted at time t+At, ft^ 
width of wax cell , ft 


1 This assumes that some finite thickness must have melted before it 
would be detectable on the film. The nunter DFIN is the height 
correspondi ng to a node with any amount melted. 

2 This was included to be general but has not been used to date as 
bottom plate temperatures were specified as input data. 

This has been included but set at a small value to essentially 
correspond to the top being insulated. 
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TABLE C-1. FORTRAN COMPUTER PROGRAM FOR MELTING 


1 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


DIMENSION RV(5,31) , RH(5,31) , C(5,31) , (JS(5,31 
DIMENSION T2(5,31) , F(5,31) , Q(5,31) , QRAT(5,31 
DIMENSION TM1(9), TM2(9) , TM3{9) , TM4(9) , TIM(9) 
DIMENSION QRA2(5, 31) 

COMPUTATIONAL PARAMETERS 
N=27 
M=4 
AM=M 
MM=M+1 
NN=N+2 
NJ=N+3 
NI=N+4 
ND=9 
riDP = 8 
TAU=1 . 1 
KC0UNT=1 
MC0UNT=1 
MFIN=500 
EPS=1 .E-06 
KCHK=1000 
J0E = 1 

DT=1 .OE-04 

TAU2=(10.*DT)+(DT/3.) 

PHYSICAL PROPERTIES 
TAMB=80 . 

TIN=73.5 
HT0P=1 .E-08 
HB0T=5.0 

G=( 32. 2*3600. *3600.) 

WAX 

DEN=47.2 

TK=0.087 

CP=0.5 

TTR=73.04 

HTR=22.108 

TMELT=89.8 

HMELT=73.357 

TREF=50.0 

BETA=0. 00045 

VIS=14.3 

CL=CP 

TKL=TK 

BOTTOM PLATE DENOTED BY 1 
DEN1=171.0 
TK1=93.0 
CP1=0.22 

FIN DENOTED BY 2 
DEN2=171 .0 
TK2=93.0 
CP2=0.22 

TOP PLATE DENOTED BY 3 
DEN3=72.5 
TK3=0.09 
CP 3=0. 33 


. Tl(5,31) 
, T3(1 ,31) 


166 



TABLE C-1. (Continued) 


55 GEOMETRY PARAMETERS 

56 W=0. 75/12.0 

57 H=2. 625/12.0 

58 6=5.0/12.0 

59 S=W/(2.0*AM) 

60 51=0.032/12.0 

61 52=0.008/12.0 

62 53=0.25/12.0 

63 F2JM=1 .0/(32.0*5*12.0) 

64' VERTICAL RESISTANCES 

65 RV(1 ,2)=(Sl/(TKl*S2*B))+{2./(HB0T*S2*B)) 

66 DO 10 1=2, MM 

67 10 RV(I,2)=(S1/(2.*TK1*S*B))+(1./(HB0T*S*B)) 

68 RV(1 ,3)=(S/(TK2*S2*B))+(S1/(TK1*S2*B)) 

69 DO 20 I =2, MM 

70 20 RV(I ,3)=(1 ./(2.*TK*B) )+(Sl/(2.*TKl*S*B)) 

71 DO 30 J=4,NN 

72 30 RV(1 ,J)={(2.*S)/(TK2*S2*B)) 

73 DO 40 J=4,NN 

74 DO 40 1=2, MM 

75 40 RV(I,J)=(1./(TK*B)) 

76 RV(1 ,N+3)=(S3/(TK3*S2*B))+(S/(TK2*S2*B)) 

77 DO 50 I =2, MM 

78 50 RV(I ,N+3) = (S3/(2.*TK3*S*B))+(1./(2.*TK*B)) 

79 RV(1 ,N+4)=(S3/(TK3*S2*B))+(2./(HTOP*S2*B)) 

^8 RV(2,N+4)=(S3/(2.*TK3*B*(S+(S2/2.))))+(1./(HT0P*B*(S+(S2/2.)))) 

81 *-D0 60 1 = 3, MM 

82 60 RV(I,N+4)=(S3/(2.*TK3*S*B))+(1./(HT0P*S*B)) 

83 HORIZONTAL RESISTANCES 

84 RH(2,N+3)=((S2+S)/{2. *TK3*S3*B ) ) 

85 ^DO 70 I =3, MM 

86 70 RH(I,N+3)=(S/(TK3*S3*B)) 

87 DO 80 J=3,NN 

88 80 RH(2,J)=(S2/(2.*TK2*S*B))+(1 ./(2.*TK*B)) 

89 DO 90 J=3,NN 

90 ►DO 90 1 = 3, MM 

91 90 RH(I,J)=(1./(TK*B)) 

92 RH(2,2)=((S2+S)/(2.*TK1*S1*B)) 

93 ^DO 100 1 = 3, MM 

94 100 RH(1 ,2) = (S/(TK1*S1*B)) 

95 RH(2,N+3)=RH(2,N+3)+RV(l,N+3) 

96 NODAL CAPACITANCES 

97 C(1 ,2) = ((DENl*Sl*S2*B*CPl)/2.) 

98 DO no 1=2, MM 

99 no C(I ,2)=(DEN1*S1*S*B*CP1) 

100 DO 120 J=3,NN 

101 120 C(1 ,J)=((DEN2*S2*S*B*CP2)/2.) 

102 C(1 ,N+3)=((DEN3*S2*S3*B*CP3)/2.) 

103 C(2,N+3)=(S3*B*(S+(S2/2.))*DEN3*CP3) 

104 ^DO 130 1 = 3, MM 

105 130 C(I ,N+3)=(DEN3*S3*S*B*CP3) 

106 DO 140 J=3,NN 

107 DO 140 1=2, MM 

108 140 C(I,J)=(DEN*(S**2)*B*CP) 

109 ►Q1 = (C(3,4)*(TTR-TREF)) 

no Q2=gi+((DEN*(S**2)*B*HTR)) 
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TABLE C-1. (Continued) 


111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 


— *~Q3=Q2+((C(3,4))*(TMELT-nR)) 

Q4=Q3+ ( ( DEN*( S**2 ) *B*HMELT ) ) 
INITIALIZATION OF PERTINENT QUANTITIES 
TIME=0.0 
QWAX=0.0 
QBW=0.0 
QSW=0.0 
QTW=0.0 
QT0P=0.0 
QFTR=0.0 
QTTR=0.0 
V1=0.0 
V2=0.0 


DMO=0 .0 


DO 150 J=2,NJ 


DO 150 1=1, MM 
150 T1(I,J)=TIN 

— ^^IF{TIN.LT.TTR)QSIN=(C(3,4)*(TIN-TREF)) 

—— IF(TIN.GT. TTR.ANO.TIN.lt. TMELT) QSIN=Q2+((C(3,4))*(T1N-TTR)) 
— ► IF(TIN.GT.TMELT)QSIN=Q4+((C(3,4))*(TIN-TMELT)) 

DO 160 1=1, MM 
160 T1(I,N+4)=TAMB 
DO 170 J=3,NN 
DO 170 1=2, MM 
170 QS(I,J)=QSIN 
DO 180 1=1, MM 


0(I,2)=0.0 
180 Q(I,N+3)=0.0 
DO 190 J=3,NN 
190 0(1 ,J)=0.0 
DO 199 J=2,NI 
DO 199 1=1, MM 


F(I,J)=0.0 
199 QRAT(I,J)=0.0 

READ(5,n) (TMId), 1 = 1 ,ND) 

READ(5,11) (TM2(I), 1=1, ND) 

READ(5,11) (TM3{I), I=1,ND) 

READ{5,11) (TM4(I), 1=1 ,ND) 

READ(5,11) (TIM(I), 1=1 ,ND 

11 F0RMAT(8F10.0) 

WRITE{6,22) TIME,W,H,N,M 

22 F0RMAT(1X,5HTIME=,E15.8,10X,2HW=,E15.8,10X,2HH=,E15.8,10X,2HN=,12, 
15X,2HM=,12) 

WRITE(6,33)Q1,Q2,Q3,Q4 

33 FORMAT(1X,3HQ1=,E15.8,10X,3HQ2=,E15.8,10X,3HQ3=,E15.8,10X,3HQ4=,E1 
5.8) 

WRITE(6,44) 

44 F0RMAT(2X,4HI J ,5X,19HVERTICAL RESIST ANCE,6X,21HH0RIZ0NTAL RESIST 
ANCE,6X,17HN0DAL CAPACITANCE,6X.11HTEMPERATURE,5X,7HQS( I ,J) ) 

DO 200 J=2,NI 
DO 200 1=1, MM 
IF(J.EQ.N+4) GO TO 1 
GO TO 2 
1 QS(1,J)=0.0 
RH{I,J)=1.E08 
C(I,J)=0.0 
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TABLE C-1. (Continued) 


167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 
-179 

180 

181 

182 

183 

184 • 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
^21 


COMPUTATION SECTION COMPUTATION SECTION 


2 IF(I.EQ.1)RH(I,J)=1.E08 
IF(J.EQ.2.0R.J.EQ.N+3)QS(I,J)=0.0 
IF(I.EQ.1)QS(I,J)=0.0 

WRITE(6,55)I,J,RV(I,J).RH(I.J).C(I.J),T1(I,J).QS(I,J) 
onn ^0^f^(1^>I2,lX,l2,5X,El5.8.l0X,El5.8.l0X.E15.8,6X,El6.8.4X.El5.8) 
200 CONTINUE 

DO 889 J=3,NN 
889 T3(1,J)=T1(1,J) 

COMPUTATION SECTION 

3 TIME=TIME+DT 
DO 210 J=3,NJ 

WT(MM,J) = ((T1(M,J)-T1(MM,J))/RH{MM.J))+((T1(MM,J-1)-T1{MM.J))/RV 

210 CONTINUE 

DO 211 J=3,NJ 
^DO 211 I=2.M 

j.>l))+((Tl(I,J-l)-Tl(I,J))/RV(I,J))+((Tl(I,J+l)-Tl(I,J))/RV(I,J+l) 

211 CONTINUE 

212 DO 220 J=3,NJ 
DO 220 1=1 ,MM 

220 QS( I ,J)=QS(I ,J)+(QRAT(I,J)*DT) 

DO 240 1=2, MM 

240 T2(I,N+3)=Tl(I,N+3)+((QRAT(I,N+3)*DT)/C(I,N+3)) 

T2(l,N+3)=T2(2,N+3) 

DO 250 J=3,NN 
DO 250 1=2, MM 

IF{QS(I,J).LT.Q1)T2(I,J)=TREF+(QS(I,J)/C(I.J)) 

IF(QS(I,J) .GE.Ql .AND.QS(I,J).LE.Q2) T2(I,J)=TTR 

IF(QS(I,J) .GT.Q2.AND.QS(I,J).LT.Q3) T2( I ,J)=TTR+( (QS{ I ,J)-Q2)/C( I , 

J j ) 

IF(QS(I ,J) .GE.Q3.AND.QS(I,J) .LE.Q4) T2( I ,J)=TMELT 
WHEN GOING FROM MELT TO FREEZE OR VICE-VERSA CHANGE THE FOLLOWING CARD 

F(I.J) = (QS(I,J)-Q3)/(DEN*(S**2 

) *d*HMELT ) 

IF(QS(I,J)GT.Q4) T2(I,J)=TMELT+({QS(I.J)-Q4)/C(I,J)) 

WHEN GOING FROM MELT TO FREEZE OR VICE-VERSA CHANGE THE FOLLOWING CARD 
IF(QS(I,J).GE.Q4)F(I,J)=1.0 
IF(QS(I,J).LE.Q3) F(I.J)=0.0 
250 CONTINUE 

SPECIFICATION AND/OR DETERMINATION OF FIN TEMPERATURES 

THE FOLLOWING DO LOOP ASSUMES FIN TEMPERATURES FOR ITERATION 
DO 255 J=3,NN 
255 T2(1,J)=T3(1,J) 

DO 260 L=1 ,NDP 

IF(TIME.GE.TIM(L) .AND. TIME. LE.TIM(L+1)) GO TO 4 
260 CONTINUE 

4 FAC=(TIME-TIM(L))/(TIM(L+1)-TIM(L)) 

T2(l ,2)=TM1(L) + ((TM1(L+1)-TM1(L) )*FAC) 

T2(l ,8)=TM2(L) + {(TM2(L+1)-TM2(L))*FAC) 

T2(1 ,15)=TM3(L) + ((TM3(L+1)-TM3(L))*FAC) 

T2(l ,21)=TM4(L) + ( (TM4 L+1)-TM4(L )*FAC) 

DO 270 1=1 ,MM 
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TABLE C-1. (Continued) 


2^3 SspIciRErFiSVEMPERATURES DETERMINED BY STEADY STATE EQUATIONS 
224 256 MC0UNT=MC0UNT+1 

22b DO 280 J=3,MN 

226 T3(l ,J)=T2(1 ,J) 

227 IF(J.EQ.8.0R.J.EQ.15) GO TO 5 

229 l2O;j^4(T2nJ-l)/SV(l,0))*(T2(2,0)/«H(2,J))*(T2(1J*1)/RV(l,^^ 

230 )))/((! ./RV(1 ,J)) + (L/RH(2,jy)+{1 ./RV(1 ,0+1) H 

231 5 CONTINUE 

232 280 CONTINUE 

233 IF(MCOUNT.GT.MFIN) GO TO 8 

234 DO 281 0=3, NN 

235 DIF=T2(1 ,0)-T3(l ,0) 

236 IF(ABS(DIF).6T.EPS) GO TO 256 

237 281 CONTINUE 

238 IF(OOE.EQ.2) GO TO 285 

239 OOE=OOE+1 

240 DO 888 0 = 3, NO 

241 DO 888 1=1 ,MM 

242 888 QS(I,0)=QS(I,0)-(QRAT{I,0)*OT) 

244 Sm2(Ij)^’(T2(2,J)-T2(1.J)1/RH(2.J))4((T2(1,0-1)-T2(1.0))/P.V(1,0) 

247 (MM,0)) + ((T2(MM,0+1)-T2(MM,0))/RV(MM.0+D) 

248 282 CONTINUE 

249 DO 283 0 = 3, NO 

251 ^0RA2(I 0 ) = aT2(I-l,0)-T2{I.0))/RH(I.0))+((T2(I+l.J)-T2(I,0))/RH(I+ 

253 ) 

254 283 CONTINUE 

255 DO 284 0=3, NO 

256 DO 284 1=1 ,MM , „ 

257 284 QRAT(I,0)={QRAT(I,0)+QRA2(I,0))/2.0 

258 GO TO 212 

259 285 MC0UNT=1 

260 DO 286 1 = 1 ,MM , , , 

261 286 QRAT(I,2)=({T2(I,2)-T1(I,2))*C(T.2)/DT) 

262 DO 287 1=1 ,MM , , 

263 287 QS(I,2)=QS(I,2)+(C(I,2)*(T2(I,2)-T1(I,2))) 

264 DO 290 1=2, MM 

265 QBW=(T2(I,2)-T2(I,3))/RV(I.3)+QBW 

266 290 QTW=(T2(I,N+3)-T2(I,N+2))/RV(I,N+3 +QTW 

267 qBT=QBW+((T2(l,2)-T2(l,3))/RV(1.3)) 

268 DO 300 0=3 ,NN 

269 300 qSW=QSW+{{T2(l,0)-T2(2,0))/RH(2.J)) 

270 DO 310 1=1, MM ™ 

271 310 QTOP=((T2(I,N+3)-TAMB)/RV(I.N+4))+QTOP 

272 DO 311 0=3, NN 

273 311 QFTR=QFTR+QRAT(1 ,0) 

274 DO 312 1=1 ,MM 

275 31 2 QTTR=QTTR+QRAT ( 1 ,N0 ) 

276 QWAX=(QBW+QTW+QSW)*DT+(QWAX) 

277 RATIO=QSW/QBW 
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278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 


TABLE C-1. (Continued) 

ERR0R= ( ( QBT - ( QBW+QSW+QTW+QTOP+QFTR+QTT R ) ) *1 00 . 0 ) /QBT 
DO 930 0=3, NN 
DO 930 1=2, MM 
V2=V2+(F(I,J)*(S**2)*B) 

930 CONTINUE 

DAVG=((2,0*V2)/(W*B))*12.0 

QMELT=((V2-V1)*DEN^HMELT)/0T 

V1=V2 

V2=0.0 

DELT=ABS(T2(1 ,2)-TMELT) 

PR=(VIS*CL)/TKL 

RA= ( ( DEN**2 ) *G*CL*BETA*DELT*( DAVG*^3 ) ) /( V I S*TKL*1 728 . 0 ) 
IF(T2(I ,2),LE.TMELT) RA=0.0 

I F ( RA . GT . 1 . EO 5 ) T K= ( TKL*0 . 1 04* ( RA**0 . 305 ) * ( PR**0 .084)) 

I F(RA.GE. 3500.0. AND. RA.LE.1.E05) TK=(TKL*0 .229*( RA**0.252) ) 

I F(RA.GE. 1700.0. AND.RA.lt. 3500.0) TK=(TKL*0 .00238*( RA**0 .816) 

TKR=TK/TKL 

DO 945 1=2, MM 

IF(F(I,3).LT.0.25) GO TO 945 

RV(I ,3) = (1 ./(2.*TK*B)) + (S1/(2.*TK1*S*B)) 

945 CONTINUE 

DO 946 J=4,NN 
DO 946 1=2, MM 

IF(F(I,J).LT.EPS) GO TO 946 
RV(I,J)=(1./(TK*B)) 

946 CONTINUE 

DO 947 1=2, MM 

IF(F(I,NN).LT.0.75) GO TO 947 
RV(I,N+3)=(S3/(2.*TK3*S*B))+(1./(2.*TK*B)) 

947 CONTINUE 

DO 950 J=3,NN 

IF(F(2,J).LT.0.25) GO TO 949 
RH(2,J)=(S2/(2.*TK2*S*B))+(1 ./(2.*TK*B)) 

949 CONTINUE 

— ^DO 950 I = 3, MM 

IF(F(I,J).LT.EPS) GO TO 950 
RH(I,J)=(1 ./(TK*B)) 

950 CONTINUE 

IF(KCOUNT.EQ.KCHK) GO TO 6 

KC0UNT=KC0UNT+1 

GO TO 7 

6 WRITE(6,66) TIME, QWAX, ERROR 

66 F0RMAT(lX,5HTIME=,E15.8,10X,5HQWAX=,E15.8,10X,6HERROR=,E15.8) 
DO 313 J=3,NN 
AJ=J-3 
BJ=J-2 

IF(F(2,J) .GT.F2JM) DFL0=(BJ*S*12 .0) 

DFIN=(AJ*S*12.0) 

IF{F(2,J).LT.EPS) GO TO 314 

313 CONTINUE 

314 DO 315 J=3,NN 
AJ=J-3 

IF(F(MM,J).LT.1.0) DMID=((AJ*S)+(F(f^,J)*S))*12.0 
IF(F(MM,J).LT.1.0) GO TO 316 

315 CONTINUE 

316 WRITE(6,67) DFIN, DFLO, DAVG, DMID 
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TABLE C-1. (Concluded) 


334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 


67 F0RMAT(1X,5HDFIN=,E15.8,10X,5HDFL0=,E15.8.10X,5HDAVG=,E15.8J0X,5H 

DMID=,E15.8) 

DD0T=(DMID-DM0)/DT 


DM0=DMID 

WRITE(6,68) RA, TKR, QMELT , DDOT 

68 F0RMAT(1X,3HRA=,E15.8,10X,4HTKR=,E15.8.10X,6HQMELT=,E15.8,10X,5HDD 

0T=,E15.8) 

WRITE(6,77) QBW, QSW , QTW, QBT , RATIO 

77 F0RMAT(1X,4HQBW=.E15.8.3X,4HQSW=,E15.8,3X,4HQTW=,E15.8,3X,4HQBT=,E 

15.8,3X,6HRATI0=,E15.8) 

88 FORMATUx! , 2 X , 1 HJ , 10 X, 11 HTEMPERATURE, 10 X J5HFRACTI0N MELTED JOX 
,11HENERGY RATE , lOX J3HENERGY STORED) 

DO 320 0=2 ,NJ 

DO 320 1 = 1 , MM ^ ^ 

WRITE(6 ,99)1 ,0 ,T2{ 1 ,0) ,F( 1 ,0) ,QRAT ( 1 ,0) ,QS( 1 ,0) 

99 F0RMAT(1X,I2,1X,I2,8X,E15.8,8X,E15.8,8X,E15.8,8X,E15.8) 


320 CONTINUE 
KC0UNT=1 


7 QBW=0.0 
QSW=0.0 
QTW=0.0 
QT0P=0.0 
QTTR=0.0 
QFTR=0.0 

DO 330 0 = 2, NO 
DO 330 1 = 1 ,MM 
330 T1(I,0)=T2(I,0) 

00E=1 

IF(TIME.LT.TAU) GO TO 3 

8 WRITE(6,111) MCOUNT 
111 F0RMAT(1X,I3) 

STOP 

END 
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TABLE C-2. DISCUSSION OF COMPUTER PROGRAM FOR MELTING 


In the following discussion, references are made to line nuirbers 
corresponding to those designated on the copy of the program. 


LINES 
1 - 4 


5-23 
24 - 54 
55 - 63 
64 - 82 

83 - 95 


96 - 108 
109 - 112 

113 - 144 


145 - 150 


DISCUSSION 


required dimension statements for subscripted 
variables; values should be (MM, NI) for all 
double subscripted variables except T3 for which 
they should be (1, NJ); values should be ND for 
single subscripted variables; NOTE: Q(MM,NI) 

is superfluous and can be omitted with lines 
136 - 140. 

specification of computational parameters 

specification of physical properties 

specification of geometry parameters 

computation of all vertical thermal resistance 
values RV (I.J) 

computation of all horizontal thermal resistance 
values RH(I,J); note that line 95 is a special 
definition which amounts to bypassing node 
(l,N+3) which was done to overcome stability 
criterion required by this small corner node 

computation of all nodal capacitance values 
C(I.J) 

computation of energy stored by a wax node relative 
to TREE for the start and end of phase transition 
and the start and end of fusion, respectively 

initialization of pertinent quantities; the initial 
value of the stored energy depends on the relationship 
of the initial temperature to the reference temperature 
TREE; note that lines 136 - 140 are superfluous 
and can be omitted; some initialization of certain 
parameters is done in the DO loop between lines 
160 and 172 which are set primarily to avoid 
random print-out and are not essential to the 
computation done in the heart of the program 

input data values for measured fin and bottom 
plate temperatures and corresponding time values 
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TABLE C-2. (Continued) 


151 - 172 
173 * 174 

175 

176 

177 - 188 


189 - 191 


print out of initial values for checking purposes 
and print out of certain computed quantities for 
informational purposes 

initially defines T3(1,J) for all fin nodes and 
sets these equal to the initial temperatures 
Tl(l.J) 

beginning of main computation scheme 
time is stepped forward by At 

loops which compute and store the net rate of heat 
transfer to nodes (I.J) based on old temperatures 
T1(I,J); the rate of heat transfer to node (I.J) 
is given by 



+ 


T(I.J-l)-Td.J) ^ 
RV(I.J) 


T(I.J+1)-T(I.J) 

RV(I.J+1) 


this expression must be modified accordingly for 
nodes near a boundary which are not surrounded 
by four neighbors 

computation of total energy stored by node (I.J) 
since the start which is given by 


192 - 193 


^stored 


^ qAt 
time 


computation of new top plate temperatures T2(I,N+3) 
from the expression 


T2(I.J) =T1(I.J) 


194 sets the corner top plate node (l,N+3) temperature 

equal to that of the second node (2,N+3); omission 
of the corner node in the computation scheme was 
done to avoid stability problems due to its small 
size 

195 - 208 computation of new wax node temperatures from 'the 

energy stored by the nodes end their cepeclty and/or 
phase change enthalpy values; when the stored 
energy lies between Q1 and Q2 the new temperature 
is forced to be the transition temperature and 
when It lies between Q3 and Q4 the new temperature; 
is forced to be the fusion temperature; also 
the fraction of the node which has undergone 
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TABLE C-2. (Continued) 


209 - 210 
211 - 212 

213 - 220 

221 - 222 
223 


224 - 237 


238 

239 


phase change is calculated from the relationship 
of the stored energy to Q3 and Q4; note that 
certain designated cards need to be changed when 
running the program for freezing as contrasted 
to melting 

beginning of determination of fin temperatures 

all new fin temperatures are set to T3(1,J) 
which simply represents an assumed value always 
corresponding to the previously computed value 
except at the very beginning at which time it is 
set as the initial temperature 

interpolation scheme which assigns new temperatures 
to the three nodes on the fin and one on the bottom 
plate corresponding to positions where temperature 
measurements were made; the new temperatures are 
linearly interpolated from the input data 

assigns all nodes along the bottom plate the same 
value of new temperature 

beginning of iteration process to determine unspecified 
fin temperatures from steady state equation; fin 
nodes were not treated as transient cases due to 
their extremely small capacitances that would impose 
a severe stability criterion 

iteration process used to determine unspecified 
fin temperatures; in each iteration, T2(1,J) 
is computed from steady state equations and then 
compared with T3(1,J) which corresponds to the 
calculated temperatures during the previous 
iterative step; the iteration is continued until 
the differences between computed fin temperatures 
and their corresponding values for the previous 
iterative step are all acceptably small; should 
the iteration exceed MFIN counts the program is 
directed to stop 

when counter JOE equals 2, the new temperatures at 
all nodes are considered to be the solution at the 
particular time and the program advances to line 259 

increase of counter JOE to 2 


240 - 242 the energy stored at each node (I,J) is reset 

back to its original value; this Is to allow for an 
improved computation of the net heat transfer rate 
to each node to be made and then a recomputation of 
the stored energy and the corresponding new temperatures 

243 - 254 computation of net heat transfer rate to each node 

using new temperatures T2 
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TABLE C-2. (Continued) 


255 - 257 
258 


259 - 278 

279 - 282 

283 

284 

285 - 286 

287 

288 

289 - 290 
291 - 293 

294 

295 - 315 
316 


317 

318 


calculation of net heat transfer rate to each node 
as the average of that based on old temperatures 
T1 and new temperatures 12 

return to line 189 which consists of redirecting 
the computation through that of computing improved 
new temperatures, energy storage values, and ^ 
fractional melted values by using the improved 
(averaged) heat transfer rate (Lines 255 - 257); 
this corrective technique is only employed once 

computation of various heat transfer quantities 
from the new temperatures obtained at time t+At 

computation of volume of melted wax 

computation of average liquid depth from the 
volume melted 

computation of energy required to melt the wax 
which has melted ^ 

resetting of VI and V2 for next time step 

determination of absolute value of temperature 
difference between bottom plate and the fusion 
temperature 

computation of the Prandtl number 

computation of the Rayleigh number 

determination of effective liquid conductivity 
due to convection by using correlations of 0‘Toole 
and Silveston 

computation of ratio of effective liquid conductivity 
to actual value 

recomputation of thermal resistances in the liquid 
by using the effective thermal conductivity rather 
than the actual value 

counter check which controls printing out of 
desired results as well as computation of 
interfacial position 

counter advance 

by-pass of printing results except when line 316 
is executed 
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TABLE C-2. (Concluded) 


319 

321 

333 
336 • 

338 . 
352 - 
359 - 

362 

363 

364 - 

366 

367 


- 320 

- 332 

- 335 
■ 337 

- 351 
' 358 

361 


365 


write statement for printing results 

computation of interface position for columns 
next to fin and next to centerline 

write statement for printing results 

calculation of interfacial velocity and renaming 
interfacial position to provide for determining 
its change at the next time step 

write statements and corresponding formats for 
printing results 

reinitialization of pertinent quantities for next 
time step. 

setting new temperatures for current time step to 
be old temperatures for the next time step 

reinitialization of counter 

comparison of time to upper limit value which 
when exceeded results in stopping the program 

printing out of value of counter used in fin 
temperature iteration 

STOP 

END 
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Figure C-1. Skeleton flow chart. 
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TABLE C-3. FORTRAN COMPUTER PROGRAM FOR SOLIDIFICATION 


C T£bT ^30-xb TtSf INCH CELL 

DI'^ENSIUN MV(5«31)i RH(5i3l)« C(5«31)/ Qb(5>31)« TllS 
1# 31 ) 

DInENSlON T2(bi31)i F|5>31>« Q^AnS«31>/ T3(1#31 >«QRa2 
1 ( 31 ) 

DI'iENSlON Trll25l« TM2(25I* T*13(25)# TMA125)/ TIH125) 

3I1-NSIQN vR(Si31)< h«»15/31I 

C COmPuTaTI JnAl PARAnETE*?S 
n-27 
A 

A'1»'1 

1 

ns»\+2 
N J»N-*-3 
Nl -NAA 
NO-23 
nOP-22 
TAu-1 *5 
KCOjM- 1 
lCOJNT-1 
iF I N-50C 
EPS-1 t£-06 
kCm<- 1 UuC 
jOE- 1 

DT-0.0C005 

TAlj2-i lOi-DT )♦(DT/3. ) 

AF-*tC« 

UN-1 «Q«EPS 

C PHYSICAw PROPERTIES 
TAna-ac. 

TIN-100. 

HTOP-I.E-08 
MBOT -5 • 

Q-(3£>17A.^3600i.^3600*) 

C WAX 

DEn-a 7 i 2 
T<-0.087 
CP-0.5 
TTR-73.0A 
hTR-22. 108 
T'AE^T-SS *8 
h'1EwT»73.357 
TREP-50. 

BETA-0.00CA5 
Vis-IA.3 
CL-CP 
T<l - T < 

C rtOTTDI Plate DENOTED BY I 

OEnI-171 . 

T<1-93. 

CPl-0.22 
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TABLE C-3. (Continued) 


C ns DiNJTEO BY 2 
DEN2-171 . 

T<2-93. 

CP2-0.22 

C TOP PLATE DEnCTEO BY 3 
DEn3»72«5 
TK3-0»09 
CP3«0*33 

c geometry parameters 
w-0t75/12*0 
H»?»625/12*0 
B»5t0/12»0 
S»w/ ( 2* ♦AM ) 

Sl-0>032/12«0 

S2>0*008/12*0 

S3«0*25/12*0 

F2JM-1 .0/( 32*0»S*12»0» 

c vertical resistances 

RV(1>2)-(S1/(TK1»S2*B) ) ♦ ( 2 • / ( HaOT*S2*B ) ) 

DO 10 I»2>MK 

10 RV( Ii2l-(Sl/(2t*T<l»S»Bl )♦< li/(HB0T*8*B) > 
RV(li3)«IS/lTK2AS2AB) ) ♦ ( Si/ I TK1*S2*B ) I 
DO 20 I»2/MM 

20 RVlIi3)-(l./(2*»T<*B) ) ■*• ( SI/ ( 2 • »Tk1 »S*B ) ) 

DO 30 j«A/NN 

30 RV l l> J ) • ( I 2. »S 1 / ( TK2»S2*8 ) » 

00 YO J»4<NN 
DO 40 1»2<MM 
40 RV( I#J)«( 1*/(TK»B) ) 

RV( WN^3)-lS3/tTK34S24B) I 4 ( S/ { T<2482*B ) ) 

DO 50 I-2/MM 

RV(r*N44) ■ (S3/( 2.4TK34S*S) > ♦ < 1 • / ( HT0P»84B >J 
50 RVt I#N43)-(S3/(2.»TK34S*B» I ♦ ( 1 . / ( 2 1 4TK*B > ) 

RV(liN4 4)>(S3/(T<34S24B) ) ♦ ( 2* / ( HT0P»62«B ) ) 

RV(2iN + 4)-(S3/(2»*Ti<3*B*(S+(S2/2* ) ) 1 ) ♦ ( 1 • / < HTDP*B4 ( 84 
1 (S2/2* ) ) ) ) 

C horizontal resistances 
DO 70 1-2/MM 

RH(I/2) - ( S/ ( TK14S14B ) I 
70 RH( I/N43)-(S/(T<34S3*Bn 
DO 90 J-3/NN 
DO BO I-2/MM 

BO RH ( li j ) ■ ( 1 •/ I TK4B n 
DO 90 J-3/NN 

90 RH(2/J)« ( 52/ ( 2»»Ti< 2*S4B ) ) 4 ( 1 . / ( 2 • *T<*B ) ) 
RH(2/2l«MS24b)/:2t4TKl*Sl4B)» 

RH< 2/ N43 )■( I S24S ) /( 2«4TK3483»& I J 
RH(2/N43)-RM(2/N43)4RV{1/N43) 

c nodal capacitances 

Cl l/2)a( I DEN14S14S24B4CP1 l/2i ) 

DO 110 1-2/MM 
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TABLE C-3. (Continued) 


llO C ( I « ^ B ( OEn 1 «S 1 1 I 

DO 120 j»J>NN 

120 C(l«J)a( ( 0 En2<»S2«S*B«CP2 )/2* ) 

C « 1* Nt-3 )■ ( ( DEn3*S2«S3»B«CP3)/2* ) 

DO 130 lB2>hN 

130 C ( I » N*3 I ■ ( DEn3»S3*S*B*CP3 ) 

C(2#N + 3)»(S3¥tJ»(S^(S2/2»l )*DEn3*CP3) 

DO 140 j«3iNN 

00 140 iB2«nn 

140 C I I * J ) ■ < DEN* ( S*»2 ) 4B4CP J 

01 ■ (C(2<3}*(TTR-TRED ) 

Q2«0l4( (0EN*IS**2)*B*HTR) I 

03 ■ 02 ♦ ( (C(2<3) J^ITHELT-TTR) ) 

Q4»33+( (0EN*(S*¥2) *B*hHELT I ) 

c Initialization of pertinent quantities 
T I ^E*0 • 0 
QWAXbOcO 
Q3w*0«0 
05 «< *0 • 0 
QTw»0 » 0 
OTOP-0.0 
UF TR«0 • 0 
QTTR«0»0 
Vl*OtO 
V2-0.0 
DflO-O.O 
DO 150 J-2/Nj 
DO 150 

150 T1 1 1/ J J-TIN 

1F(TIn»LT»TTR)QSIN« ( C( 2a3)«( TIn-TREF) ) 

IF( TiNtQT •TTRiAND*TIN*LT.TMELT ) QSIN*Q24 ( C ( 2/ 3 ) * ( TI N 
l-TTR ) ) 

IF ( TiNtQT • TNELT ) QS1N-Q4+ i ( C ( 2> 3 ) J* I TIN-TMELT ) » 

DO 160 IbI/MM 
160 Tl I I jN+A I-TaMB 
DO 170 j» 3*NN 
DO 170 I«2iMh 
170 uS ( I« j ) aQSlN 
03 199 J-2/M 
DO 199 I»1/PM 
F( I>J|aO*0 
199 ORAT < 1/ j )«0»0 

READI5/11) ITm1(I)> 1b1«nO) 

REAOtS/11) (TM21I), I-1<N0» 

REA0(5>11) (TM3(D* I-1jND) 

REaD(5/11I (TN4(I), IbIjNO) 

REA0(5^11) (Tl^dl# Ib1«NO) 

11 F3RnaT(8F10i0I 
DO 910 KTC"1* nD 
910 Til", < <TC » -TIM ( KTC I /360C* 
wRlTE(6i22) T I nE a w / h > N/ M 
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o r> 


TABLE C-3. (Continued) 


22 FOM'iAT ( IXiSm I1t"/E15»B*10X/?K**"4E15«8#10X#2HH«#E15*8 

12/ 15x>2mM»/ 12) 

w I T E I 6 / 3 3 ) C 1 / G 2 / 3 3 / 0 * „ 

33 F0«1AT(lX/3MQl«/E15»Bi 1 OX# 3H02»# E lb • 8/ lOX/ 3HQ3"/ E 15 • 8 
I# lOX# 3h04«#E 15»B ) 
w'?IT£(6#4‘») 

Hh F0^<1AT t EX# ^hl j# 5x# ISHVERTICAw RES 1ST ANEE# 6X 

l#2lr(H0RIZ0NTAL RES1STANCE#6X# 17MNOOAL CAfAC I T ASCE# 6X 
2#llHT£nPERATbR£#5X#7H0S(I#J) ) 

CD 200 j»2#M 
DO 200 I■l#^'.^l 
IF ( j.EU»N + 4 ) GO TO 1 
GO TO 2 

1 QS(I#J)>C«0 
RH( I#JI"1»E08 
C( 1# J)«0*0 

2 IF( I.EGtl ) RH( I# j)"l*EOB 
IF(J»EQ.2«0R»J'EQ»N+3) QS( I#J)-0»0 
IF(I.Euil) 3SII#J)>0*0 

wRITEt6#55)I#j#Rv(I#J)#WH(I#J»#C(I#J)#Tl(I#J»#QS(I#J) 
55 F0R'iaT(1x#I2#1X#I2#5X#E15»B#1CX#E15»8#10X#E15»8a6X 
1#E15»8#AX#E15*8) 

200 continue 

DO AOO J«3#NJ 

DO aOO I«2#^*N 

VR( I#J)»RV( I#j) 

HOP hr ( I # J ) ■FiH I I# J ) 

00 889 j"3#NN 
889 T3 ( I# J I -Tl ( I# J ) 

conputation section computation section computation 
SECTION 

3 TIME-TIME+DT 
DO 210 J>3#NJ 

ORAT( I# J )•( ( Tl ( 2# J )-Tl ( l#U) l/RH(2# J) )♦( ( Tl( I# J-l )-Tl ( I 

l/Dlj/RydiJIl-HlTlIl/J + D-TKl/J) )/RV( I# J4>1 ) ) 
QRAT(Mn#j)-( (T1(M#J)«T11MM#J) l/RH(rin#J) )♦( (Tl(MM#J-l) 
1»T1(MN#J) )/RV(MM#J) )♦( (TUMM/J^li-TUrtM/J) )/RV(MM#J 
24-1 ) ) 

210 continue 

IF(MtEQ*l IGO TO 212 
DO 211 J-3#NJ 
00 211 I»2#M 

QRaT(I#j)-((TI(I- 1#J)-T1(1#J) I/RH(I#J) ) + MT1(I + 1#JI-TI 
llI#J))/RM(Itl#j>)»((Tl(I#J«ll*Tl(I#J»)/RV(I#J))*(<Tl(I 
2#J^l)-Tl(I#JI I/RV(I#J + IM 

211 continue 

212 00 220 J-3/NJ 
DO 220 l«l#Mn 

220 QStl#j)»QSll#JI+l3RAT(I#J)*DTl 
DO 2 a 0 I»2#Mri 

2*»0 T2tl#N + 3)«TltI#N43)<>(( QRAT ( I # N t3 I *0T I /C « 1# N + 3J ) 
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TABLE C-.‘L (Continued) 


T2( N+3 ) «T2 ( S+3 ) 

D3 250 J*3iNN 
00 250 I«2/rM 

IPt3S(Iij).LT«Jl) T2( I* JI»TREF^(0S( I# J)/C( J) ) 
ir(aS(I/j)>CE>31«*SOOS(IiJ)«LE*Q2l T2{I#J)>TTR 
IFOSIIiJ)<GT<02*ASD«QS(l#J)*LT*Q3) T2(I>J)>TTRt((QS(I 
l>J)-02)/C( 1/ J) ) 

IF ( 35 ( 1« J ) •GE'ua* *^0«QS( J) «LE*Q4 ) T2 ( 1/ J ) >Ti£LT 
C WHEN GOING FROh ^ElT TO FREEZE OR ViCE-VERSA CHASQE THE 
C FCuLOwiNQ CA 

l'=‘l3S(]/j)iGT*33.ANC*0S(IiJ)»LT»0A» F( 1/ J)-(04-0S< I 
1 /J))/(DEn»(S»* 2)*3»HMELT1 
IF( 3S I 1/ j ) tGE'OA J F(I«j)aO«0 

IF(3S(lij)iGTtQ4) T2( I# J (•TMELT+C (QSC I # J ) -04 I /C 1 1 # J I ) 
WHEN going from mElT TO FREEZE OR VICE-VERSA CHANGE THE 
following CA 

IF I 35 ( j ) tLEtOS ) F(Iij)-ltO 
250 CONTINUE 

SPECIFICATION ANO/JH OETERMINaTICN OF FiN TEMPERATURES 
The following DO LOOP ASSUMES FIN TEMPERATURES 
FOR ITERATION 
DO 255 U-3^ NN 

255 T2(l/ J)-T3( li Ji 
00 260 l-I/NDP 

IF( TIME.GE*TIM(L I .ANDiTImE»LE«TIM(L+1 > ) 30 TO 4 

260 continue 

4 FAC-(TlME-TlM(L))/(TlM(L4ll-TlM(L)) 
T2(l/2l-TMl(Ll-M(TMl(L4ll-TMl(Ln4FAC» 
T2(l/8)-TM2(L)f((TM2(L4l I-TM2 (l ) I4FAC) 

T2«1/15)-TM3(L1+( (TM3(L41 >-TM3(L) J^FAC) 
T2ll/21)-TM4(L)4((TM4(L4l I-TM4(L) )*FAC ) 

00 270 I-liMM 

270 T2l I>2)-T2( 1#2) 

C unspecified fin temperatures OETERmINEO 8Y SfEAOY STATE 
c equations 

256 MCOUNT-MCOunT+I 
DO 280 U-3^NN 
T3( J |-T 2 ( i>u I 

IF t J»E0*8«0R* J«E3t 15 ) OO TO 5 
IS- ( j.EQ*21 ) GO TO 5 

T211/U)-( IT2(1/j-1)/Rv( 1>U) ) ♦ ( T2 ( 2/ U > /RH I 2# J ) )♦( T2( 1/U 
l+l)/RVIi>U4llll/l(lt/RV(l#U) )♦( l«/RH(2*U> )♦( l»/RV{ I# J 
241 )) ) 

5 continue 

280 continue 

IF (MCOunT«QT*MFIN ) QO TO i 
DO 281 J>3iNN 
DIF-T2I 1/U)-T3( liU) 

IF ( ABS ( DIF ) tQT*EPS I 00 TO 256 

281 continue 

IFiu0EiEQi2) QO TO 285 
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TABLE C-3. (Continued) 


JDE-JOE+I 
D3 868 j>3«Nj 
DO 888 I«l«hn 

0 88 QS ( J ) aQS ( I« J ) * t 3RAT ( J l«OT I 
DO 282 JaO^Sj 

(T2(2> J)-T2JIaJ)»/R 8C2#J) !♦«( 121 1# J-1 )-T2< 1 
l/J)l/RV(liJ)) + (lT2(l#J*l »-T2( 1# J) »/RV( Ij J+1 1 ) 

QRA 21 j I ■ ( ( T2( Ml J )»T2(riH/ J J )♦( CT2lM.ii J-l ) 

1- T?(K'i,j) |/Rv(MMiJ) )♦( ( T2( MMi J*l l-T2(Mi# J » )/RV ( MM# J 

2- f 1 ) ) 

282 CO-'iTIwuE 
irti‘EQ.1 IGO Tj 283 
DO 283 ja3iNJ 

DO 283 I«2iM 

ORA2(IiJ)a( (T2il-liJ)»T2(liJ)» /R mC I iJ))+l(T2(I^liJ)aT2 
1 I I* j ) )/R«( l + li J ) ) f ( (T2Ui J-1 1-T2C Ii J) >/RV( Ii J) )♦( (T2(.I 
2i J+l ) -T2 ( Ii j » ) /RV 1 li J*l ) I 

283 CONTINUE 

DO 28^ j«3iNJ 
DO 28<* I-lihM 

28^ OR at I Iij)-IORAT| I i j)tQRA2( IiJI )/2*0 
GO TO 212 
285 MCOdNTal 

DO aOI j»3iNj 

DO 401 I»2ihM 

IFCFC Ii J) •GT«EP 5*ANDiF< I iJI«LT«UNI QO TO 402 
IF(F( Iij) .GT«UN) 30 to 403 

GO TO »01 ^ 

‘lOa RV( Ii J)"VR( IiJ) 

Ri(IiJI«hS(IiJ> 

GO TO API 

«»02 RVCliJI"VR(IiJI/Ag 

RHC Ii JI»HRC I iJ»/AF 

RVC Ii J* 1 )«VRC li J»1 )/AF 

RHII-»liJI«riRCl»liJ)/AF 

»01 CONTINUE 

DO 286 I-liMM 

206 QR.ATt Ii2}-( i T2( I# 2 I -T1 ( Ii2i t«C( Ii2>/0T) 

DO 287 laliMfl 

287 0S( Ii2)-QSCIi2)4-IC( 1i2)«(T2(Ii2)-T1CZa2) > ) 

DO 290 I»2iMM 

QBw-(T2( Ii2)-T2l Ii3) ) /RV ( 1 1 3 1 ♦3BW 
290 0T.<«(T2( I>N*3)-T2( IiN*2» I /RV ( 1 1 N+3 ) ♦QTw 
oaT ■UBW+ ((T2(li2)-T2(li3l )/RVI li3> ) 

DO 300 Ja.IiNri 

300 QS^^-QSmx (T2I 1ij)-T2(2iJI )/RH<2iJ) ) 

CO 310 I«liMM 

310 qT0»-(«T2( IiNi-3)-TAMBI/RV< IiN^*) I +070? 

DO 311 jaJiNN 

311 QFTR-QFTR+QRaTC li J) 

DO 312 l-liTM 
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TABLE C-3. (Continued) 


J12 orTR-UTTR + URAT ( I> NJ ) 

0»<AX«(iiBw^OTw>QSH)*DT»(awAX) 

RATIO-QSm/OBk 

ERR0R«( (C3BT-(QSw + iSH^QTw + Q.TOP» 3 FTR^(jTTR) )»100»0)/3BT 
DO 930 j«3/Nn 
DO 930 1»2>MM 
V2»V2^(F(I/J)»1S»*8)*B) 

930 continue 

DavS«( (2*0»V2I/(W*B) I*12»0 

0NELT«( ( V2-V1 1 aOEN¥HMELT)/DT 

OFREZ-UMELT 

V1»V2 

V2«0 • 0 

IF( <C0uNT»EC*KCM< ) GO TO 6 

KCOjNT«KCCUNT-*>1 

GO TO 7 

6 w«ITE(6/66) TI^lEi OwAX# ERROR 

66 FOR NAT I 1X/5HTInE-/E 15.6< ^OX/6h3WAX•/E^5•^^/lOX#6HERROR■ 
l/E15•8 ) 

DO 313 J-3>NN 
A J» J“3 
BJ-J-2 

IF(F(2>j),qt.F2jN) 0FL0-( BJ*S»12«0) 

OFIn-( AJ«S*12*0) 

lF(F(2i J) .LTtEPS) GO TO 31A 

313 CONTINUE 

314 DO 315 J-3>NN 
A JB j«3 

IF(F(MNij)*LT*l«0) DM10 b((AJ*S)+{F(HM^j)* 8>)»12»0 
IF(F(Nr.i J) •LT*lt0) GO TO 316 

315 CONTINUE 

316 wRITE(6i67) DFlN> OFLOi DAVQi OHIO 

67 FORNAT( lXi5HDFlNB«E15i8> lOXi S hOFLO*/ E 15 • 8/ lOX* SHOAVQ b 
1/E15t8/ 10Xi5HDNl0>«E15>8 I 

DD0T-(DMID“DM0)/DT 

DN0«0MID 

WRITEI6/69) QFREZ# DCOT 

69 FORNAT( lXi6M0FRE2«*El5*Bil0X#5HDD0T»<El5»8) 

WRITE(6i77) QBw> 3 Sh« QTw* OBTj RATIO 
n FOR NAT ( lXi4H0Bw-<£l5»8>3x/4MQS«-AE16*8i3Xi4M0Tt<B/E15»8 
1/3 x#4HQBTb>E15»8* 3X*6HRATI0b#E15*8) 

WRITE ( 6i88 ) 

86 FDRNAT(2Xi lhI#2Xi iHJi lOX^ llHTENFERATuRE# IOX a 15HFRACTI0 
IN NELTEOi lOXi I IMENERGY RATEj lOXi 13 hENERQY STORED) 

DO 320 jb 2/ N J 
DO 320 lal/hri 

WRITE (6#99) 1/ J# T?( I# j ) # FI 1/ J)# 3RAT( J)/Q8( I> J) 

99 FDR NAT I I X/ 1 2« IXi 1 2« 8X« E 15 *8/ 8Xi EIS* 8« 8x« E15 • 8« 8X 
1«E15>8) 

320 continue 
KCOUNT- l 
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TABLE C-3. (Concluded) 


QSW-0*G 
QTm«0 • 0 
qTC®»0 • c 
QT7^»0«0 
QFTR«0«0 
DO 330 jaSiNj 
DO 330 

330 T1 ( li J ) «T2 ( 1/ J 1 
jOE-i 

IF ( TI^L »LT *Tau I GO TO 3 
8 wRITEI6<ni) MCOunT 
111 FOW'IAK IXi 131 
STOP 
END 
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NOMENCLATURE 


A 

ADI 

Bi 

C(l ,J) 
C 

P 

GH 

GV 

9 

H 

h 

I 

J 

k 

L 

M 

N 

Nu 

PCM 

q 

1 1 
q 

Ra 


Area 


Alternating direction implicit numerical method 

r.. .1 hX 

B 1 ot modul us = \ — 
k 

Thermal capacitance of node (l,J) 

Heat capacity 

Horizontal Thermal Conductance 
Vertical Thermal Conductance 
Acceleration of gravity 
Height of PCM in Cell (Figure 3l) 

Convective heat transfer coefficient 
Latent heat of fusion 

Designation of nodal location (See Figure 31) 
Designation of nodal location (See Figure 31) 
Therma 1 conduct i vi ty 

Distance between bottom surface of cell and liqui 
interface at center of cell 

Number of horizontal nodes in PCM (Figure 31) 

Number of vertical nodes in PCM (Figure 31) 

1 hx 

Nusselt number = i — 

Phase change material 

Heat transfer rate 


Heat transfer rate per unit area 

„ . . , . gBx^AT 

Rayleigh number = ^ 

' va 


d-sol id 
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NOMENCLATURE (Concluded) 


RH Horizontal thermal resistance 

RV Vertical thermal resistance 

S Nodal spacing PCM 

Thickness of bottom plate 
$2 Thickness of fin 

S3 Thickness of top plate 

T Temperature 

t time 

W Width of PCM in Cell 

X Significant length 


GREEK SYMBOLS 

a Thermal diffusivity » 

A Denotes a finite increment 

u Dynamic viscosity 

V Kinematic viscosity 

p Density 


SUBSCRIPTS 

f fusion value 

g Glass 

5 Surface 


SUPERSCRIPT 

Denotes calculated value at time t + At 
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TABLE C-4. CINDA COMPUTER PROGRAM FOR FREEZING 


rtCD ^THfcPMAL LPCS 

BCD 9 THFhMAL CmPACUOK y,AX hODF 3/9 FIN SPACING STUDY 
flCD 9 TFST 4(;?3u- 7) POF FLUID Ir^LLT NONAUtCAML PLASTIC oOX 
BCD 9 FIRST CLLL- 3 UU/Hj/HR FLUID FLOWH ATL-K-PLfcX IGLAS= • 09 
i»CD 9 K->rAX=.U8 »9X13 NODLS IIM wAX-IUpUf KIN IEMPS AMD BASE 
BCD 9 TtwPS FROH test DaTA 
END 

BCD SNOL.F data 
31»100. » .00^ 

GtN 

GEN 

GEN L0»E » 1 » lliD. »-.2 » ^ * i • » 1 . 

GEN »1U(). »-l . » 1 . » i . » 1 . 

GEN t la ^ l a f la 

-IrlOOarOa DbBOrrO 

~t*»10{la»0a SFIl'J 

-?ri(JU.»0a IFI 

-14»]U0a»Ua iFIU 

”ls»100a»0a 'tFIi-i 

-lfl»10U.»Oa tFlD 

“19»lU0.»0a 'fcFiM TENiP 

Gt’N » S2 »l»lUOa»Ua»la»la»la 

END 

BCD ^CONDUCTOR DATA 

GEN ii»2»l»l»0»2*l»3a30*la»la»l« 

GEN 2lil»4»l»l»0»201»l»aii5»la»l«*l« 

GEN 4 f 24 » 1 »2 »l»4*l»la S8 »la*la*l* 

GEN 20Sr4Br 1 »201 » l»20b» i » a 031 3 » 1 a » 1 . » 1 . 

301»?»20l^ aU71 

GEN 307 » 2 r 1 » 204 »l»3»l»a071*la*l«»la 
GEE' 312»2»1» 20rt »l»b»l» a*t7l»la»l»rla 
GEN 317»2»l»212»l»7»l»a m 7 I r 1 a » 1 . » i a 
GEN 322 »2»i»2lh»lf9»i»aU71»la*l»»la 
GEN 327»2» l»220r 1 » 1 1 * 1 » a l)7 1 » 1 a » 1 a » 1 a 
GEN 332»?»1»224»1 » lo» 1» a l) / 1 » 1 a »1 i » 1 a 
GEN 337»2» 1 »228» 1 » 13 » 1 » a 07 1 » 1 a » 1 a » 1 . 

GEN 342r2»l»232» 1 » 1 7 » 1 r a 0 / 1 » 1 a » 1 a » I a 
GEN 3b7»2»l »244»1 »23» 1» a 071 » 1 a » 1 a » I . 

GEN 347»2» 1 »236» 1 » 19 » I » a 07 1 » 1 a » 1 a i 1 . 

GEM 362»2 » 1 » 248r 1 »23» lra0 71»l.»l.»la 
GEN 3b2»2» l»240r 1 »21 » 1 » a 071 » 1 a » i a » I a 
J^b7r?b2»27f a07l 
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TABLE C-4, (Continued) 


3EN 

J)02p2p 1 

»;-'Ui 

* 1 p2U2 p 1 p 

.1)31 

f 1 . » 1 . 

f 1. 


<5U6r?D3 

p20‘^ 

p .U31 




i3tN 

309f 1 

p2i]b 

rl»2Ut:» Ip 

.031 

f 1. f 1 . 

f 1 . 


-d14»3» 1 

p 2.1)9 

»lr21U» If 

.031 

f 1 . f 1 . 

f 1 . 

vjEN 

319»3» 1 

p2lb 

f 1 f 214 f 1 f 

.031 

» 1 . f 1 . 

f 1 . 

Gt.N 

6ii4 » .3 r 1 

p2\1 

f 1 f 2 i 6 • 1 f 

. 031 

f l.f 1. 

f 1 . 

3E(n' 

J>29p:^,1 

p221 

f 1 f 222f 1 f 

.031 

f 1 . f 1 . 

f 1. 

GLIn) 

334 f 3 f 1 

p22b 

p i f 22b p 1 f 

.031 

. i. f 1. 

f 1 . 

Gt N 

339 » 3 » 1 

p2?9 

f 1 f 231) f 1 f 

.031 

f 1. r 1 . 

f 1 . 

OtN 

344 » 3» 1 

p?b2> 

f 1 f 234 p 1 p 

.031 

f l.f 1. 

ft. 

(it N 

349»3» 1 

p2b7 

f 1 f 2 3 b f 1 f 

. «j.3 1 

f 1 . f 1 . 

ft. 

Gh f\i 

334 » 3 » 1 

p ^>41 

p 1 » 2.42 p 1 p 

.031 

f 1 . f 1 . 

f 1. 

GhN 

3b9»3» 1 

p f^4b 

f 1 f 2 4 f 1 f 1 f 

. 031 

f 1 . f 1 . 

f 1. 

GF.N 

364 » 3 » 1 

p299 

f 1 ' 2 b t.! f 1 f 

.031 

f l.f 1. 

f 1. 

(ihN 

^;b3»?» 1 

p2^9 

f 1 » 2 8 f 2 f • 

0294 

f i. f 1 . 

ft. 


<rbb»?bl 

p51p 

. u294 





<r.b6 p 2b2 

P 29 P 

. 0294 




ht 

3 ft p 2 p 1 p 

22 p 1 

r ,5 U r i r * n r :> r 1 • f 

1 . f 1 . 



r <51! f J5i » • 0 t) 

I'J tip 1 t 2b p 1 p tip . (I r I , » 1 . » i , 

h’ui ) 


1 » I) . 

2 pU. 

i • vJ , 

4. ■» 1) , 

^ I 0 # 

f > » 1 1 , 

7 »U » 

hpO . 

9p0. 

llrO. 

J ? » 0 . 

1 3 » ft , 

1 b » n , 
16 r (I > 


190 



17»0. 
1 6 » 0 » 
(j. 
^0 f 0 . 

2i»n. 
22»0. 
2i r 0 . 
24»0. 
25»0. 
26r0. 
27r0. 
28f u. 
29f n. 
Jo»n. 
31^0. 
32 » 0 . 
33 » 0 . 
34»n. 
35 » 0 . 
3b»0. 
37rfi, 

38 r 0 • 

39 f 0 . 
4U»0. 
41 »n. 
42»P. 
43 r 0 . 
44»(J, 

4 5 » r < * 
46 » 0 , 
47»n. 
48»[t. 
49»n. 
bO » 0 « 
51 »f, 
52 » 0 . 
53 » n • 
54 » 0 . 
55»0. 
56 » 0 • 
57» 0 » 
58 » 0 . 
59 » 0 . 
60 » 0 . 
61 » 0 • 


TABLE C-4. (Continued) 


62 » 0 . 
63» 0 . 
64»0, 
65» 0 • 
66 » 0 • 
67 » 0 . 
68 » 0 . 
69 » n » 
7 0 M) . 

71 #0. 

72 Ml. 
7 3 » 0 . 
74f 0. 
75 » 0 . 

7 6 » n « 
77*0. 
78 » 0 . 
79*0. 

8 0 * 0 . 

81 *u. 

82 * () . 
83*0. 
84 * 0 . 
85*0. 
86 * 0 « 
87*0. 
88 * n . 
89*0. 
90*0. 
91*0. 
92*0. 
93*0. 
94*0. 
95*0. 
96 * 0 . 
97*0. 
98*0. 
99*0. 

100*0. 

101*0. 

1 02 * 0 . 
103*0. 
104*0. 
105*0. 
106*0. 


107*0. 
108*0, 
109*0. 
110 * 0 . 
1 11 * 0 . 
112 * 0 . 
113*0. 
114*0. 
115*0. 
116*0. 
117*0. 
118*0. 
119*0. 
120*0. 
121*0. 
122 * 0 . 
123*0. 

124 * 0 . 

1 25 * I) . 
126 * 0 . 
127*0. 
128 * 0 . 
129*0. 
130*0. 
131*0. 

1 32 * 0 , 

1 33 * 0 , 
134*0. 
135*0. 
136*0. 
137*0. 
138*0. 
139*0. 
140*0. 
141*0. 
142*0. 
143*0. 
144*0, 
145*0. 
146*0. 
147*0. 
148*0. 
149*0. 
150*0. 
151*0, 
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Ib2r0. 

153»0. 

154 Ml. 

155 » 0 • 

1 56 » 0 » 

1 5T ^ n . 
I58»n. 
159»P. 
ibn»o. 
16l»0. 
I62»n. 
163»0, 
164»n. 
it>5 » r» » 

166 r 0 . 

167 MU 

168 MU 

169 MU 
I7n»n. 
I7l»n. 
172»0. 
173MU 
I74r0. 
175»0. 

176 MU 

177 MU 

178 MU 
179»0. 
18»M). 

181 MU 

182 MU 
183MU 
184 HU 
1B5»CU 
186»n. 
1 8 7 » n , 
188»0, 
189»P. 
19n»0. 

191 MU 

1 92 » 0 . 

193 MU 

194 MU 

195 MU 
196»P. 
197»n. 


TABLE C-4. (Continued) 

198 »CU 
199»0. 

200 » 0 . 

201 » 0 . 

202 » 0 . 

203»n. 

204»0. 

205»0. 

206»n. 

2U7 » 0 • 

208»n. 

209»0. 

210»0. 

2ll»0. 

212 * 0 . 

213*0. 

214*0. 

215*0. 

216*0. 

217*0. 

218 * 0 . 

219*0. 

220 * 0 . 

221 * 0 . 

222 * 0 . 

223*0. 

224*0. 

225*0. 

226*0. 

227*0. 

228*0. 

229*0. 

230*0. 

231*0. 

232*0. 

233*0. 

234 *(U 
235*0. 

23b * 0 . 

237 * 0 . 

238 * 0 . 

239*0. 

240*0. 

241*0. 

242*0. 

243*0. 


244*0. 
245 * 0 . 
246*0. 
247*0. 
248*0. 
2‘*9*ru 
250*0. 

251 * 0 . 

252 * 0 . 

253 * 0 . 

254 * n . 
255*0. 
256*0. 
257*0. 
258 *'0. 
2-59 » 0 • 
260 * 0 . 
261 * 0 . 
262 * 0 . 
263 » f) . 
264*0. 

265 * 0 . 

266 * 0 . 
267*0. 
268 * 0 . 
269 *»U 
270*0. 
271*0. 
272*0. 
273*0. 
274*0. 
275*0. 
276 * 0 . 
277*0. 
278*0. 
279*0. 
280 * 0 . 
281 * 0 . 
282 * 0 . 
2183 * 0 . 
284*0. 
285*0. 
286 * 0 . 
287*0. 
288*0. 
289*0. 


192 



TABLE C-4. (Continued) 


<:90»n. 

292 M). 

293»0. 

294 » 0 . 

295»n, 

296 t 0 , 

297r0. 

298 » 0 • 

299»n. 

i00»0. 

301 » f) , 

302»n. 

303 » 0 • 

304»0. 

306f 0, 

306»n, 

307»fi. 

3U8»0. 

309f 0. 

31 0 » 0 • 

31 1 » 0 . 
312r0. 
313» .617 
314» .617 
315» .617 
316f .617 
317» .617 
318» .617 
319» .617 
32Df .617 
321» .617 
322» .617 
323» .617 
324» .617 
32b» .617 
326» .617 
327» ,617 


328» .617 
329» .617 
330» ,617 
331 r .617 
33? f .617 
333» .617 
334» ,61? 
335» ,617 
336» .617 
337 ^,61/ 
338» .617 
339r .617 
340f .61/ 
341r .617 
342» .617 
343f .617 
344 » .617 
345r ,617 
346* .617 
347» .617 
348» .617 
349r .617 
350» ,617 
351r .617 
35?» ,617 
353 f ,617 
354f .617 
355 » , 617 
356» .617 
357» .617 
358> .617 
359» .617 
360 » .617 
361 ,. 61 7 
362r .617 
363*^A17 
364» .617 
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TABLE C-4. (Continued) 

‘tAHlTh NODL OAMPiNb 

(jAMHO » • i 
l'iL0OPf41 SA 

AKLXCA'.lJl “bKhLAXATiON CRlTtPiA 
UKLXCA ► .01 
uTIMFI • .UUbi^5 
CSGFAC*10. 

1 IMt>llJ »4 . 0 

OUTHHT» 0.1 

END 

bC[) OARHAY UATA 

1 S TE''^P qUIIU) uE wax RcF -bO-F 

(•. »bn. f .f^49iir 73. f . lo7<:?f73. » .2U3b»90. f .bV52» yo . 

.7206»lbU. »hNlj ^ REK NOuE v/ERSUS TEMP 
-2 % LABEL ARRAY FOR WAX 
G-WAX »DO*TlRlr Q*TlMf I WAX »ENU 
-3 i LABEL I EMPb EXCEPT WAX 
T1 fT2Fl3f TH»Tb»T6» l7fTb»T9»Tin»Tll 
Tl^: »T13. Tl4f T) bf Tlo» T17» Tltt» Tl9f T2U»T21 » T22 
T2o»T24»T2bf 1 2t> » T2 / » T?B » T29 » T30»T31»END 
4r0. .lOU. » .139 » 94.0 » .27B»9l.6» ,417 » 91. 3 r ,bbb»91.b» . 69b » 9 1.0 
. 834 , 9U . 6 » . 97 3 f 90 . 2 » 1 . 1 1 2 ► 88 . b » 1 . 2 d 1 f 8b . b » 1 . 390 » 84 . 3 » 1 . 529 
82.1 »1.8 74»80,U» I.f5l3»78.8»l. 94b »77,2»2.u8b»75.0»2.224»72.8 
2 . 3b3 » 7 1 . 6 » 2 . bU2 » 7(i . 0 » EMO 
REM MIPHLE FIr NODE - 1 EST TEMP T?.6 

b»0. » lOu. Of . 139»96.4» .278»-9o.2r .417 » 91 .7 » .bbb»91 ,0» . 695 » 90.6 
. 834, 90. 3 f .973 r 90.4. 1.112 f 90.3 f 1 .261 »90.bf 1.390 » 90.4 f 1 .529 
90 .6. 1 .674 .90 ,U r 1 .nl3.84i.9» ] .946 , 89 , 4 » 2 . 08 b » 88 . 2 . 2 . 224 . 85 . 7 

2.363»B3.6.2.bU2»fil .3»EflU 
REM middle fin NUDE I EbT TEMP 127 

D,o. . lOU.Of .139»98.b» . 278 » 9b . 5 , . 41 7 » 93 . 2 . . bbb r 92 , 1 » . 695 » 91 . 0 
,834»9u.7..937.9l.l»1.112»90,9fl.2bl.91.3fl.390»91.2»l. 529 
91 .6.1 .674 .91.0. 1 .cl3.90 .8. 1 ,946»9u.b»2. U8b»90 .8. 2.224 .90.3 
2.363»9U.7»2.bu2.9u.0»Ef'D 
HEM TOP FXM NODE TEbT TE'-'P 28 

7.0. . 100.0* . 139.80 .6. .278.73.6, .417 .70.3. ,5bb»67 .7. . 695 .65.8 
.834,83.8. .973.62.9. 1. 1] 2.80.9. 1.281.59,4.1.390.58.4. 1.529 
57, 8. 1 . 6 74. bo, 4. 1.6 13. 55. 7. 1.946. 54, 7. 2. 085. 53. 6. 2. 224. 52. 3 
2,363.51 .3»2.5U2.50.4.tMU 
KEM PLATE NOuE TEMP TEbT DATA 133 
EDO 
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C Sou MR 
CNFwHK 

tr4D 
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TABLE C-4. 

BCD 5VAHTAHLPS i 

i')METfH(Tl»T201»G^01»Kl) 
GMETFK(T2rTiiUl»G301 »K2) 
GMF.TFK{T2Ub»T201»G2u'5»K^) 
'vMETFK (T2U2f T201 »G3 u2»K4) 

1 FK n 1 » T 2 0 2 » G 2 1 1 2 » K 5 ) 
4MFT FH ( T 20 1 » T202 » G302 » Ko ) 
wNpFlFF(T20b»T202»G2bh»K7) 
(vMFT FR (1 2 0 3 » T2 U2 » <;3u 3 » K h ) 
GMFTFR(T1 »T^;n3»b2n3rKy) 
GMETFR(T202»T20o^G3u3»Kin) 
ivi^Fl f K ( T2f)7 f T203 rG2u 7 f k 1 1 ) 
GMFT(-K(T2(J4»T20o»G3<i6»K12) 
u^1E^^:K(Tl »T2U4»G204»K13) 
4MFlFKn2()3rT204»G3iibfK14i 
(vMFTFKn20B»T204»G2b4»KlS) 
OHFI FR ( T3 » T2(»4 r G307 » K lb ) 
‘^FETFR n^iOl f r20b»G2ub»Ki7) 
OMFTFR n4fT20b»Go0H»K18) 
(.vFI- I F R ( 1 2 ()9 r T 2 U b » G2 U 9 > K 1 B ) 
' O''! F I F H ( 1 2 0 6 r T 2 U b » 3 0 b » K 2 n ) 
I F R ( 1 2U2 ► Tii Ub » C’2 1»6 f K 2 1 ) 
tyFiE lFRnPObf T2 Ud»G3o 9 rK22) 
f.'KF TFR ( 1 2 1 0 f r 2(jt) f G2 1 0 f K23 ) 
uFF TFR n 2 07 f T 2 Oo » G3 i 0 f K24 ) 

‘v!«.FTF R ( T 2 03 » T 2 0 / . (,,2 u 7 » K 2 G ) 
CN;E I F R ( T20b » 120 / » G3jl 0 » K2b ) 
<>MFTFRn211 » T2i; /rb2il »K27) 
FRFiFRn208. 120 /»G3i 1 r k^M} 
t-W'FlrRn204»T20b»G2b«»K2b) 
p'j F F i F R ( 1 2 IJ 7 » 1 2 ( i fS r 3 i 1 » k 3 0 ) 
oFElF R(T212» 12db » G2 L2 f K3 1 ) 

< >HF If R { T b » T2 Ub » G3 1 2 » K 32 ) 
NR.F rFR(T20b»1209 rG2i;bfK33) 
•/F'F;T F R (Ifif 1209 » b313 f K34 ) 
GFE I F R (1 2 I 3 » T209 r G? 1 3 » K3‘> ) 
ijFFTFR ( 12 1 0 » I 209 » G3 1 4 » K 3^,. ) 
i^MF: I F R ( T20h » T2lU r G2 1 0 » KbV i 
' . F F. 1 F R { 1 2 0 9 » T 2 1 U » G 3 X 4 f K 3 ' ' ) 
Gf^'F IFR ( T 2 14 r T2 1 U r G2 i Li » K 3 ' ) 
gF F I FR ( T2 1 1 » T21 U r G31 b r KUf, ) 
Giv'iETFR (T207» 121 i rG2 H f k4] ) 
OMF TFR n 2 1 0 » T2 1 i » G3 I S » K 42 ) 
GMF lER n 2 1 b » 1 2 1 1 r G2 i b f h 4 3 ) 
GMFlFR(T212f 1211 »G3ib»K44) 


(Continued) 

0METFR( r206k 1212rG2i2»K4b) 
bi''iF TFR(T211»T2l2»G31b»K4b) 
tvF!ElFRn21b» 1212»G2lb»K47) 
'.-iN L 1 FK 1 7 7 » T2 12 » G 3 1 7 » K 4 f ) 
G'MFTFR(T2U9»T2l3f G2x3» K49) 
G'NpFTFR (TbrT2l3rG316f kbO ) 
0NpFTFK(T217»T213rG21 7»Kbl ) 
G ML 1 F R n 2 1 4 r 12 1 3 » G 3 1 9 » K b2 J 
GF F 1 FR ( T2 1 U » T2 14 ».g2 x 4 » k :j3 ) 
G F, F T F R (1 2 1 3 » T 2 1 4 » G 3 i b , K b 4 ) 
^v^ t If R ( T21H » I 214 r G2 J H » K bb ) 
uF-' £1 FR ( T 2 1 b » T2 1 4 f G32 n f K bb ) 
<JMFTFH ( T 2 1 1 » 12 1 3 » G2 I b » k r.7 ) 

< ■££! F H n 2 1 4 r T 2 1 b hG 32 It » K 30 ) 
x'k F 1 FR ( I 2 1 9 » 12 1 3 r G2 X 9 r k 3b ) 
^■MtrFR(T2lDrl21b»G32l»KF'U) 
(.■MF7rR{T2l2»T21brG2lb»Kbl ) 
GFF7F R (T21b»T21o»G32l »Kb2 ) 
v»V.LlFR(T22U*721u»G2y-0»Ke,3) 
GMF lFR (T9»Ty:;lb»b322»K.b4) 

G4 F 1FR(1213»T217» G2 1 7 » K b 3 ) 
gFiF.T FR ( T 1 0 » I 2 1 7 » ii32b f kbn ) 
G-ME 1 EK ( T22 1 » T2 1 7 » G2-- 1 » Kb 7 ) 
gFETE I' (T218»I217» G3yr 4 , Kb.'. ) 
•iNElFRn214»T21b»G2j GfKbb) 

''WETFR n 2 17 » T2 lb » G3.r:4 f k /.) ) 
GMFTFR(l222f T21<irG222»K71 ) 
(>METFR ( 1 2 1 9 > T2 1 b » G3.' G r K 7 2 ) 
x/FFl EP n 2 lb » T219 r G2 xb , K 73 ) 
gMF;TF R ( 12 1 8 ► 12 1 V » G3b b » K 7^ ) 

< kOi lFR ( 1223 r 1219 »G2G3»K YbJ 
i.'iVFl F (7 { 1220 f T21 9 » G3bb f K 7b ) 
<‘MF fkR{T216»7220»G2y^i)»K77) 
GFE l FR { 72 19 » T220 r G32b f K 74 ) 
<.MF7FR<T224» 122u»G2.-4rK7b) 

< ;|v.F7f R ( 7 1 1 r 1 220 H.32 / » Kbli ) 

( . F F 1 F R i T 2 1 7 » 7 2 2 1 ► G 2 y 1 » K 1 ) 
■..''''FTFK(712» 1221 rb32-»K82) 
Gk[ 7 F R ( T22b» 7 22l f G2< b r kb “i ) 

' iMFl FR ( 7 222 » 1221 t G32b r Kb4 ) 
'.■FF IfR ( 7218 » 1222 * G2.^2 » Kbb ) 
|>F1FR(7221 ► 1222»G3b9»KFb) 

' k'l If H ( f 22G r f 222 » G2, b » K t>7 ) 
Gf.it iFR (7223F1222»G3onfK84) 
GMF7FR(T2l9»T223»G223»kb9) 
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(. Mfl fR (1 

('NiFTfc H n ?-?4 » T22:i » G3.il » K9? ) 

( MEl FRll^'-^O# r^24»G2^4f K93) 

( 4/.F 1 F R n 223 r T224 » G32l » K94 ) 
^•.^.F^FR^22n»T2?4»G22H^K9‘)) 

'vNiF 1 F R ( T 1 3 FI 224 f G 3.32 » K96 ) 
wN'F TF R (12?! »T22D»G225rK97 ) 
uNFlF kU 14f t22b»G333»K9h) 

OF F I F K n 229 » 1223 » G2<i 9 » K99 ) 

i.KEl P R t T 22b t T22b t G3.^4 » K 1 G U ) 
uFif- l FK ( I 222 » T22o » G22b » K 1 G1 ) 
^,,F;F1FR (T22b» T22o»G3o4#K102) 
fjl^F jf. ( I 231 ) , T22t)» G23n f K1D3) 

uMt 1 FR n 2 27 ► T22o r G3 » K 1 04 ) 

< > HF_ I F K n 2 2 3 f V 2 2 7 » G2 2 7 » K 1 n D ) 

ONif I F K ( T22b » T 22 V » G32b » K 1 06 ) 

u>N.fc 1 F K ( I 2 3 1 » T 2 2 7 r G 2 1 f K 1 n / ) 

GKF.TFH I T22ri » T227 » G3.ib » K 1 Ott J 
.iMFT F K 1 1224 » T22d » G22R » K1 09 ) 
wRFTFR < 1 227 » r228 » G30b » K 11 U ) 
i'F t TFR 1 1232. » T22b » G?.i2 » K 1 1 1 ^ 
(-(N.E i F K ( T 1 "3 » 1228 » '^33 i » K 1 123 

,N F IF H ( T 22 h » T 229 f G2 c9 » K 1 1 2 ) 

UNiF 1 F H n I -3 » I 229 f G330 » K 1 I '* ) 
UMFll FK ( T233 » 1229 » G233 » K 1 1 b ) 
wMF" T FK n 230 » 1229 » G339 » K 1 1 D ) 
CjMF T F H ( T22b » T 23U » G2b0 » K 1 1 / ^ 
WMF IFR IT229 » 123U t G3 J9 r K 1 I h ) 
wN’FlF Hi [254»T23U»G234fKll9) 
GR.ElFH ( 1 231 » 123 u » G340 »K 12U ) 
oNT'l F K ( 1 227 r 1231 » G23 1 » K 1 2 1 ) 

O'l^FlFK (1230 » 1231 * G34 0 r K 1 22 ) 
lFH(T?3b»1231*G2ob»K123) 

OMF I F K ( 1232 » 1231 » G34 I » K1 24 ) 

< N F IF R (T226» 1232 » 02.32 »F 12b) 
GMIIFR (1231 » T232»G34l »K 12t>) 
uFF 1 FR ( 1 23b » T232 » 023b » K 12 7 ) 

ivNiFlFHll 1 7»1232*G342»K12A) 

* jNiFI FF((1229»r233»G233»Ki29) 

(.MFl FK (1 13 » 1233 » b343 f K130 ) 
ol FlF'R ( 1 2 37 » T 23o * G2.37 » K 1 31 ) 

, ,|v'F lFR(1234»12.33»G.344»Ki32) 

1 F R ( 1 2 3 n » 12 34 » G234 f K 1 33 ) 
GF''t' TFR ( 1233 » 1234 » G344 » K 134 ) 


(Continued) 

«MElFR(123H»T234»G2J8,»K13b) 

UMFI FH (1 2 3b r 1234 . G34b . K 1 36 ) 
uf>p1fR( f231» 123b»G2.3b»K137) 
GMFTFR ( 1234 rl 23b »G3nb#K 13b) 
wNiFTFR ( 1239 » 123b » G? .39 r K 1 39 ) 
(yf.'Fl FR ( 1 236 » T23b » G34b » K 14 0 ) 
wMF If R ( 1232 » 123o » G236 » K ) 

,.(yiF I F R (T23brl23orG3'+b»K142 ) 
^m-^FTF K 1 1 240 f 1 23o » G24(l » K 143 ) 
^.FiF T FR ( T 1 9 » 1236 » G34 / * K 144 ) 
URFI FR ( 1 233 f 123 7 r G237 » K. 14b ) 
wRF: TFR (1 20 » 1 237 * b34M * K14G ) 

,y(,,f.lP>in24lfl237»G24l»K147) 

GMFTF (T2.iB 1 1237 » G349 » K 14« ) 

. iN'FlFR ( 1 234 » 1236 » G2 o8 » K 1 49 ) 

>vMF 1 1 K ( 1 2 3 7 » 1 2 .36 » 034 9 » K 1 ) ) 

OMt 1 F R ( 1242 » 123ti » G242 » K .1 SI ) 
UI^Fll F'R ( 1 2.^9 » 1236 » G3b0 » K 1 b2 ) 
(JMF.l FK ( 1 23b » 1239 * G2-»9 » K J ‘i3 ) 

( -MF 1 F R ( 1 2 36 f T239 » G3‘J n » K 1 b4 ) 
WMF! FR (1243 » 1239 *G?43»K]bb) 
wf^FlF K( 1240 » 1239 rG3sl rKlSo) 

j p K ( 1236 » 1240 r G24 0 » K 1S7 ) 
6F:F TF"R (1239fl24U»f^3bl»Klb6) 
vMFIFR ( 1244 » I 24 0 » G244 f K iS9 ) 
i.'MFl FK (I 2 l » 1 240 » G3b<r: » MbO ) 
(.Ft TF K (1237 » 1241 »G24l »Klbl ) 

i.MF 1 FK ( 12.2 » l'24 1 » 0.353 » K 162 ) 
t MFl F K ( 1 24b » 1 24 1 » 024 5 » K 1 03 ) 
(-jVF'IFH (1 242 » 124 1 * (^3t)4 r K 164 ) 

( ;FiF.lFR ( 1236 » T242 » <724 2 t K 10b ) 
oFh 1 F R ( 1 24 1 » 1242 » G3s4 » K 1 Ou ) 
(.RE IF R ( 1 246 f 1 242 » G?hO » K 1 07 ) 
GME IFR ( 1 243 » 1 242 » 03bb » K 1 n6 ) 

•jMF I F R ( 1 2 39 » 1243 » 024 .5 » K 1 ) 

ORE IF R ( 1242 » 1 243 » G3bb » K 1 7 0 ) 
(7RFlFR(1247»l243»G247fK ) 71) 
>iMF I F R ( 1 244 » 1 243 » G3b6 » K 1 72 ) 
uMF IFR (1240 t T244 r G2 h4 » K J 7o ) 
6NiFl FR ( 1 243 » 1244 » G3bb » K 1 74 ) 
GRfclFK(124Hrl?44M,246fK17b) 
(.■mpl P R I r;^3» V244 »G3b 7 » K170 ) 
OFF 1 F R ( 1 4 1 r T 2 h b r G2 H 5 » K 1 7 7 ) 

6MFTFK ( 124 * 1 24b » 1 74 ) 
GMF IFK ( 1 249 1 1 24b » G249 t K 179 ) 
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uMETFK(T?4b»T24b»G3:39fKl Au) 
<wMETFRn2‘+2»T24o»G2‘+6»KlAi ) 
lJ^ltTFK(T24^J^,T24t^»G3b9^Kl^£;) 
OMbTFRn2bU»T24b»G2oO»KlA3) 
WMF! FR(T247fT24t»»G3o0»K18H) 
V/MFlFRn243rT24 7rG247.K18b) 
OVFTFR(T246»T247fG3t3n»Kl Pfa) 

WMF TFRn2bl »T247 »G?d1 »K1R7) 
GMFTP?K ( T248 » T24 7 f G3o I » K 186 ) 
GWE7 FH n 244 » T24 tit 6246 »KJ h 9) 
CvNiETFKn24 7rT246t636ltKlMfj) 
Gi^FTFR ( 1 2b2 1 T246 1 G2b2 1 K 191 ) 
GI«FTFR(T2btl248tG362tK192) 
<vMtTFR(T24bt 1249 1 6249 1 K 193 ) 
^N^ETFR (T26t 1249 tG36btK 194 ) 
uMElFRn 23t T249t02b3tK19S) 
6METF R(T2bU»T249t G3o4 t K 1 9b ) 
ttMElFR(T246t T2b0tG2b0tK197) 
RN’FTF R ( T 249 1 T25U 1 63b4 r K 1 9h ) 
GMKl FR ( T 30 t I 2b0 f 6254 t K 1 99 ) 

OMElFR(I2bl t l2bUtG3D5tK20U) 
9N-ETFR(T247tT2bltG2bl tK2nl) 
wMETpkdPbOt 12bl tG3b5tK2(t2) 
9VETFR(T31 tl25l tG25btK203) 
GMElFR ( 1252 1 1251 1 G3b6 1 K204 ) 
wr^ETFRn248tT252tG2b2»K2n5) 
6METFR(T251 t 1252tG3bGtK2nb) 
OMETFK(T29t T2b2tG25btK2U7) 
«METFR(T27t I 2b2t(;3h7 tK2U6) 
AL)D(Kl»K2tK3tK4tK209) 
AD0tK5tK6tK/tKBtK21b) 
AL)O(K9tKl0tKlitK12tK211 ) 

ADn(Kl3»H4tK15tKl6tK212) 

ADD(K 17rK lHtKl9tK2n tK213) 
mDD(K 21 tK22tK23»K24tK214) 
ADDlK25tK2btK27tK28tK21b) 

ADD 1 K 29 r K 3U t K 31 1 k 32 t K 2 lb ) 

ADD ( K 33 » K 34 t K,35 1 K 36 1 K2 1 7 ) 
ADD(K37tK38tK39tK4n tK2lH) 
ADD(K41f K42tK43tK.44tK2l9) 
ADD(K45tK46tK47tK48tK220 ) 
ADD(K49»K50tK5l tK52tK221 ) 

ADD ( K53 » K 54 t Kbb t K 56 tK222) 
ADD(K57»K56tK59tK60 tK223) 

ADD ( K61 » K62 1 K63 1 K64 1 K224 ) 
ADD(K65tk66tK67tK68tK225) 


(Continued) 

ADD (Kb9tK7UtK7l»K72tK22b) 
ADD(K73t|<74tK75tK76t K227) 

ADD ( K 77 1 K 78 1 K 79 1 K 80 t K 226 ) 
aDI)(k« 1 tKH2tK83tK84tK229) 
ADD(KH5tK6btKH7tK88 t K23U ) 
ADD(K89tK90f K91 tK92 tK231 ) 
ADD(K93tK94rK9btK9b»K232 ) 

ADD ( K 97 r K9H t K99 1 K 1 0 0 t K23 ^ ) 
ADIKkIOI tKlU2tKlU3tr,104rK234) 
A0D(K10brKliit>tK107th 108tK23j) 
ADD(Kl09tKll0tKll 1 tM12tK23o) 
ADD (K113tK114tK115t^llfo,^23/) 
AUD(K117tKll«tKll9tKl2n tK23b) 
ADD(k 121 tK 122 tK 123 t^ 124 tK239) 
ADD(K12btKl2btK127t(<. 128rK24u) 
At)D(‘tl29tKl3i)tK131 t»^ 1 3?tK24 I ) 
ADD(K133tKlo4»K135t,^13btK242) 
ADO (k 1 37 1 K138 t Kl39t M40 t K243 ) 
AUD(K141tKl42tK^43t^l44rK244 ) 
AnD(K14b»Kl4btKl47trvl48tK24-j) 
ADD(Kl49tKlb0tKial , 1 b2rK24b) 

ADD(K153tKlb4tKl55tr. lb6tK24 /) 
ADDtKlb /rKlbBtKlbPt,^ loD tK24b) 
AD0(K161 tKlb2tK lb3tMb4f K249 ) 
ADD(Klb5tKlb6tKlb7tt\lb8rK2bo) 
ADD ( K 169 1 Kl /O t K 17 1 1 t.l72 1 K25 1 ) 
A0D(K173*K1 74tK175tM7btK252) 
ADD (K 177 tK 178 1 K179 1 MrtO tK2bo) 
ADD(KlttltKl82»Klb5tr.l84tK254) 
ADDiK185tKla6tKl67tM8btk25b) 
ADU(K189tKl90tK191 tKl92 tK25b) 
ADD (Kl9j»K194»K195r^l96tK25/) 
ADDtKl97tM98tK199t^2U0 tK2b6) 
ADO i K201 1 K2U2 » K203 t»s 204, K259) 
ADD(K2UbtK2U6»K207tt\2U8,K26o) 
DlDEGl (K313tAltl201 ) 

DlDEGl (K3l4t A1 , 1202) 
l.JDtGl (K315t A1 1 1203) 

DlDEGl (K316t Alt 1204) 

UlDEGl (K 31 7 1 Alt 1205) 

DlDEGl ( K3 18 tAltTEOb) 

DlDEGl (K319tAl t 1207) 
ltlDE6l<K320t AltT208) 

I'lPEGl (K321 1 A1 1 1209) 

OlDEGl (K322t Alt 1210) 

DlDEGl (E323tAl tT211 ) 
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DlDtel f Al » 12 ) 

LintGl (K32br Al »1 213) 

I IDtGl (K326» Al »T214) 

(K327 t M » T21b) 

L lDt-ftl (K32H» Al » 1216) 

LUDbGl (k 329» Al f T217 ) 

I IDtGl (K33(I»A1»T21H) 

L/lOtGl (k 331»A1»T219) 
l..ll)tGl(K332» Al» T22n) 

GlOLGl (K333» Al r 1221 ) 
niDh-Gl (K334»A1» 1222) 

[JlDIr Gl (K33b » Al » 1223) 
L>1 DLg1(k336»a1»T224) 

UlDLGl (K337»AlrT22b) 

L'lDt Gl (k 338» Al » 1226) 

[ IDLf.l (k 339» Air 1227) 
ulDEGl lK340r A1*12?B) 

DlDfcGl (K341 » Al fT229) 

DlDbGl ( k 342» Al » 1230 ) 

I'l DtGl (K343» Al » 1231 ) 

UlDtGl (K344r Air 1232) 

DlDtGl (K34br Al »T233) 

UlDtiGl 1 k 346 r Al r 1234 ) 

UlDtGl (K347rAlrT235) 

OlDtGl (K348r Al rl236) 

OlOLGl (K349 r Al r 1237 ) 

UlOLGl (K3b0 r Al r 1238) 

L'lDtGl (K3bl r Al r T239) 

UlDtGl {K352rAlr 1240) 

UlUtGl (K353r Al r T241 ) 
DlDtGl{K3b4r Air 1242) 

UlDtGl (K355rAlrl243) 
UlDtGl(K356rAlrT244) 

UlDtGl (K357rAl r I 245) 

UlDtGl (K3b8r Al rT246) 

UlDtGl (K359r Air 1247) 
l>lDtGl (K360 r Al r 1248) 
DlDtGl(K36lrAlrT249) 
DlDtGl(K362r Alrl2b0) 

UlDtGl (K363r Al r 1251 ) 

UlDLGl (K364r Air T252) 
NiLTPL.Y(K209rDTlMEUrK261) 
MLTPLY ( K2 10 r DT IMtU Ms2b2 ) 
MLTPLY(K211»DT IMEU»t263) 
MLTPLY ( K2 12 r DT I MtU r K264 ) 
MLTPLY(K213rDTlMEUrK265) 


(Continued) 

MLTPLY (K214rOTlMtUrK26b) 
MLTPLY (K215rD^MEtJr^267) 

MLTPl.Y (K216rDTlMtUrK26P ) 
MLTPLY (K2l7rUT IMhUr is269 ) 
ML TPL Y ( K2 1 8 r 1) 1 IMhUrK270) 
MLTPLY (K2 19 r dT IM tU r k271 ) 
MLTPL Y ( k 2?0 r dT IML l I r k 272 ) 
MLTPLY(K221rDTlMEUrK273) 
MLTPLY(K222rDT I,vltUr^274) 
MLTPLY (h?23rDTlMtUrK27b) 
MLTPLY (k 224 r IjllMLU r K?.7b ) 
MLTPLY <K225 ml Iivlt llrK2 77) 
MLTPLY (K226rf)T I MhU r r27b ) 
MLTPLY ( K227 r QT I MEU r K279 ) 
MLTPLY (K228rDTlMtUr ^2B0 ) 
MLTPLY(K229rDTIMEUr^281 ) 
MLTPLY (K230rUT I MtU rrv282) 
MLTPLY(K231rDTlMEUrK283) 
MLTPLY (K232rDTlMEUr (-,284) 

MLTPLY ( K233 r UT I MEU r K 285 ) 
MLTPLY (K234 r UT I MEU r -v?8b ) 
MLTPLY (K235rDl IMtUris?87) 
MLTPLY (K236rDT IMEUrK288) 
MLTPLY lK237rUT lMtUr^2«9) 
MLTPLY(K238rDl lMFUrK290 ) 
MLTPLY(K239rUTlMtUrK29l ) 
MLTPLY ( K240 r DT IMt ( I r k 292 ) 
MLTPLY(K241 rOTlMEUr ^293) 
MLTPLY ( K242 r QT IMEU r K294 ) 
MLTPLY (K243rDTlMLUrK295) 
MLTPLY ( K244 r DT I MEU r k296 ) 
MLTPLY (K245rUT IMEUrK297) 
MLTPLY(K246rOT IMEUrK298) 
MLTPLY (K247rOT IMEU rK299) 
MLTPLY ( K249 r Dl IMhU r K30 1 ) 
MLTPLY (K248rOT IMLUrK300) 
MLTPLY ( K250 r OT I ML U r k30? ) 
MLTPL Y(K251rOTlMtUrK3U3) 
MLTPLY (K252rDTlM£UrK3U4) 
MLTPLY ( K253rDT I MtU rK305) 
MLTPLY (K2b4»DTlMLUrK306) 
MLTPLY ( R255rDTl MEU rK3U7) 
MLTPLY (K256»DT IMEU rK308) 
MLTPLY (K257»DT IMEU rK309) 
MLTPLY (K258»DTlMEUrK310) 
MLTPLY (K259rDTlMEUrK311) 
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TABLE C-4 


(Continued) 


MLTPLY ( K260 p DT iMhU » ^312 ) 
ADD(K?fol »K313 »kj13) 
/\n 0 (K 26 £^»K 314 »K 314 ) 
ADri(K2ti3»K315»K3lS) 
ADD(K264f K3l6rK31h) 
AL>D(K2bbrK3i 7 »K317) 
Ar)D(K26t>f K3iy »K314) 
A[lDiK2b7»K319»K3i9) 
AUri(K?68fK320»K32 0) 
At)D(K269»K3^1 »K321 ) 

ADD (K27U rK3^i2 » K322 ) 
AUP(K271 ►K3<i3»K323) 
Af'D(K272»K324»K024) 
aDD(K273»K3^;5»K32‘S) 
ADD<K274»K326»K326) 
ADD(h27b»i<3^7»K327) 
ADD(K27b»K328»K32M) 
A[iD(K277»K329»K329) 
ADU<K278fK330»KO30) 
A0U(K279»K3bl »H33l ) 

ADD ( K280 ► K 332 »K 33?) 
An0(K281»K3o3»K333) 
ADL)(K282»K334»K334) 

ADO (K283»K3o 5»K338) 
A0D(K284»K33b»K338) 
AUD(K28b»K337»K337) 

ADD ( K28b » K 338 » K338 ) 
A[)D(K287rK339»K339) 
aDD(K288»K34(J p K34 0 ) 
ADD(K2B9»K341 »K34l ) 
ADDtK290 »k342 »K34? ) 


AD0(K2'^1 »K343»K343) 
Af)D(K292»K344 »K3 a 4) 
Alni ( k293 p K3nb » K 34b ) 
ADD(K294»K34 o»K 346) 
ADn(K29b»K347»K3'4 7) 
AUD(K29D»K34d»K348) 
AUD(K297»K349»K3AQ) 
ADD(K298rK3bO»K3bn ) 
ADD(K299»K3bl »K3bl ) 
AnO(K300»K332»K3b?) 
ADD (K 301 »K3o3»K3b3) 
Ai)D(K302rK3b4»Kob4) 
ADD ( K 303 » »<3bb » K3bb ) 
ADD(k304f K33b»|< 3bb) 
ADD<K30b»K3r>7»K3b7) 
ADD ( K 30t> » K3bfl f K 3b8 ) 
ADD(K307»K3b9»K3b9) 
ADD(K308rK3bU»K3b(J ) 
ADD(K309»K3 o 1»K3M ) 
ADD ( K 31 0 » K 3o2 ' K 3c? ) 
>\DD ( K31 1 » K3o3 ► K3c3 ) 
AD0(K312»K3 o4^K3c4) 

D I Dt.G I { T I D » A 4 » r b ) 
i-lDLGl (TlviEDf A4» 17) 
lilDLGl ( TlWEWf Abf I 14 J 
l-lntGl (TI^ENr ASr r lb) 
UlDt G 1 ( T I MEIj » Ab » 1 lb) 
DlDKGl (TI7tD»Ab» H9> 
DlDtGl n iMEixr A7f 11 ) 

END 
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TABLE C-4. (Continued) 


BCD 3VAKIAHLES 2 
END 

BCD 30UTPUT CALlS 


STMijRD 

PRINTL1a2»K209»k 261 »K313» T201 ) 
PRINTLt A2»K210 »k 262»K314» T2U2) 
PRIMTH A2»K211 »K263»K31b» T2U3) 
HRINTL(A2»K212»K264»K316»T2U4) 
PR1NTL(A2»K213 »k 265»K317» T20b) 
PRINTL( A2»K214»K266»K318# T2U6) 
PRINTL C A2»K21b»K267»K319» T207) 
PRINTL(A2»K216»K26fl»K320»T2u8) 
PRIIMTL(A2»K2l7rK269»K321» r2UP) 
PRINTL( A2f K218»K270 »K322» T210) 
PRINTL(A2»K2l9»K271»K323r r2ll> 
PRINTLi A2»K220»K272»K324»T212) 
PRINTL( A2»K22l»K273»K325»T213) 
PRINTL ( A2»K222»K274rK326»T2i4) 
PRINTL(A2»K223»K275»K327»T21b) 
PRINTL(A2»K224.f K276»K328» T218) 
PRINTL(A2»K225rK277»K329f F217) 
PRINTL(A2»K22b»K27fl»K330»T2lH) 
PRINTL(A2»K227 »k 279»K331 »T219) 
PRINTL( A2»K228»K2fl0»K332»T220» 
PRINTL(A2»K229 »k 281 »K333»T22D 
PRINTL ( A2 » K230 » K282 » K334 » T222 ) 
PRINTH A2»K231»K283»K335»T223) 
PRINTL { A2 » K232 » K284 » K33b » T224 ) 
PRINTL ( A2 » K233 » K285 » K337 » f 22b ) 
PRINTL ( A2 » K234 » K286 » K338 » T226 ) 
PRINTL(A2»K238 »k 287»K339* T227) 
PRINTL(A2rK236»K288»K340f i2^8) 
PRINTH A2»K237»K289»K341» I22B) 
PR I NTL ( A2 * K238 » k 290 » K342 » T23d ) 
PRINTL(A2»K239fK291»K343» 1231 1 
PRINTL ( A2 » K240 » K292 » K344 » T 232 ) 
PRINTLi A2»K241 »K293»K34b»T233) 
PRINTH A2»K242rK294»K346r 1234) 
PRINTL( A2»K243»K29h»K347» 1235) 
PRINTL (A2»K244 »k 296»K34H» T23b) 
PRINTL(A2»K24b»K297»K349»T237) 
PRINTL(A2»K24b»K298»K3bO» T23«) 
PRINTL(A2»K247»K299»K3bl» T239) 
PRINTL (A2»K248»K300»K352» 1240) 
PRINTL(A2»K249»K301»K353»T241) 
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TABLE C-4. (Concluded) 
f-'KINTL{ Ai?«K25U»K:jO?»K Jb4» T?4?J 

PKlNTL(A?»Ki?5l »K.^03»K35b» ( 243) 

PKINTL ( A2»K2b2f K304»K3!b6f F2H4) 

PRINTL(A2»K2b3»K3n5»K3b7. T24b) 

KRIf-iTL (A2fK2b4,K306»K3ba» T246) 

HRINTL( A2rK25b»K307»K3b9» 1247) 

PR I NIL ( A2 f Ki^bR r k 3(18 » K36(! f T24B) 

(•'RK'NiTL ( a 2 » K2b7 » k,3 09 » K381 » I 249) 

PRINTL(A2»K25H»K3l0f K36k» (200) 

HRlHTL(A2»K2b9»K3H fK363» F2bl ) 

' RRIi')TL{A2»K2bU»K312’»K364f T2b2) 

PRINTL( AbfTi »Tk;» T3» I 4 » Tb » Tb » T / » TH » T9 » T 1(J r 11 1 r T12 » T 1 3 » T14 
Tib»T16»Ti7f TlHf Tlb,T20»721 »T22»T23»T24»T25»T26^T27 
T28»T29»T3 (i»T31) 

tMU 
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APPENDIX D 


TEST DATA 
Introduction 

Test data for twenty phase change tests run at Marshall Space Flight 
Center's Test Division are included herein. These data includes temperatures 
and phase chaise position time histories. Data for 1.9— cm (3/4-in.) and 
0.635-cm (l/4-in. ) cells for twenty tests are given. Of these, nine tests are 
melting runs and eleven tests are freezing runs. Data are divided into four 
sections; Section I, freeze temperatures; Section II, freeze front position 
height; Section III, melt temperatures; and Section IV, melt front position 
height. Figure D-1 gives a presentation of temperature instrumentation loca- 
tion for all runs. Only data from the last 1.9-cm (3/4-in.) and the center 
0.635-cm (l/4-in. ) cells are presented. 

General Data Description 

Temperature plots for test 230-5 through 230-15 were hand plotted, and 
all subsequent tests were computer plotted. Although numerous additional 
temperature data were acquired, only the fin and cell temperatures for the 
center 0.635-cm (l/4-in. ) cell and the last 1. 9-cm (3/4-in. ) cell(i.e., at 
fluid outlet) are presented. Thermocouple T31 through T33 represent fluid 
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temperatures not used in analyses and T-34 represents the lower plate temper- 
ature. The phase change position data were visually read from filmed data 
using a film analyzer. Readings were recorded and stored on computer cards. 
These cards were input to a 4020 computer plotter to generate position versus 
time plots. Position data for center of cell readings are given for all tests. 

In soi. e cases, the phase change position on the upstream fins are also given, 
these data are represented by starred and unstarred lines, respectively. In 
the case of the 0.635-cm ( l/4-in. ) cell, when only one set of data is given 
for both fin and center cell, this indicates that these two heights coincide. In 
some cases the data were truncated when the cells were completely melted, 
in others the data plotting was continued. In the later case, this can be easily 
seen by the flatting of the data plots. 

Across the top of some data plots the following information parameters 
are given; test type, internal test designation (official test designation), test 
duration, date the test was run, initial paraffin temperature and fluid inlet 
temperature. 

In test 230-52 only 2000 seconds of data is available because of 
improper camera positioning. In test 230-7 the 1.9-cm (3/4-in.) cell data 
were inadvertently lost. For this case the smoothed data are presented in 
lieu of raw data. 
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Figure D-1. Temperature instrumentation. 





SECTION I 


FREEZE TEMPERATURE DATA FOR 
FINNED THERMAL CAPACITORS 
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Figure D-2. Freeze temperature data for finned thermal capacitor (test no. 230-7) . 
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Figure D-5. Freeze temperature data for finned thermal capacitor (test no. 230-10) . 




Figure D-6. Freeze temperature data for finned thermal capacitor (test no. 230-11) . 
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Fio'ure D-8. Freeze temperature data for finned thermal capacitor (test no. 230-15) 




Figure D-9. Freeze temperature data for finned thermal capacitor (test no. 230-15) . 
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Figure D-10. Freeze temperature data for finned thermal capacitor (test no, 230-51) 
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Figure D-11. Freeze temperature data for finned thermal capacitor (test no. 220-51) 
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Figure D-12. Freeze temperature data for finned thermal capacitor (test no. 230-52) 
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Figure D-13. Freeze temperature data for finned thermal capacitor (test no. 230-52) 
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Figure D-14. Freeze temperature data for finned thermal capacitor (test no. 230 53) 
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Figure D-15. Freeze temperature data for finned thermal capacitor (test no. 230-53). 
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Figure D-17 . Freeze temperature data for finned thermal capacitor (test no. 230-54) . 
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Figure D 19. Freeze temperature data for finned thermal capacitor (test no. 230-55), 
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no. 230-56). 
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Figure D-23. Freeze temperature data for finned thermal capacitor (test no. 230-58) 




SECTION II 

FREEZE FRONT POSITION DATA FOR 
FINNED THERMAL CAPACITORS 
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Figure D-24. Freeze front position data (test no. 230-7) 
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Figure D-25. Freeze front position data (test no. 230-7) . 
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Call and Fin RaadiD9, in. 



Figure D-27. Freeze front position data (test no. 230-10) . 


232 











Figure D-29. 


Freeze front position 


data (test no. 230-11). 
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Figure D-30. Freeze front position data (test no. 230-15). 
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Cell and Fin Reading, 























Figure D-,'32. Freeze front position data (test no. 230-51) . 
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Cell and Fin Reading, in. 
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Figure D-33. Freeze front position data (test no. 230-51). 
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Cell and Fin Reading, in. 
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Figure D-34. Freeze front position data (test no. 230-52) . 
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Ceil and Fin Reading, in. 
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Figure D-35. Freeze front position data (test no. 230-52) . 
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Fioure D-37. Freeze front position data (test no. 230-57). 
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Figure D-38. Freeze front position data (test no. 230-54). 
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Call and Fin Raading. in. 



Figure D-41. Freeze front position data (test no. 230-55). 
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Cell and Fin Reading, in. 
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Figure D-42. Freeze front position data (test no. 230-56). 
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C«ll and Fin Reading, in. 
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Figure D-43. Freeze front position data (test no. 230-56). 
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Figure D-44. Freeze front position data (test no. 230-58). 
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Figure D-45. Freeze front position data (test no. 230-58) . 
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SECTION III 


MELT TEMPERATURE DATA FOR 
FINNED THERMAL CAPACITORS 
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Figure D-46. Melt temperature data for finned thermal capacitor (test no. 230-5). 
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Figure D-47. Melt temperature data for finned thermal capacitor (test no, 230-5), 
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Figure D-49. Melt temperature data for finned thermal capacitor (test no. 230-6) . 
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Figure D-50. Melt temperature data for firmed thermal capacitor (test no. 230-8) . 
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Figure D-51. Melt temperature data for finned thermal capacitor (test no. 230-8) . 
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Figure D-52. Melt temperature data for finned thermal capacitor (test no. 230-9) 



Figure D-53. Melt tcmperatvire data for firmed thermal capacitor (test no. 230-9). 
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Figure D-55. Melt temperature data for finned thermal capacitor (test no. 230-49) 
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data for finned thermal capacitor (test no. 230-59) 
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Figure D-63. Melt temperature data for finned thermal capaeitor (test no. 230-60) 




SECTION IV 


MELT FRONT POSITION DATA FOR 
FINNED THERMAL CAPACITORS 
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Call and Fin RaadinB.in. 
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Figure D-64. Melt front position data (test no. 230-5). 
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C«ll and Fin Reading 
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Figure D-65. Melt front position data (test no. 230-5) . 
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Cell and Fin Reading, in. 
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Figure D-66. Melt front position data (test no. 230-6). 
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Call and Fin Reading, in. 
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Figure D-67. Melt front position data (test no. 230-6) . 
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Call and Fin Raading. in. 
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Figure D-68. Melt front position data (test no. 230-8) . 
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Figure D-69. Melt front position data (test no. 230-8) . 
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Cell and Fin Reading, 
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Figure D-70. Melt front position data (test no. 230-9) . 
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Cell and Fin Reading, in. 
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Figure D-71. Melt front position data (test no. 230-9) . 
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Ceil and Fin Reading, 
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Figure D-72. Melt front position data (test no. 230-49) . 
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Cell and Fin Reading, in. 



Figure D-73. Melt front position data (test no. 230-49) . 
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Call and Fin Reading, i 
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Figure D-74. Melt front position data (test no. 230-50). 
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c«ll and Fin Raading,in. 


HtLT TEST 24 llJO-30) 


9733 SEC DUR *DaTE W19/7J ♦♦ llJf 


110112 
Oil Oil 



Figure D-75. Melt front position data (test no. 230-50) . 
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Cell and Fin Reading, in. 



Figure D-76. Melt front position data (test no. 230-57) . 
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Figure D-77. Melt front position data (test no. 230-57). 
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Ctll and Fin Roading, in. 
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Figure D-78. Melt front position data (test no. 230-59) . 
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Call ind Fin Raadinf, in. 



Figure D-79. Melt front position data (test no. 230-59). 
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Cell and Fin Reading, 



Figure D-80. Melt front position data (test no. 230-60) . 
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Figure D-81. Melt front position data (test no. 230-60) . 
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APPENDIX E 


DERIVATION OF INTERFACE BOUNDARY CONDITION* 

0S 

An observer on the interface at S(x', t) moves with the velocity — in 

3t 

9 S 

the y direction. He observes p— 6 x 6t mass change phase during the time 

9 1 

dx dx 

interval 6t in the distance interval x^ •“”5“ < x x^ + — . The heat 

dx 

liberated because of the change in phase is pAH ^ and is equal to the net heat 

dx 5x 

conducted away from the interface between x' - — and x^ + — during 6t. 

2 2 

In the y direction, the heat conducted away is given by 


- K 


L ay 


{- 


6 X 6t 


for 


y = s(x' ,t) 

In the X direction the heat conducted away is given by 


(E-l) 


9T 9T 

-K, - (-K, ® 


L 0x 


s ax 


) (-i '**) 


(E-2) 


for 


y = s(x- ,t) 

Equating the net heat conducted away from the interface and the heat librated 
due to phase change gives 


K. A. Rathjen; and L. M. Jiji: Heat Conduction with Melting or 
Freezing in a Corner. Trans, of ASME, Journal of Heat Transfer, February 
1971. 
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9T 


K 


S ^ 8x8t + K 


8T 


L as 


s ax ax 
ax 


L ax 

3T. 


ax 


• axat 


_ K — — 6x5t + K„ ^ 6x6t = pAH ^ 6t6x 

L ay s ay at 


for 


y = s(x' ,t) 

rearranging and dividing by 6x6t we have 


ax ax 

1y “ ~ ^ ^ 


as 

L ax ax “ s ax ax 


for 


y = S(x' ,t) 

We know the total differential is given by 

ax ax ax 

dx„ = dx + dy + dt 

s ax 8y at 

However, along the interface Xg = X^, a constant, so that 
^ aXg 

" ' "57 ^ ^ * "aT * • 

for 

y = S(x ,t) 


AX. 

pAH - , 


(E-3) 


and we also know 


, as , as 
'■5' = ta * at * 


(E-4) 


(E-5) 


for 


y = s(x,t) 
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Putting (E-5) into (E-4) gives 




a. 

ot ay at / 


y = s(x,t) 

Since x and t are independent 

ax 3T 
s ^ s as 

ax ■ ay ax ’ 


y = s(x,t) 


ax ax 
s s as 

at ” ay at 

from equation (E-6), for y = S(x ,t) 
Also by the same method 


as 

ax ay ax 


y = S(x,t) , 


L as . . , 

— at • y = s(!c,t) 


Substituting (E-7) and (E-9) into (E-3), we obtain 
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9T^ 9T^ 

K — — 

S 9y L 9y 





+ K = 

9y V9x/ 


Axi 9S 

pAH - 


for 


y = S(x,t) 


Since x' is a general coordinate 


9T 


K 


S 


9T, 


S 9y 


- K 


L 9y 


1 + 


= 5t 


for 


y = S(x,t) , 

If the dimensionless temperatures and are introduced, 

boundary condition form is obtained, 
where : 



T = Fusion temperature, 

F 

Tg = Solid Temperature, 

T = Liquid temperature, and 
L 

= Wall temperature. 


the desired 
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